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PREFACE 


This  report  describes  the  work  performed  under  Task  2  of  the  DOT/FAA 
High  Velocity  Jet  Noise  Source  Location  and  Reduction  program  (Contract 
DOT-OS- 30034) .  The  objectives  of  the  program  were: 

•  Investigation  of  the  aerodynamic  and  acoustic  mechanisms  of  various 
jet  noise  suppressors,  including  scaling  effects 

•  Analytical  and  experimental  studies  of  the  acoustic  source  distri¬ 
bution  in  such  suppressors,  including  identification  of  source 
location,  nature  and  strength  and  noise  reduction  potential 

•  Investigation  of  inflight  effects  on  the  aerodynamic  and  acoustic 
performance  of  these  suppressors. 

The  results  of  these  investigations  lead  to  the  preparation  of  a  design 
guide  report  for  predicting  the  overall  characteristics  of  suppressor  con¬ 
cepts  from  models  to  full-scale  static,  to  inflight  conditions,  as  well  as  a 
quantitative  and  qualitative  prediction  of  the  phenomena  involved. 

The  work  effort  in  this  program  was  organized  under  the  following  major 
Tasks,  each  of  which  is  reported  in  a  separate  final  report: 

Task  I  -  Activation  of  Facilities  and  Validation  of  Source  Location 
Techniques. 

Task  2  -  Theoretical  Developments  and  Basic  Experiments. 

Task  3  -  Experimental  Investigation  of  Suppression  Principles. 

Task  4  -  Development  and  Evaluation  of  Techniques  for  Inflight  Investi¬ 
gation. 

Task  5  -  Investigation  of  Inflight  Aero-Acoustic  Effects  on  Suppressed 
Exhausts. 

Task  6  -  Preparation  of  Noise  Abatement  Nozzle  Design  Guide  Report. 

Task  1  was  an  investigative  and  survey  effort  designed  to  identify 
acoustic  facilities  and  test  methods  best  suited  to  jet  noise  studies.  Task 
2  was  a  theoretical  effort  complemented  by  theory  verification  experiments 
which  extended  across  the  entire  contract  period  of  performance. 

This  volume  is  a  supplement  to  the  Task  2  report,  documenting  a  general 
method  for  predicting  the  aerodynamic  and  acoustic  behavior  of  turbulent 
jets.  The  objective  of  the  report  is  to  provide  users  with  a  description  of 
the  method  and  associated  computational  procedure  in  sufficient  detail  that 
it  can  be  implemented  and  utilized  as  a  useful  engineering  tool. 


Task  3  represented  a  substantial  contract  effort  to  gather  various  test 
data  on  a  wide  range  of  high  velocity  jet  nozzle  suppressors.  These  data, 
intended  to  help  identify  several  "optimum"  nozzles  for  inflight  testing 
under  Task  5,  provide  an  extensive  high  quality  data  bank  useful  to  prepara¬ 
tion  of  the  Task  6  design  guide  as  well  as  to  future  studies. 

Task  4  was  similar  to  Task  1,  except  that  it  dealt  with  the  specific 
test  facility  requirements,  measurement  techniques  and  analytical  methods 
necessary  to  evaluate  the  inflight  noise  characteristics  of  simple  and  com¬ 
plex  suppressor  nozzles.  This  effort  provided  the  capability  to  conduct  the 
"flight"  effects  test  programs  of  Task  5. 
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1.0  SUMMARY 


This  report  represents  a  supplemental  Task  2  effort  under  Contract 
DOT-OS-30034  documenting  the  computerized  jet  noise  prediction  method.  A 
complete  description  of  the  computer  program  is  provided,  including  examples 
of  input  preparation,  output  cases,  and  a  listing  of  the  FORTRAN  computer 
code.  The  mathematical  model  is  briefly  summarized  (it  appears  in  detail  in 
the  Task  2  report  proper) . 

A  unified  aerodynamic/acoustic  prediction  technique  has  been  developed 
for  assessing  the  noise  characteristics  of  suppressor  nozzles.  The  technique 
utilizes  an  extension  of  Reichardt's  method  to  provide  predictions  of  the  jet 
plume  flow  field  (velocity,  temperature  and  turbulence  intensity  distribu¬ 
tions).  The  turbulent  fluctuations  produced  in  the  mixing  regions  of  the  jet 
are  assumed  to  be  the  primary  source  of  noise  generation,  as  in  the  classical 
theories  of  jet  noise.  The  alteration  of  the  generated  noise  by  the  jet 
plume  itself  as  it  propagates  through  the  jet  to  the  far-field  observer 
(sound/flow  interaction  or  fluid  shielding)  is  modelled  utilizing  the  high- 
frequency  shielding  theory  based  on  Lilley's  equation. 

These  basic  modelling  elements  (flow  field  prediction  turbulent  mixing 
noise  generation,  and  sound/flow  interaction)  have  been  coupled  together  in 
a  discrete  volume-element  formulation.  The  jet  plume  is  divided  into  ele¬ 
mental  volumes,  each  roughly  the  size  of  a  representative  turbulence  cor¬ 
relation  volume  appropriate  to  that  particular  location  in  the  plume.  Each 
volume  element  is  assigned  its  own  characteristic  frequency,  spectrum,  and 
acoustic  intensity.  The  sound/flow  interaction  effects  for  each  volume  ele¬ 
ment  are  evaluated  from  the  flow  environment  of  the  element.  The  individual 
volume  elements  are  assumed  to  be  uncorrelated  with  each  other,  so  that  the 
total  contribution  to  the  far  field  is  simply  the  sum  of  the  individual 
volume  element  contributions. 

The  aero-acoustic  model  discussed  herein  is  directed  toward  prediction 
of  high  velocity  jet  noise  (1500  -  3000  feet  per  second),  for  arbitrary 
nozzle  shapes,  including  sound  pressure  level  spectra  at  any  observer  loca¬ 
tion. 


The  model  in  its  present  form  does  have  certain  limitations.  For 
multielement  configurations,  all  elements  must  be  parallel  and  non¬ 
impinging.  In  addition,  the  flow  field  calculation  assumes  constant  static 
pressure  mixing,  so  that  multichute  or  multitube  nozzles  with  significant  base 
pressure  variations  will  not  be  properly  simulated.  This  limitation  is  important 
for  closely-spaced  nozzle  elements.  Finally,  the  shock-cell  noise  portion  of  the 
prediction  is  likely  to  be  inadequate  for  multielement  and/or  multiflow  con¬ 
figurations. 


2.0  INTRODUCTION 


rv 


Many  jet  noise  suppressor  nozzles  have  been  designed  utilizing  intuitive 
notions  of  how  to  suppress  jet  noise,  and  have  demonstrated  substantial  noise 
reduction,  although  often  at  the  expense  of  considerable  thrust  loss, 
increased  engine  weight,  manufacturing  cost  and  complexity.  Seemingly  minor 
changes  in  suppressor  nozzle  design,  for  the  purpose  of  improving  thrust  per¬ 
formance,  often  result  in  substantial  loss  of  noise  suppression.  It  is  there¬ 
fore  highly  desirable  to  have  available  a  quantitative  prediction  technique 
for  estimating  the  aerodynamic  and  acoustic  characteristics  of  suppressor-type 
nozzle  configurations,  so  that  design  and  optimization  studies  can  be  made 
prior  to  construction  and  testing  in  order  to  minimize  the  time  and  cost  of 
development.  This  technique  should  ideally  be  sensitive  to  the  controllable 
design  variables,  and  should  be  based  on  a  minimum  of  empiricism.  Any 
empiricism  employed  should  be  tied  more  or  less  to  physical  (flow  and  acoustic 
propagation)  characteristics  rather  than  geometric  parameters,  i.e.,  it  should 
be  "universal"  in  a  normalized  (but  perhaps  restricted)  sense. 

Previous  work  on  modelling  of  jet  aero-acoustic  characteristics  has  been 
confined  to  simple  round  and  coannular  jets.  One  of  the  first  attempts  at 
developing  a  comprehensive  aerodynamic  acoustic  jet  model  was  published  by 
Lee,  Kendall,  et  al..  Reference  (1),  and  Grose  and  Kendall,  Reference  (2). 

This  approach  utilized  an  extension  of  Reichardt's  method  for  predicting  the 
jet  flow  properties  (this  method  is  adopted  herein  and  in  Reference  3)  for 
round  and  lobe-type  nozzles.  Acoustic  power  spectra  are  predicted  in 
Reference  (1)  and  (2),  based  on  a  "slice-of-jet"  model  wherein  the  power  per 
axial  slice  is  computed  and  related  to  a  certain  frequency  band  by  means  of 
empirically-derived  frequency  versus  axial  distance  relations.  Successful  pre¬ 
dictions  are  confined  to  low  Mach  numbers.  A  volume-element  summation  ("lump- 
of-jet")  approach  was  first  developed  by  Benzakein  et  al.,  Reference  (4),  for 
round  and  coannular  jets.  A  finite-difference  turbulent-kinetic-energy  model 
was  developed  therein  to  predict  the  jet  flow  field,  while  the  classical 
Lighthill  (5)  and  Ffowcs-Williams  (6)  formulations  were  employed,  with  suit¬ 
able  empirical  modifications,  to  predict  the  noise  from  each  volume  element. 
Extensions  of  the  method  of  Benzakein  et  al.  (4)  to  distinguish  between  self¬ 
noise  and  shear-noise,  proposed  by  Jones  (7)  and  Ribner  (8),  were  developed  by 
Knott  (9)  and  Moon  (10).  Recently,  Chen  (11)  has  applied  Kendall's  method 
(2)  to  predicting  power  spectra  of  coannular  jets. 

The  above  methods  either  ignore  the  sound/flow  interaction  effects  com¬ 
pletely,  or  recognize  only  source  convection  in  absence  of  a  shrouding  flow, 
which  has  been  shown  to  give  incorrect  simulation  (see  Mani,  Reference  12)  for 
all  but  the  lowest  jet  velocities,  especially  when  predicting  sound  pressure 
level  spectra  at  observer  angles  close  to  the  jet  axis.  The  aero-acoustic 
model  discussed  in  this  report  is  directed  toward  prediction  of  high  velocity 
jet  noise  (1500  -  3000  fps),  for  arbitrary  nozzle  shapes,  including  sound 
pressure  level  spectra  at  any  observer  location. 


Section  3  summarizes  the  mathematical  model  (which  is  presented  in  de¬ 
tail  in  Reference  3)  and  Section  4  describes  the  computer  program.  A  descrip¬ 
tion  of  inputs  to  the  computer  program  is  presented  in  Appendix  A  and  a  des¬ 
cription  of  the  outputs  is  presented  in  Appendix  B.  Appendices  C  and  D  con¬ 
tain  a  sample  output  and  a  computer  source  listing  for  the  program,  respec¬ 
tively. 


3.0  ANALYSIS 


3.1  OUTLINE  OF  METHOD 


Consider  a  jet  plume  exhausting  from  a  nozzle  of  arbitrary  shape,  as 
shown  in  Figure  1.  Utilizing  the  modified  Reichardt  theory  described  in 
Reference  (3),  the  mean  velocity,  temperature,  density  and  turbulent  shear 
stress  distributions  can  be  computed  throughout  the  jet  plume.  The  required 
inputs  are  nozzle  shape,  nozzle  exit  plane  total  pressures  and  temperatures, 
and  ambient  total  and  static  pressures  and  temperatures. 

The  jet  plume  is  subdivided  into  elemental  volumes  which  are  approxi¬ 
mately  the  size  of  a  typical  turbulent  correlation  volume  or  "eddy  size".  The 
modified  Reichardt  theory  can  provide  the  aerodynamic  properties  at  any  arbi¬ 
trary  point  in  the  plume  because  of  the  closed-form  solution  formulation  (the 
calculation  is  not  a  finite-difference  method  whose  grid  points  are  estab¬ 
lished/  dictated  by  the  accuracy/stability  requirements  of  the  numerical 
procedure).  The  flow  properties  are  therefore  computed  at  the  geometric 
centers  of  these  eddy  volume  elements. 

The  noise  generated  by  each  of  the  volume  elements  is  estimated  from  the 
classical  Lighthill  (5)  expressions  for  noise  produced  by  free  turbulence, 
assuming  that  the  turbulence  can  be  modeled  as  locally-isotropic ,  convecting 
quad rupole  sources  of  sound,  as  proposed  by  Ribner  (8).  The  turbulent  struc¬ 
ture  parameters  (intensity,  length-scale,  characteristic  frequency,  spectrum) 
required  for  computing  the  generated  noise  are  derived  from  the  calculated 
mean-flow  distributions  using  previously  established  empirical  similarity 
relations,  developed  by  Davies  et  al.  (13). 

For  each  volume  element,  the  effect  of  convection  and  fluid  shielding  on 
the  emitted  sound  of  that  element  is  computed.  The  flow  properties  in  the 
vicinity  of  the  element,  i.e.,  mean  velocity  profiles  and  temperature  profiles, 
determine  the  amount  of  flow  shrouding  or  fluid  shielding  seen  by  that  ele¬ 
ment.  From  the  generated  noise  spectrum  and  the  predicted  convection  and 
fluid  shielding,  the  net  emitted  noise  level,  at  each  observer  angle  and  1/3- 
octave  band  frequency  of  interest,  is  calculated.  The  contributions  from 
each  volume  element  in  the  jet  are  summed  on  a  mean-square  pressure  basis, 
assuming  that  individual  volume  elements  are  uncorrelated  with  each  other,  to 
provide  the  total  1/3-octave  spectrum  observed  in  the  far  field  sound  pres¬ 
sure. 


A  typical  example  of  how  a  jet  plume  is  subdivided  into  "eddy"  volume 
elements  is  illustrated  in  Figure  2.  The  subdivisions  are  very  small  close  to 
the  nozzle  exit  plane  where  the  turbulence  length-scales  are  small,  and  the 
volume  elements  are  made  progressively  larger  in  the  downstream  direction, 
simulating  the  increasing  length-scale  with  downstream  distance. 
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Figure  1.  Typical  Jet  Plume  Exhausting  from  a  Nozzle  of  Arbitrary 
(Noncircular)  Planform  Shape. 
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In  addition  to  the  calculation  of  turbulent  mixing  noise,  the  shock-cell 
broadband  noise  is  also  computed  using  a  modification  of  the  Harper-Bourne/ 
Fisher  method.  Reference  (14),  as  described  in  Reference  (3).  The  mixing 
noise  spectra  and  the  shock-cell  noise  spectra  are  summed  on  a  mean-square 
pressure  basis  to  yield  the  total  jet  noise  spectra,  at  each  far-field  obser¬ 
ver  angle. 

A  block  diagram  of  the  computation  sequence  for  the  jet  noise  prediction 
just  described  is  shown  in  Figure  3.  Note  also  that  atmospheric  attenuation 
corrections,  using  the  method  of  Reference  (15),  are  made  to  the  predicted 
spectra  to  account  for  air  attenuation  at  the  appropriate  far-field  distance. 
From  the  far-field  sound  pressure  spectra,  the  overall  sound  pressure  levels 
(OASPL),  perceived  noise  levels  (PNL)  and  sound  power  spectrum  (PWL)  are  also 
evaluated . 


3.2  FLOW  FIELD  PREDICTION 


The  jet  plume  flow  field  is  computed  using  the  extension  of  Reichardt’s 
method  developed  in  Reference  (3).  The  method  basically  consists  of  pre¬ 
dicting  the  diffusive  transport  of  momentum  flux  and  enthalpy  flux  from  a 
specified  exit  plane  distribution  to  various  axial  stations  along  the  plume. 

In  addition,  the  various  components  of  turbulent  shear  stress  are  also  calcu¬ 
lated,  being  related  to  directional  derivatives  of  the  axial  component  of 
momentum  flux.  From  these  distributions,  the  mean  axial  velocity,  density  and 
turbulence  intensity  distributions  are  estimated. 

A  typical  suppressor  nozzle  planform  shape  at  the  nozzle  exit  plane  is 
shown  in  Figure  1.  The  nozzle  contour  can  be  defined  by  coordinates  (vD,  a). 
The  jet  nozzle  exit  plane  conditions  are  denoted  by  subscript  "j",  and  ambient 
field  (external  flow)  conditions  by  subscript  "a".  The  flow  conditions  at  any 
flow  field  point  (x,  r,  j>)  are  computed  from  the  following  equations,  taken 
from  Reference  (3) : 


Momentum  Transport: 


2  2 

pir  -  p  u 

a  a 


=  ~  /( P  u2  -  p  U2) 

2  tt  J  j  j  a  a 


,.2.  I  .  -R2/b2  , 

-  P  il  )  I  1  -  e  o  m  |  d$ 
a  a 


Heat  Transport: 

pUH  = 


7“  /(pjUjV  t1  '  e'Rp/bh] 


di> 


(1) 


(2) 
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Figure  3.  Block  Diagram  of  Jet  Noise  Aeroacoustic 
Prediction  Model. 


Shear  Stresses: 


where 


denotes  a  contour  integral  around  the  nozzle  planform  boundary. 


In  the  above  equations,  U,  p,  and  H  are  the  mean  velocity,  density  and 
enthalpy,  respectively,  and  x  ,  xr  and  are  the  x,  r  and  <J>  -  components 
of  turbulent  shear  stress.  The  coordinates  R0  and  9  are  defined  by  the 
relations 


o 


(r 


2  r  v  [1 
o 


cos(tj>-a)  ] 


and  R  cos<J> 
o 


r-v  cos(<J>-a) 
o 


(6) 


The  turbulent  mixing  length  parameters  bm  =  Cmx  and  bh  =  Chx  are  determined 
by  the  momentum  spreading  rate  parameter  Cm  and  enthalpy  spreading  para¬ 
meter  Cft.  Empirical  expressions  for  these  have  been  developed  in  Reference 
(3),  and  are  as  follows: 


r  =  _ Q..  02 .5 _ 

m  (1  +  CjVR)  (1  +  C2Mj  ) 


1. 15C 

m 


(7) 
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where  VR  =  Ua/U.  and  Mj  -  Uj/Cj,  the  jet  exit  Mach  number.  The  constants 
Cx  and  C2  were  determined  frotn  calibrations  with  conical  nozzle  flow  field 
velocity  measurements,  and  values  of  Ci  =  0.25  and  C2  =  0.08  were  found  to 
give  the  best  agreement  with  experiments. 


The  velocity  and  density  are  determined  from  the  distributions  of 
(pU^)  and  (pUH)  from  the  following  expressions: 


U 


where 


and 


and 


pUH 
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H 
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+  ♦ 

1  2 

=  CpT  +  j  pVZ 

Hp  +  1  (pU 
(pu2)  .  T 

U  2  •  T 


CpTa(pU2) 


-  C  T 
P  a 
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(9) 
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(11) 


In  the  above  expressions,  Cp  and  Rg  are  the  specific  heat  at  constant 
pressure  and  universal  gas  constant,  respectively  (Rg  =  yC  /(y-1)),  y  is 
the  ratio  of  specific  heats,  and  p  is  the  jet  static  pressure,  assumed  to 
be  equal  to  the  ambient  static  pressure.  The  local  temperature  is  T,  and 
Ta  is  the  ambient  static  value. 

The  axial  turbulence  intensity  u'  is  computed  from  the  shear  stress 
components  utilizing  the  following  expression: 

(u’)2  =  ^/(tr/p)2  +  (t^/p)2+  et  (10  tx/p)2  (12) 

where  6t  and  8t  are  empirically  -  determined  constants. 

Modifications  to  the  above  computation  procedure  are  developed  in 
Reference  (3)  for  nozzles  with  an  axisymmetric  centerbody  or  plug.  This 
modification  consists  of  a  coordinate  transformation  of  the  variables  (Rq, 
4>)  as  follows: 
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2  2 

R  =  (Ar-Av)  +  2ArAv  [  1  -  co  s  (  d>-a  )  ] 
o 

R  cosi  =  Ar  -  Av„  cos(4>-a) 

o  _ 0  , - 

where  Ar  =  ^/r2  -  r^  ( x )'  and  Av  ~  yj (o) 

and  r{,(x)  is  the  centerbody/plug  contour  coordinate  specification.  For 
rj,(x)  =  0,  Equation  (13)  reduce  to  Equation  (6,). 

Equations  (1)  through  (13)  completely  define  the  flow  field  calculation 
procedure.  As  discussed  in  the  previous  section,  the  flow  properties  are 
evaluated  at  all  field  points  corresponding  to  eddy  volume  element  centers 
(x,r  ,<t>) 


3.3  TURBULENT  MIXING  NOISE  PREDICTION 


Expressions  are  developed  in  Reference  (3),  for  the  far-field  noise  of 
converted  quadrupoles  imbedded  in  a  parallel  shear  flow,  utilizing  high- 
frequency  asymptotic  solutions  to  Lilley’s  equation.  These  expressions, 
for  a  source  of  unit  volume  and  unit  strength,  when  applied  to  a  collection 
of  sources  distributed  throughout  a  parallel  shear  flow  model  of  the  jet 
plume,  yield  the  following  equation  for  the  far-field  noise  spectrum: 

p2  (R.e.n)  =  /  a(a  +  4a  +  2a  +  2a  )dy  (1^ 

v  *«»“.'  j  xx  xy  yy  yz 

T 

where  /  (  )dy  implies  integration  over  the  entire  jet  plume.  The  factors 
axx,  axy»  ayy  and  ayz  are  the  directivity  factors  for  each  of  the  contri¬ 
buting  quadrupole  types  contained  in  each  turbulent  eddy  volume  element. 

The  factor  a  is  given  by 


a  -  (— )2(^)2  (1-H  COSfl )  ”hl-M  co 

16*ZR2Ca4  V  po>  Xc0  ’  0  c 


s  0 ) 


(15) 


where  1(D)  is  given  by 

1(D)  =  p2  (u')2(Dt0)^  exp  £-■ j |  (ftt0)2J 


(16) 


and  t0  is  the  two-point  velocity  correlation  characteristic  time  delay. 
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Subscript  "o"  refers  to  the  volume-element  or  eddy-center  location  conditions. 
The  parameters  M0  and  Mc  are  defined  as 


M  =U  /C  M  =U  /C 

o  o  a  c  c  a 


where  Uc  is  the  effective  convection  velocity  of  the  eddy. 

The  directivity  factors  aXx,  aXy,  3yy  and  a^z  have  different  forms, 
depending  on  the  location  of  the  source  and  the  velocity  and  temperature 
profiles  in  the  vicinity  of  the  source.  As  discussed  in  Reference  (3),  these 
factors  depend  explicitly  upon  a  shielding  function  g2 ,  which  has  the  form 


2,  N  (1-M  cos0)2(Ca/C) 2  -  cos^d 

8  (r)"  (1-M  cos8)2 

c 


(17) 


where  C=C  (r)  and  M=M(r)/Ca.  Given  a  velocity  U(r)  and  temperature  profiles 
C(r)  ~  /T(r) ,  the  shielding  function  profile  can  be  computed  per  Equation 
(17).  Defending  on  the  observer  angle  6  and  the  profile  shapes,  the  pro¬ 
file  of  g  (r)  may  have  both  positive  and  negative  "zones",  or  may  be  posi¬ 
tive  for  all  values  of  r.  If  a  negative  region  exists,  then  the  possibility 
of  fluid  shielding  of  the  sound  scarce  exists,  depending  on  the  location  of 
the  source  relative  to  the  negative  or  shielding  zone.  The  location  r=rD 
where  gZ (r)  crosses  zero  is  termed  a  turning  point;  practically  speaking, 
more  than  one  turning  point  can  occur,  although  more  than  two  turning 
points  is  very  rare. 

A  maximum  of  two  turning  points  is  considered  herein.  This  allows  six 
possible  situations  regarding  the  source  location  relative  to  a  shielding 
zone.  These  are  illustrated  qualitatively  in  Figure  4.  This  figure  shows 
g2(r)  plotted  versus  radial  distance  in  the  jet  plume  for  six  cases.  In 
case  (a),  g^(r)  is  positive  for  all  values  of  r,  so  no  shielding  of  the 
source  occurs  no  matter  where  it  is  located  radially  in  the  jet.  In  case 
(b),  the  source  is  located  outboard  of  the  region  where  g^(r)  is  negative, 
so  still  no  shielding  occurs.  In  case  (c) ,  however,  the  source  is  located 
inboard  of  the  turning  point,  rG  <  r0-^,  and  lies  inside  the  shielding  zone. 
The  source  will  therefore  be  shielded  to  some  extent. 

In  cases  (d,e,f)  shown  in  Figure  4,  there  are  two  turning  points  roy 
and  ra2  >  rai.  In  case  (d),  the  source  lies  outside  of  both  turning  points 
and  therefore  sees  no  shielding.  In  case  (e) ,  the  source  lies  between  the 
two  turning  points  and  thus  Is  shielded  by  the  fluid  layer  between  r=ro  and 
r=ro2. 
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No  Turning  Points 


One  Turning  Point 


One  Turning  Point 


Dimensionless  Radial  Coordinate  r/a 


Figure  4.  Possible  Solution  Types  for  a  Maximum  of  Two  Turning 
Points.  Shaded  Areas  Denote  Shielding  of  Source. 


The  acoustic  radiation  solutions  for  all  of  the  six  cases  illustrated  in 
Figure  4  are  worked  out  in  detail  in  References  (3),  and  yield  different  forms 
for  the  directivity  factors  axx,  &xy  >  ayy  and  ay2  for  each  case.  These 
directivity  factors  have  the  following  forms: 
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2  2 
g  cos  6 
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2  2 

where  gQ  is  the  value  of  g  (r)  at  the  source  radius  r=rQ.  The  shielding 
coefficients  Bxx»  &xy>  8yy  and  Byz  depend  upon  the  case  encountered  in 
Figure  4.  If  the  parameters  Bqi»  $02  and  812  are  defined  as 
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exp 

-2k 
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|g2(r) 

1/2  dr 
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where  k=P./Ca,  the  shielding  coefficients  are  then  determined  from  the 
following  table: 
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Table  1.  Shielding  Coefficients 


Case 


^xx 

1 

1 

601 

1 

B02 


e 


12 


6xy 

1 

1 

0 

1 

0 

e12 


6 


yy 

l 

l 

o 

l 

0 

312 


6 


yz 

1 

1 

0 

1 

0 

312 


Note  that  a  value  of  Bxx,  etc.,  of  unity  indicates  no  fluid  shielding,  as  in 
cases  (a),  (b)  and  (d) .  When  the  source  is  imbedded  within  the  shielding 
zone,  as  in  cases  (c)  and  (e),  only  the  x-x  quadrupole  contributes,  for 
reasons  explained  in  Reference  (3). 


3. A  SHOCK-CELL  NOISE  PREDICTION 


As  discussed  in  Section  3.1,  the  shock  cell  broadband  noise  pre¬ 
diction  is  based  on  a  modification  of  the  theory  of  Harper-Bourne  and 
Fisher  (14) .  A  thorough  discussion  of  this  theory  and  its  application  to 
noncircular  nozzles  can  be  found  in  Reference  (3). 

Although  the  analysis  of  Reference  (3)  demonstrated  that  noncircular 
nozzle  shock-cell  noise  exhibits  the  same  scaling  of  noise  level  with 
operating  conditions  as  that  for  a  conical  nozzle,  the  influence  of  nozzle 
shape  on  noise  level  and  spectrum  shape  was  not  quantified  to  the  extent 
that  a  verified  prediction  method  could  be  established.  A  method  for 
predicting  the  shock  cell  structures  (number  of  cells,  spacing,  etc.)  is 
still  required  before  a  general  shock  noise  prediction  procedure  can  be 
established.  An  interim  shock  cell  noise  prediction  method  was  therefore 
adopted  for  incorporation  into  the  unified  aeroacoustic  prediction  model, 
with  the  expectation  of  replacement  by  a  more  general  method  at  some  future 
dace . 


The  interim  method  is  essentially  that  of  Deneuville  (16),  with  some 
modifications  to  simplify  the  calculation  and  incorporate  some  of  the  ideas 
developed  in  Reference  (3).  The  method  consists  of  modeling  the  shock  cell 
noise  component  of  the  spectrum  by  two  straight  lines,  as  illustrated  in 
Figure  5.  The  primary  variables  required  are  peak  sound  pressure  level 
SPLp  and  the  frequency  fp  at  which  this  occurs. 
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SPL,  dB 


Given  a  nozzle  operating  at  pressure  ratio  PR,  jet  velocity  Vj, 
having  total  flow  area  A  j  ,  equivalent  diameter  Dg  =  /AAj / n ,  and  hydraulic 
diameter  =  4Aj/Pw  (where  Pw  is  the  nozzle  boundary  wetted  perimeter), 
the  parameters  SPLp  and  fp  can  be  computed  as  a  function  of  observer  angle 
0  ^ .  The  computation  proceeds  as  follows: 


(1)  compute  shock  strength  parameter  8  from 


8 

-Ah2?  -*] 


(20) 

(21) 


(2) 


compute  average  shock-cell  spacing  L 


avp 


from 


L  =  1.1  6  D 

avg  eq 


(22) 


(3)  Compute  peak  noise  frequency  from 
U 

fp  =  7 —  a  +  Mc  cos  e.) 

‘"avg 

where  U  =  0.7  U.  and  M  =  U  /C  . 

c  j  c  c  a 


(23) 


(4)  Compute  peak  SPL  from 

SPLp  (0i)  =  152.6  +  40  log1Q(0)  +  10  log1Q  (Aj/R2) 

+  10  log10  (Dh/Deq)-40  loglO(l-Macos0i)  (24) 


+  10  log1()  (N/8) 

where  Ma  =  Ua/Ca,  the  flight  Mach  number,  and  N  is  the  number  of  shock 
cells  (usually  N  =  8  is  assumed) .  The  spectrum  shape  is  then  calculated 
from  the  assumed  two  straight-line  model  shown  in  Figure  5.  The  equations 
are  as  follows: 


SPL(f,ei)  =  SPLp(9i)  -  10  log10(f/fp),  f  >  fp 
SPL(f,91)  =  SPLp(01)  +  70  log10(f/fp),  f  <  fp 


(25) 
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Equations  (20  through  25)  completely  describe  the  shock  cell  noise  predic¬ 
tion  method  for  obtaining  1/3-octave  spectra  at  any  observer  angle  0^. 


3.5  AEROACOUSTIC  MODEL  INTEGRATION 


An  outline  of  the  basic  prediction  method  and  a  description  of  the 
component  building  blocks  have  been  presented  in  the  preceding  sections. 
This  section  describes  how  these  building  blocks  are  tied  together.  Addi¬ 
tionally,  some  practical  guidelines  are  given  which  were  found  from  exper¬ 
ience  to  be  helpful  in  producing  reasonably  accurate  predictions  while 
maximizing  computational  efficiency. 

The  mixing  noise  spectrum  is  computed  from  equation  (14),  with  the 
integration  over  the  jet  plume  }  (  )  dy  replaced  (or  approximated)  by  a 
summation  over  all  eddy-volume  elements.  Since  equation  (14)  represents 
the  narrowband  spectrum  in  terms  of  source  frequency  SI  (emitted  frequency 
in  a  reference  frame  moving  with  the  eddy) ,  a  conversion  is  made  to  1/3- 
octave  spectrum  based  on  observed  center  frequency  f.  It  is  assumed  that 
the  1/3-octave  band  level  can  be  approximated  by  the  narrowband  level 
evaluated  at  the  center  frequency,  multiplied  by  the  bandwidth,  rather  than 
integrating  the  narrowband  level  distribution  over  the  bandwidth.  The 
standard  1/3-octave  center  frequencies  are  used.  For  each  eddy-volume 
element,  at  each  observer  angle  0^  =  180°  -  0,  the  source  frequency  is 
calculated  for  each  1/3-octave  center  frequency  from  the  relation 

fi  =  2rf(l-Mc  cos  0)  (26) 


The  spectrum  shape  of  a  given  eddy,  equation  (16) ,  determines  the  amount 
contributed  by  a  given  eddy  at  that  1/3-octave  frequency.  Theoretically, 
all  eddies  contribute  some  amount  at  all  1/3-octave  frequencies,  but  the 
rather  "peaky"  nature  of  I  (SI)  given  by  equation  (16)  dictates  that  the 
major  contribution  of  an  eddy  will  be  in  the  vicinity  of  fl  ~  4/tq. 

Computationally,  the  integrand  of  equation  (14)  can  be  expressed 
explicitly  in  terms  of  observed  frequency,  since  the  source  frequency 
always  occurs  in  the  combination  SI/ (1-Mccos0) .  The  one  exception  is  in  the 
exponent  of  the  spectrum  function  I(S2),  equation  (16).  For  example,  the 
product  aaXx  can  be  rewritten  as 

_  Pa2(u')7(2irfTO)^6xxcos40  (Ca/C0)2  r  ,  2 

aaxx  - j— - - - - - -  exp  -  qr  (S2t0) 

16it  R 4  C  (1-M  cos0)  (1-M  cos0)  L 
a  o  c 

Similar  expressions  can  be  derived  for  otaxy,  etc. 
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The  eddy  convection  factor  (1-Mccos0)  has  a  singularity  at  Mccos0  =*  1. 

To  circumvent  this  computational  difficulty,  it  is  replaced  by  a  modified  con¬ 
vection  factor  as  suggested  by  Ffowcs-Williams  (6)  and  Ribner  (8),  as  follows: 


C  =  1-M  cos0  -*/( 1-M  cos0)2  +  (a  u'/C  )2  (27) 

c  \  c  t  a 

where,  as  Ffowcs-Williams  and  Ribner  have  shown,  the  term  (atu’/ca)  accounts 
for  the  finite  life-time  of  the  eddy  as  it  is  convected  downstream.  The 
constant  at  was  determined  (from  comparison  of  prediction  with  experiments) 
to  be  approximately  0.5,  independent  of  source  location,  jet  operating  con¬ 
ditions  and  nozzle  geometry.  An  additional  assumption  was  made  that  the  flow 
convection  factor  (1-Mocos0)  can  be  replaced  by  the  modified  eddy  convection 
factor  C  given  by  equation  (27). 

The  eddy  convection  Mach  number  must  be  a  function  of  the  local  flow 
Mach  number  of  the  eddy-volume  being  considered.  Several  expressions  for  Mc 
were  tried;  the  best  results  (over  a  wide  range  of  nozzle  operating  condi¬ 
tions  and  geometries)  were  obtained  with  the  following: 

M  =  (M  +  0.65  M.)  (28) 

c  2  o  j 

where  Mj  =  Uj/Ca,  the  nozzle  exit  acoustic  Mach  number.  Equation  (28)  repre¬ 
sents  a  simple  average  of  the  classical  assumption  Mc  =  0.65  Mj  and  the  local 
Mach  number  M0  =  M(r0).  For  suppressor  nozzles,  equation  (28)  works  best  if 
Mj  is  replaced  by  the  postmerged  potential  core  Mach  number  Mm  =  Um/Ca. 

This  can  be  evaluated  from  the  results  of  the  flow  field  calculation  des¬ 
cribed  in  Section  3.2. 

The  acoustic  theory  summarized  in  Section  3.3  applies  only  for  axi- 
symmetric  jets.  For  nonaxisymmetric  nozzles  (multilobe,  chute,  tube,  etc.), 
an  assumption  is  therefore  made  that  a  representative  average  radial  profile 
at  each  axial  station  can  be  derived  which,  when  inserted  in  the  acoustic 
calculation,  will  adequately  model  the  acoustic  characteristics  of  the  asym¬ 
metric  jet.  The  mass-averaged  values  of  U  and  p  are  calculated  from  the 
azimuthal  average  of  the  quantities  pU^  and  pUH.  The  resulting  distributions 
of  U(x,r)  and  p(x,r)  are  then  employed  in  the  acoustic  calculation  described 
in  Section  3.3. 

For  axisymmetric  jets,  the  empirical  similarity  relation  [tg  -  (3U/3r) 
derived  by  Davies  et  al.(13)  can  be  used  to  calculate  the  local  value  of  td 
from  the  radial  velocity  profiles  U(r).  For  a  multichute  or  lobe  nozzle, 
however,  there  are  volume  elements  close  to  the  nozzle  exit  plane  which 
have  negligible  radial  gradients  9U/9r  but  large  azimuthal  gradients  3U/3ij>. 

It  is  therefore  assumed  that  t0  is  a  function  of  both  3U/9r  and  3U/ 3<J> .  From 
Reichardt's  hypothesis  (see  Reference  3),  the  transverse  shear  stresses  are 
related  to  9U/3r  and  3U/ 3 <p  by  the  approximate  formulae 
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Il 

p 


and 

P 


\  9U2 

r  3  4> 


2 

where  X  =  Cm  x/2.  Thus  the  transverse  derivatives  of  U  can  be  expressed  in 
terms  of  Tr  and  as  follows: 


3U  ,  (tr/pi 

3r  U  C?  x 
m 


and 


1  3U 
r  3<J> 


IWp 2 

U  C2  x 


m 


A  new  transverse  derivative  3U/3n  is  therefore  defined  in  terms  of  the  above 
r  and  derivatives  as  follows. 


3U 

3n 


(29) 


such  that 


(30) 


The  parameter  yt  is  an  empirical  constant  which  must  be  evaluated  from  com¬ 
parison  with  experiments. 

For  prediction  of  jet  noise  in  flight,  the  mean  square  sound  pressure 
level  should  be  multiplied  by  the  dynamic  amplification  factor  (1  +  Macos0)  , 
where  Ma  =  Ua/Ca.  In  addition,  the  convection  and  flow  Mach  numbers  are 
replaced  by  (Mc-Ma)  and  (MQ-Ma) ,  respectively,  where  Mc  and  M0  are  evaluated 
in  a  reference  frame  fixed  to  the  nozzle. 

Finally,  in  all  predictions  of  1/3-octave  spectra,  the  atmospheric 
attenuation  corrections  given  by  Reference  (15)  are  applied,  using  standard- 
day  corrections  (70%  relative  humidity  and  59°  F  dry-bulb  temperature) 
evaluated  at  the  center  frequency. 
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4.0  COMPUTER  PROCRAM  DESCRIPTION 


4.1  INTRODUCTION 


The  central  theme  of  this  report  has  been  the  development  of  a  general 
method  for  predicting  the  aeroacoustic  characteristics  of  turbulent  jets  for 
a  wide  range  of  nozzle  flow  conditions  and  nozzle  geometries.  The  preceding 
section  describes  the  analytic,  algebraic  relationships  which  have  been 
developed  to  represent  the  various  physical  aspects  of  the  mathematical 
model.  The  purpose  of  this  section  is  to  describe  the  computational  algo¬ 
rithms  and  associated  computer  program  which  provides  the  necessary  link 
between  the  symbolic  representation  of  the  model  and  the  actual  numerical 
results  of  the  prediction  method. 

The  computer  program  is  written  in  FORTRAN  IV  language.  It  has  been  run 
on  both  the  GE/Honeywell  6080  and  GDC  7600  computers,  and  can  easily  be  modi¬ 
fied  for  running  on  other  systems.  The  program  subdivides  the  jet  plume 
utilizing  a  built-in  grid  system  which  requires  minimal  input  for  specifi¬ 
cation.  This  grid  system  can  be  superseded  by  the  user  through  more  complex 
input  if  desired.  The  nozzle  geometry  is  input  through  discrete  point  coor¬ 
dinates  for  each  nozzle  element  boundary,  and  up  to  109  elements  can  be  input 
for  a  given  case.  A  maximum  of  24  axial  stations  along  the  jet  plume  is 
permitted,  and  up  to  200  radial  points  per  axial  station  can  be  accommodated. 

These  limits  can  be  changed  if  so  desired  by  modifying  the  appropriate 
DIMENSION  and  COMMON  statements  in  the  program  logic. 

The  limiting  assumptions  made  in  developing  the  method  have  been  dis¬ 
cussed  in  Section  3.0,  but  it  is  appropriate  to  summarize  them  here  to  warn 
against  indiscriminate  violation  of  these  limitations.  They  are  as  follows: 

1.  The  exhaust  nozzle  elements  should  be  co-planar,  i.e.,  each  tube  or 
chute  of  a  multielement  configuration  should  exhaust  at  the  same  axial 
plane.  However,  nozzle  element  exit  planes  can  be  staggered,  provided 
that  the  mixing  layer  of  a  given  element  does  not  impinge  on  the  wall 
of  another  element. 

2.  The  jet  exhaust  gases  must  all  be  of  the  same  constituent,  because 
the  calculation  cannot  accommodate  gas  mixtures  or  species  concen¬ 
trations. 

3.  Within  any  nozzle  element,  the  flow  is  assumed  to  be  uniform  at  the 
exit  plane. 

4.  The  time- averaged  static  pressure  is  assumed  to  be  constant  and 
uniform  throughout  the  jet  flow  field  and  surrounding  ambient 
field. 

5.  The  exhaust  nozzle  elements  must  discharge  axially,  radial  mean 
flow  and  swirl  are  neglected  in  the  model. 

i 

1 
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6. 


The  effects  of  shock  formations  on  mixing  and  turbulence  levels  are 
neglected. 


These  assumptions  and  limitations  are  those  which  pertain  to  the  types 
of  problems  which  can  be  analyzed.  There  are,  of  course,  additional  assump¬ 
tions  which  went  into  the  formulation  of  the  model  itself  which  may  restrict 
the  accuracy  of  the  model,  but  not  the  type  of  problem  which  can  be  analyzed 
The  user  is  advised  to  consult  Reference  (3)  for  the  details  of  the  model 
formulation. 


4.2  PROGRAM  NOMENCLATURE  AND  SYMBOL  CONVENTION 


The  jet  plume  and  nozzle  geometry  coordinates  are  computed  in  the  MAIN 
routine.  The  jet  plume  is  divided  into  KX  axial  slices,  specified  by  KA 
(1  £  KA  £  KX) .  The  FORTRAN  symbol  variables  for  the  various  coordinate 
parameters  and  indices  are  shown  in  Figure  6.  Note  that  the  radial  sub¬ 
division,  specified  by  index  M  (1  £  M  £  NR),  proceeds  in  increments  DSIG(KA), 
from  SIG  =  RMIN(KA)  to  the  maximum  value  set  by  NR.  The  value  of  NR  is 
determined  during  the  calculation  from  the  location  where  the  axial  momentum 
flux  is  within  a  certain  tolerance  of  being  equal  to  the  ambient  level,  i.e., 

| RU2  -  RU2E (1) |  £  RU2M 

where  RU2M  is  a  specified  input  tolerance.  The  maximum  allowable  value  of  NR 
can  be  specified  by  the  input  variable  IQUIT.  The  program  dimension  sizes 
limit  KX  and  IQUIT  to  the  following  maximum  values: 

KX  £  24  IQUIT  £  200 

The  nozzle  geometry  itself  is  input  as  a  number  (NEST)  of  boundary  ele¬ 
ments.  Each  element  is  specified  by  coordinate  pairs  RA(I,J)  and  DALP(I,J), 
where  RA(I,J)  denotes  the  radius  and  DALP(I,J)  denotes  the  angular  increment, 
as  shown  in  Figure  6.  The  index  I  denotes  the  boundary  contour  point  number, 
and  the  index  J  denotes  the  boundary  number.  The  reference  angular  location 
for  each  boundary  is  given  by  ALPO(J).  For  each  boundary,  the  exit-plane 
values  of  total  pressure  PT(J)  and  total  temperature  TT(J)  are  also 
specified.  Boundary  number  one  (J=l)  is  always  considered  to  be  the  ambient 
field. 

The  far  field  acoustic  calculations  are  performed  on  either  a  constant- 
radius  arc  or  a  sideline  parallel  to  the  jet  axis,  according  to  whether  the 
input  variable  NUMANG  is  set  equal  to  1  or  2,  respectively.  For  NUMANG  =  1, 
the  input  DIST  is  the  arc  radius;  for  NUMANG  =  2,  DIST  is  the  sideline  dis¬ 
tance.  The  acoustic  arena  geometry  specification  in  terms  of  FORTRAN  vari¬ 
ables  is  shown  in  Figure  7.  Note  that  a  distinction  is  made  between  the 
source-to-observer  distance  RSTAR  and  the  nozzle-to-observer  distance  RADIUS. 
The  observer  angle  relative  to  the  jet  axis  THETA  is  always  in  units  of 
radians,  while  the  observer  angle  relative  to  the  inlet  axis  THETD  is  in 
units  of  degrees.  The  far-field  sound  pressure  level  SPL(I,J)  is  computed  at 
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(a)  Radial/Axial  Plane 


Figure  6.  FORTRAN  Symbol  Convention  for  Coordinates  and 
Geometric  Variables. 
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Observer 


Figure  7.  FORTRAN  Symbol  Convention  for  Acoustic  Arena  Variables. 
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every  1/3-octave  frequency  from  FMIN  to  FMAX,  at  ten-degree  increments  from 
THEt'd  =  20°  to  160°. 


A  list  of  the  important  FORTRAN  symbols  used  in  the  computer  program  is 
given  in  Table  2,  along  with  their  algebraic  equivalents  where  possible.  A 
complete  description  of  all  of  the  input  variables  and  examples  of  input 
preparation  are  given  in  Appendix  A. 


4.3  DESCRIPTION  OF  PROGRAM  AND  SUBROUTINES 


A  flow  chart  of  the  computer  program  logic  is  shown  in  Figure  8.  It 
indicates  the  sequence  of  operations,  the  interconnections  of  various  portions 
of  the  program,  and  their  functions.  A  description  of  the  main  program  and 
each  of  the  subroutines  is  given  in  the  following  paragraphs. 


4.3.1  MAIN 

The  main  program  initiates  the  computation  and  controls  the  sequence  of 
operations.  It  reads  the  input  data,  computes  the  grid  system  for  the  aero¬ 
dynamic  flow  field,  and  computes  the  various  required  nozzle  exit  plane  flow 
parameters  such  as  velocities,  Mach  numbers,  momentum  and  enthalpy  fluxes, 
etc.  The  main  program  prints  out  all  input  data,  nozzle  exit  conditions, 
nozzle  geometry,  and  coordinate  system  parameters. 

The  main  program  controls  and  executes  the  jet  plume  flow  field  compu¬ 
tation.  After  each  axial  slice  has  been  evaluated,  the  MAIN  program  calls 
subroutine  SLICE  to  perform  the  requested  acoustic  calculations.  Upon  com¬ 
pleting  the  calculations  at  all  axial  slices,  MAIN  then  calls  OUTPUT  to  per¬ 
form  some  final  calculations  and  print  out  the  far-field  noise  levels.  If 
additional  cases  are  requested,  the  entire  procedure  is  repeated,  beginning 
with  reading  of  input  data;  otherwise  the  execution  is  halted. 


4.3.2  ARCCOS(X) 

This  is  a  function  subroutine  which  computes  the  principal  value  of  the 
arc  cosine  of  the  variable  X,  It  is  used  in  MAIN  in  evaluating  certain 
angles  relating  boundary  coordinate  points  and  flow  field  location  points. 


4.3.3  ERF (X) 

This  function  subroutine  evaluates  the  error  function  of  argument  X 
using  polynomial  approximations  as  given  in  Reference  (19).  It  is  used  in 
MAIN  for  evaluating  flow  field  integrands. 


Table  2.  List  of  FORTRAN  Symbols. 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

AA 

Air  attenuation  factor 

ATMOS 

AAA 

Intermediate  variable 

LSPFIT ,  MAIN 

ABDTH 

|  A4>| 

MAIN 

ABLE 

Intermediate  variable 

MAIN 

ABPA 

I  <}>-a| 

MAIN 

ACH 

Mach  number  M 

MAIN 

ACHM 

Average  mach  number 

MAIN 

ACH2 

M2 

MAIN 

AK 

Sound  level  constant 

MAIN,  OUTPUT 

AL 

Lighthill  parameter 

MAIN 

ALFA 

Frequency  constant 

MA'N 

ALP 

Angle 

MAIN 

ALPHT 

Convection  constant  at 

SLICE 

ALP0 

Reference  boundary  angle 

MAIN 

AMUIN 

Input  turbulence  constant  pt 

MAIN 

AMULT 

Intermediate  value  for  Pt 

MAIN 

AO 

Speed  of  sound  Ca 

MAIN 

ATOTAL 

Total  flow  area 

MAIN 

B 

Intermediate  variable 

LSPFIT 

BETA 

Shock  strength  parameter  B 

SHOCK 

BETAIN 

Input  turbulence  constant  6t 

MAIN 

BETAMC 

Convection  constant  Bmc 

MAIN,  SLICE 

BK 

Intermediate  variable 

SLICE 

BKR 

Intermediate  variable 

MAIN 

BOT 

Intermediate  variable 

LSPFIT 

BUG 

Intermediate  variable 

MAIN 

C 

Constant 

LSPFIT 

CH 

Spreading  parameter  Ch/Cm 

MAIN 

CHX 

Spreading  parameter  Chx 

MAIN 

CJOCO 

Ratio  of  Cj/Ca 

SLICE 

CM 

Spreading  parameter  Cm 

MAIN 

CMAX 

Intermediate  variable 

TPNLC 

CMC 

Intermediate  variable 

MAIN 

CMMC 

Spreading  constant  C^ 

MAIN 

CMVR 

Spreading  constant  C2 

MAIN 

CNST 

Constant 

SLICE 

CO 

Ambient  speed  of  sound  Ca 

MAIN,  SLICE,  SHOCK 

COEF 

Conversion  factor 

OUTPUT 

CONV 

Convection  factor 

SHOCK 

CONVF 

Flight  dynamic  factor 

SLICE 

CON  VO 

Convection  factor 

SLICE 

CONV2 

Modified  convection  factor  C 

SLICE 

CONI 

Constant 

SLICE 
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Table  2.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 


Symbol 

Meaning 

Related  Subroutines 

C0N2 

Constant 

SLICE 

COST 

Cos  $ 

MAIN 

COSTO 

Cos  $ 

MAIN 

CP 

Specific  heat  Cp 

MAIN 

CT 

Cos  9 

SLICE,  CRD 

CTSQ 

Cos2  0 

SLICE 

CTH 

Cos  0 

SHOCK 

CVR 

Intermediate  variable 

MAIN 

DALP 

Boundary  coordinate  Aot 

MAIN 

DDTHE 

Tolerance  on  A6,  radians 

SLICE 

DDTHED 

Tolerance  on  A0,  degrees 

SLICE 

DELRA 

Transformed  boundary  radius  Av 

MAIN 

DELSIG 

Transformed  radius  Ar 

MAIN 

DELTA 

Turbulence  constant  <5t 

MAIN 

DELTIN 

Input  array  of  6t 

MAIN 

DEQ 

Equivalent  diameter  Deq 

MAIN,  SLICE,  SHOCK 

DIA 

Reference  Deq 

MAIN 

DIRECT 

Directivity  factor 

SLICE 

DIST 

Sideline  or  arc  distance 

MAIN,  SLICE 

DJET 

Reference  diameter 

MAIN 

DPHI 

A<(> 

MAIN 

DRMIN 

Ar  -  minimum  value 

SLICE 

DS 

Source  strength  amplitude 

SLICE 

DSIG 

Ar 

MAIN,  SLICE 

DSPL 

Mixing  noise  pressure 

SLICE,  OUTPUT 

DSPLl 

Intermediate  variable 

SHOCK 

DSPL2 

Intermediate  variable 

SHOCK 

DTHED 

A0,  degrees 

SLICE 

DTHM 

Maximum  increment  of  <(> 

MAIN 

DU 

Intermediate  variable 

MAIN 

DUDR 

3U/3r 

MAIN,  SLICE 

DV 

Eddy  volume  dV 

SLICE 

DX 

Axial  step  size  Ax 

MAIN,  SLICE 

EF 

Enthalpy  flux 

MAIN 

EFE 

Enthalpy  flux 

MAIN 

EM 

Mach  number 

SLICE 

EMACH 

Exit  Mach  number 

MAIN,  SLICE,  OUTPUT 

F 

Intermediate  variable 

LSPFIT 

FAC 

Intermediate  variable 

PNLC 

FC 

Center  frequency 

SLICE 

FIRSTU 

Flight  velocity  Ua 

MAIN,  SLICE 

FIS 

Intermediate  variable 

MAIN 

FM 

Mass  flow 

MAIN 
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Table  2.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

FMAX 

Maximum  observed  frequency 

MAIN,  OUTPUT 

FMIN 

Minimum  observed  frequency 

MAIN,  OUTPUT 

FO 

Observed  frequency 

SLICE,  SHOCK,  OUTPUT 

FP 

Peak  frequency 

SHOCK 

FR 

Frequency  ratio 

SLICE 

FRSQ 

Intermediate  variable 

SLICE 

FS 

Source  frequency 

SLICE 

GAM 

Specific  heat  ratio  y 

MAIN,  SHOCK 

GAMA 

Gas  constant  parameter 

MAIN 

GEXP 

Gas  constant  parameter 

SHOCK 

GM 

Shielding  function 

CRD 

GOSQ 

Shielding  function 

CRD 

G2 

Shielding  function 

SLICE,  CRD 

HF 

Spectrum  function 

SLICE 

HPSI 

Intermediate  variable 

MAIN 

HTR 

Stagnation  enthalpy 

MAIN 

I 

Index 

ALL 

IC 

Index 

LSPFIT 

I  DENT 

Title  (80-characters) 

MAIN 

II 

Index 

TPNLC 

IMH 

Index 

MAIN 

IQUIT 

Maximum  number  of  points 

MAIN 

IS 

Index 

MAIN 

ISSY 

Index 

MAIN 

I  SAVE 

Index 

LSPFIT 

ISYM 

Symmetry  indicator 

MAIN 

IT 

Symmetry  indicator 

MAIN 

J 

Index 

ALL 

TMAX 

Maximum  band  number 

OUTPUT,  SHOCK,  SLICE 

JMIN 

Minimum  band  number 

OUTPUT,  SHOCK,  SLICE 

J1 

Index 

CRD 

Jll 

Index 

CRD 

J2 

Index 

CRD 

J21 

Index 

CRD 

J211 

Index 

CRD 

K 

Index,  also  wave  number 

MAIN,  SLICE,  PNLC 

KN 

Surrounding  boundary  index 

MAIN 

KNCAS 

Case  counter 

MAIN 

KNK 

Surrounding  boundary  index 

MAIN 

KX 

Number  of  axial  slices 

MAIN 

L 

Index 

MAIN 

LAVG 

Shock  spacing 

SHOCK 

LEAF 

Number  of  boundary  leaves 

MAIN 
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Table  2.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 


Symbol 

Meaning 

Related  Subroutines 

LEAV 

Number  of  boundary  leaves 

MAIN 

LINE 

Printout  counter 

MAIN 

LPHI 

Number  of  flow  field  leaves 

MAIN 

LQ 

Index 

MAIN 

M 

Index 

MAIN,  SLICE 

MACH 

Mach  number 

SLICE 

MAXNOY 

Maximum  noy  value 

PNLC 

MC 

Convection  Mach  number 

SLICE,  SHOCK,  CRD 

MCIN 

Input  array  of  Mc 

SLICE 

MIN 

Input  array  of  M0 

CRD 

MJ 

Jet  exit  Mach  number 

SHOCK 

N 

Index,  also  number  of  shocks 

MAIN,  SHOCK,  LSPFIT 

NBREF 

Reference  boundary  number 

MAIN 

NCASE 

Number  of  cases 

MAIN 

NCBDY 

Number  of  centerbody  points 

MAIN 

NCELL 

Number  of  shock  cells 

MAIN,  SHOCK 

N COUNT 

Counter 

LSPFIT 

NN 

Acoustic  calculation  selector 

MAIN,  SLICE 

NODE 

Intermediate  variable 

MAIN 

NOV 

Minimum  number  of  points 

MAIN 

NOY 

Noy  value 

PNLC 

NPAGE 

Page  counter 

MAIN 

NPR 

Printout  counter 

MAIN 

NPRINT 

Printout  selector 

MAIN,  SLICE 

NPTS 

Number  of  points 

LSPFIT 

NR 

Number  of  points 

SLICE,  CRD 

NR1 

Index 

SLICE 

NTP 

Number  of  turning  points 

SLICE,  CRD 

NUM 

Number  of  boundary  points 

MAIN 

NUMANG 

Arena  selector 

MAIN,  SLICE 

NUMK 

Number  of  boundary  points 

MAIN 

NXC 

Index 

LSPFIT 

OAPWL 

Overall  power  level 

OUTPUT 

OASPL 

Overall  sound  pressure  level 

OUTPUT,  PNLC 

OBSTN 

Observed  Strouhal  number 

OUTPUT 

OMEGR 

Source  radian  frequency 

SLICE 

PAA 

Ambient  static  pressure 

MAIN 

PC 

Intermediate  variable 

PNLC 

PGC 

Gas  constant  parameter 

MAIN 

PHI 

Angle  <|> 

MAIN 

PI 

Constant  tt 

MAIN,  SLICE,  OUTPUT 

P 102 

tt/2 

CRD 

PI2 

2tt 

MAIN 
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Table  2.  List  of  FORTRAN  Symbols 

(Continued) . 

FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

PNDB 

PNdB 

PNLC 

PNL 

PNL 

OUTPUT,  PNLC 

PNLT 

PNLt,  tone-corrected  PNL 

OUTPUT 

POWER 

Exponent 

MAIN 

PS 

Ambient  static  pressure 

MAIN,  SHOCK 

PSQ 

Square  of  acoustic  pressure 

OUTPUT 

PSQM 

Mixing  noise  p2 

Shock  noise  p? 

SHOCK 

PSQS 

SHOCK 

PSQT 

Total  noise  p^ 

SHOCK 

PT 

Stagnation  pressure 

MAIN,  SHOCK 

PWL 

Power  level 

OUTPUT 

PWR 

Sound  power,  watts 

OUTPUT 

Q 

Intermediate  variable 

MAIN 

RA 

Boundary  coordinate  radius 

MAIN 

RAD 

Flow  integration  variable  R0 

MAIN 

RADO 

Flow  integration  variable  R0 

MAIN 

RADIUS 

Nozzle-to-observer  radius  R 

SLICE,  OUTPUT,  ATMOS 

RADX 

Argument  R0/Cmx 

MAIN 

RCBDY 

Centerbody  radial  cooordinate 

MAIN 

PRCRIT 

Critical  pressure  ratio 

SHOCK 

RCRC 

Intermediate  variable 

MAIN 

RFO 

Intermediate  variable 

OUTPUT 

RHO 

Density  p 

MAIN 

RHOE 

Ambient  density  pa 

MAIN,  SLICE 

RHOESQ 

ea 

SLICE 

RHOR 

Azimuthally-averaged  p 

MAIN,  SLICE 

RIN 

Input  radius 

SLICE,  CRD 

RJET 

Reference  jet  density  ratio 

MAIN 

RMIN 

Minimum  value  of  r 

MAIN 

RMINEX 

Exit  plane  value  of  RMIN 

MAIN 

RMINSQ 

Square  of  RMIN 

MAIN 

RMNSQE 

Square  of  RMINEX 

MAIN 

RMP 

Dummy  variable 

MAIN 

RND 

Normalized  radius  r/Deq 

MAIN 

ROOT 

Root  (zero)  of 

SLICE 

R00T2 

Jl 

SLICE 

RO 

Source  radius  ro 

CRD 

RSIG 

Turning  point  radius  ra 

SLICE,  CRD 

RSIG1 

ral 

CRD 

RSIG2 

Source  location  correction  (R*/R)2 

CRD 

RSORSQ 

SLICE 

RSTAR 

Source-to-observer  radius  R* 

SLICE 
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Table  2.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

RU 

Mass  flux  pU 

MAIN 

RU2 

Momentum  flux  pu2 

MAIN 

RU2E 

Exit  plane  value  of  pU^ 

MAIN 

RU2M 

Minimum  value  of  pu2 

MAIN 

RU2REF 

Reference  value  of  pU2 

MAIN 

R1 

Intermediate  variable 

CRD 

R2 

Intermediate  variable 

CRD 

S 

Intermediate  variable 

TPNLC 

SA 

Intermediate  variable 

MAIN 

SAC 

Intermediate  value  of  x^ 

MAIN 

SACO 

Intermediate  value  of  x^ 

MAIN 

SAR 

Intermediate  value  of  xr 

MAIN 

SARO 

Intermediate  value  of  xr 

MAIN 

SAX 

Intermediate  value  of  xx 

MAIN 

SAXO 

Intermediate  value  of  xx 

MAIN 

SBAR 

Intermediate  variable 

TPNLC 

SDU 

Intermediate  value  of  9U/3r 

MAIN 

SEFE 

Integral  of  enthalpy  flux 

MAIN 

SGN 

Sign 

LSPFIT 

SGN1 

Sign 

CRD 

SGN2 

Sign 

CRD 

SGI 

Intermediate  variable 

CRD 

SG2 

Intermediate  variable 

CRD 

SHIELD 

Shielding  integral 

SLICE,  CRD 

SIC 

Intermediate  value  of  x^ 

MAIN 

SIG 

Radius  r 

MAIN 

SIGN 

Sign 

ERF 

SIGSQ 

r^ 

MAIN 

SIGR 

Radius  r 

MAIN,  SLICE 

SINT 

Sin  0 

MAIN 

SINTO 

Sin  0O 

MAIN 

SIR 

Intermediate  value  of  xr 

MAIN 

SIX 

Intermediate  value  of  xx 

MAIN 

SPL 

SPL  array 

ALL 

SPLL 

Intermediate  variable 

TPNLC 

SPLMAX 

Maximum  SPL 

SHOCK 

SPLP 

Intermediate  variable 

TPNLC 

SPLPP 

Intermediate  variable 

TPNLC 

SPLU 

Intermediate  variable 

TPNLC 

SRU 

Mass  flux  integral 

MAIN 

SRUM 

Mass  flux  integral 

MAIN 

SRU2 

Momentum  flux  integral 

MAIN 

SRU2M 

Momentum  flux  integral 

MAIN 

Table  2.  List  of  FORTRAN  Symbols  (Continued). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

SS 

SPL  array 

PNLC 

SSPL 

Shock  noise  SPL  array 

SHOCK 

STC 

Azimuthal  shear  stress  T<f, 

MAIN 

STR 

Radial  shear  stress  xr 

MAIN 

STRFR 

Radial  coordinate  stretching  factor 

MAIN 

STRFX 

Axial  coordinate  stretching  factor 

MAIN 

STX 

Axial  shear  stress  xx 

MAIN 

SUE 

Reference  velocity 

MAIN 

SUM 

Sura 

OUTPUT 

SUMNOY 

Sum  of  noy  value 

PNLC 

SUMSPL 

Sum  of  SPL 

PNLC 

SUM1 

Sum 

CRD 

SUM2 

Sura 

CRD 

SU8 

Integral  of  source  strength 

MAIN 

SU8M 

Integral  of  source  strength 

MAIN 

SV2 

Square  of  velocity 

MAIN 

SI 

Intermediate  variable 

LSPFIT 

T 

Temperature 

ERF,  MAIN 

TA 

Intermediate  variable 

MAIN 

TAA 

Ambient  static  temperature 

MAIN 

TAO 

Intermediate  variable 

MAIN 

TAU 

Total  shear  stress  x 

MAIN 

TAUR 

Azimuthal  average  of  x 

MAIN,  SLICE 

TC2 

Intermediate  variable 

TPNLC 

TC3 

Intermediate  variable 

TPNLC 

TE 

Exit  static  temperature 

MAIN 

TEMP 

Normalized  temperature  T/Ta 

SLICE 

TERM 

Directivity  factor 

SLICE 

TH 

Angle  <j> 

MAIN 

THCR 

Critical  angle  9cr 

SHOCK 

TERM 

Directivity  factor 

SLICE 

THE 

Angle  0 

SLICE,  CRD 

THETA 

Observer  angle  9,  radians 

SLICE,  OUTPUT 

THETD 

Observer  angle  0j>  degrees 

SLICE,  OUTPUT,  SHOCK 

THO 

$o 

MAIN 

THT 

Observer  angle  0i,  radians 

SHOCK 

TI 

Intermediate  value  of  enthalpy  flux 

MAIN 

TOP 

Intermediate  variable 

LSPFIT 

TSR 

Static  temperature 

MAIN 

TSTD 

Circumferential  asymmetry  test  parameter 

MAIN 

TSTH 

Circumferential  asymmetry  test  parameter 

MAIN 

TSTHL 

Circumferential  asymmetry  test  parameter 

MAIN 

TSTL 

Circumferential  asymmetry  test  parameter 

MAIN 

Table  2.  List  of  FORTRAN  Symbols  (Concluded). 


FORTRAN 

Symbol 

Meaning 

Related  Subroutines 

TT 

Stagnation  temperature 

MAIN 

TTR 

Stagnation  temperature 

MAIN 

TURBIN 

Turbulence  intensity,  u' 

MAIN 

U 

Mean  velocity 

MAIN 

UAP 

Intermediate  variable 

MAIN 

UAVG 

Mass-average  of  U  at  x 

MAIN 

UC 

Convection  velocity  Uc 

SHOCK 

UE 

Exit  plane  velocity  Uj 

MAIN,  SHOCK 

UGLY 

Intermediate  variable 

MAIN 

UJET 

Reference  exit  velocity 

MAIN 

UMAX 

Maximum  local  value  of  U  at  x 

MAIN 

UND 

Normalized  value  of  U,  U/UREF 

MAIN 

UNITS 

Constant  for  units  conversion 

MAIN,  OUTPUT 

UR 

Azimuthal  average  of  U 

MAIN,  SLICE 

UREF 

Reference  exit  velocity 

MAIN 

U8 

Intermediate  value  of  source  strength 

MAIN 

U8I 

Integral  of  source  strength 

MAIN 

VA 

Intermediate  value  of  momentum 

MAIN 

VAO 

Intermediate  value  of  momentum 

MAIN 

VI 

Intermediate  value  of  momentum 

MAIN 

VMAX 

Maximum  of  velocities  inside  and  outside 

MAIN 

VMIN 

Minimum  of  velocities  inside  and  outside 

MAIN 

VO 

Flight  velocity  Ua 

SHOCK 

VR 

Velocity  ratio  VMIN/VMAX 

MAIN 

WITHIN 

Dummy  variable 

LSPFIT 

X 

Axial  distance  x 

MAIN,  SLICE 

XCBDY 

Centerbody  axial  coordinate 

MAIN 

XD 

Intermediate  variable 

LSPFIT 

XE 

Exit  plane  axial  coordinate 

MAIN 

XND 

Normalized  axial  coordinate  X/Deq 

MAIN 

XOR 

Variable  x/R 

SLICE 

XI 

Intermediate  variable  for  curve  fitting 

LSPFIT 

X13 

Intermediate  variable  for  curve  fitting 

LSPFIT 

X4 

Intermediate  variable  for  curve  fitting 

LSPFIT 

XA3 

Intermediate  variable  for  curve  fitting 

LSPFIT 

Y 

Intermediate  variable  for  curve  fitting 

LSPFIT 

YC 

Intermediate  variable  for  curve  fitting 

LSPFIT 

YI 

Intermediate  variable  for  curve  fitting 

LSPFIT 

Y3 

Intermediate  variable  for  curve  fitting 

LSPFIT 
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(b)  Subroutine  SLICE 


Figure  8.  Computer  Program  Flowchart  (Continued). 
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4.3.4  LSPFIT 


Subroutine  LSPFIT  is  a  curve-fitting  routine  which  utilizes  least- 
squares  polynomial  fits  of  second  order  to  perform  interpolation,  differentia¬ 
tion  and  integration  of  input  discrete-point  data.  The  calling  statement  is: 


CALL  LSPFIT (X,  Y,  N,  XC,  YC,  NC,  NF,  A) 


where  (X,  Y)  are  the  input  data  coordinates  (N  values  of  each),  XC  are  the 
values  of  X  where  output  is  requested,  YC  are  the  output  functions,  NC  is  the 
number  of  output  data  points,  and  NF  indicates  the  type  of  output  desired. 

The  coding  for  NF  is  as  follows: 


NF  =  0, 

YC  are 

interpolated  values  of  Y 

NF  =  1, 

YC  are 

derivatives  of  Y 

NF  =  -1, 

YC  is 
XC(1) 

the  integral  of  Y  from 
to  XC(J) ,  1  1  J  £  NC. 

The  parameter  A  is  the  second  derivative  of  Y.  Subroutine  LSPFIT  is  used  in 
MAIN  to  interpolate  input  plug/centerbody  geometry  coordinates  at  various 
axial  stations  in  the  flow  field,  and  to  obtain  radial  gradients  of  density 
from  the  computed  density  profiles. 


4.3.5  SLICE 

Subroutine  SLICE  directs  the  mixing  noise  calculation  for  each  axial 
slice.  The  calling  sequence  is  as  follows: 

CALL  SLICE  (X (KA) ,  DSIG(KA),  DX,  M) 

where  X(KA)  is  the  axial  location,  DSIG(KA)  is  the  radial  step  size,  DX  is 
the  axial  slice  thickness,  and  M  is  the  number  of  radial  points  in  the  slice. 
The  flow  parameters  (which  are  circumferentially  mass-averaged  values)  are 
transferred  through  labeled  COMMON  statements.  Subroutine  SLICE  computes  the 
acoustic  arena  geometry  parameters  THETA,  THETD,  RADIUS  and  initializes  SPL 
(I,J)  to  zero  during  the  first  call,  skipping  this  calculation  on  succeeding 
calls.  The  normalized  radial  profiles  of  velocity  (MACH)  and  temperature 
(TEMP)  are  evaluated,  followed  by  a  calculation  of  source  strength  amplitude 
DS  and  characteristic  frequency  FS  for  each  radial  volume  element. 
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Subroutine  SLICE  computes  the  acoustic  shielding  function  profiles  G2(J), 
the  number  of  turning  points  NTS’,  and  their  locations  RS1G.  Subroutine  CRD 
is  then  called  to  calculate  the  acoustic  shielding  exponentials  and  quadrupole 
directivity  functions.  Subroutine  SLICE  then  sums  up  the  mixing  noise  con¬ 
tributions  from  each  radial  volume  element,  factoring  in  their  individual 
source  strengths,  characteristic  frequencies,  spectrum  shapes,  directivities, 
and  shielding  factors.  The  resulting  noise  spectra  from  each  slice  is  stored 
as  the  variable  DSPL(I,J),  where  I  denotes  the  observer  angle  index  and  J  is 
the  1/3-octave  frequency  band  index.  Upon  completing  the  calculation  for  a 
given  slice,  SLICE  returns  controL  to  MAIN. 

4.  1.6  CRD 

Subroutine  CRD  computes  the  shielding  function  integrals  and  quadrupole 
directivity  factors  for  a  given  axial  slice  as  a  function  of  radial  source 
location.  The  radial  distributions  of  normalized  velocity  (MACH)  and 
temperatures  (TEMP)  and  shielding  function  (G2)  are  transferred  to  CRD  through 
labeled  COMMON  statements.  The  calling  statement  is: 

CALI.  CRD 


At  each  source  radius,  subroutine  CRD  interrogates  the  data  to  determine 
which  of  the  six  shielding  conditions  in  Figure  4  applies,  and  computes  the 
appropriate  shielding  integral  (,'oi  .  *‘02<  or  l512^>  an<*  the  appropriate 
directivity  factors.  Alter  all  radial  source  volumes  have  been  evaluated, 
CRD  returns  control  to  SLICE. 


4.3.7  OUTPUT 


Subroutine  OUTPUT  performs  the  final 
out  the  far  field  SPL  spectra,  OASPL,  PNL 
sequence  is  as  follows: 


acoustic  calculations  and  prints 
and  PNLT  directivities.  The  calling 


CALL  OUTPUT  (EMACH ,  D.1ET,  RJET,  UJET,  UNITS) 


where  EMACH,  DIET,  RJET,  and  UJET  are  the  characteristic  (usually  reference) 
jet  Mach  rumber,  diameter,  density  ratio  and  velocity,  respectively.  The 
parameter  UNITS  is  a  conversion  factor  for  converting  from  lbf/ft‘-  to 
dynes/cm-  relative  to  0.0002  dvnes/ctn^.  Subroutine  OUTPUT  converts  the 
narrow-hand  spectra  from  SLICE  into  1/3-octave  levels.  Subroutine  SLICE  then 
calls  SHOCK  to  compute  SSPL  spectra  (shock  noise)  and  adds  these  to  the 
turbulent  mixing  noise  spectra  to  obtain  the  total  noise  spectra.  The 
corresponding  power  spectrum  (PWL)  is  then  computed,  and  subroutine  ATMOS  is 
then  called  to  correct  all  SPL  spectra  for  atmospheric  attenuation.  Sub¬ 
routines  PNLC  and  TPNLC  are  then  called  to  calculate  perceived  noise  level 
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PNL  and  tone-corrected  noise  level  PNLT.  Finally,  overall  sound  pressure 
level  OASPL  is  computed,  and  all  of  these  acoustic  parameters  are  then 
printed  out.  Subroutine  OUTPUT  then  returns  control  to  MAIN. 


4.3.8  SHOCK 

Subroutine  SHOCK  computes  the  broadband  shock-associated  noise  spectra 
at  each  observer  angle.  The  calling  statement  is  as  follows: 


CALL  SHOCK 


All  parameters  are  transferred  into  and  out  of  this  subroutine  through 
labeled  COMMON  statements.  Subroutine  SHOCK  computes  the  1/3-octave  SPL 
spectra  for  each  nozzle  boundary  element  which  has  a  non-zero  shock  cell 
number  input,  NCELL  >  0.  The  individual  boundary  contributions  are  summed  on 
a  mean-square  pressure  basis  and  added  to  the  mixing  noise  spectra. 


4.3.9  ATMOS 

Subroutine  ATMOS  corrects  the  input  SPL  spectra  for  atmospheric  attenu¬ 
ation  effects  using  standard-day  atmospheric  absorption  factors  for  70% 
relative  humidity  and  59°  F  ambient  conditions.  The  calling  sequence  is  as 
follows: 


CALL  ATMOS  (SPL,  RADIUS) 


where  SPL(I,J)  is  the  sound  pressure  spectrum  array,  I  denotes  the  index  for 
observer  angle,  J  denotes  the  index  on  frequency,  and  RADIUS(I)  is  the 
nozzle-to-observer  distance  array.  The  atmospheric  absorption  in  dB  per  1000 
ft  from  reference  (15)  is  corrected  to  the  proper  distance  RADIUS (J),  and  the 
result  is  subtracted  from  SPL(I,J).  The  array  of  SPL(I,J)  returned  to  OUTPUT 
is  the  corrected  array. 


4.3.10  PNLC 

Subroutine  PNLC  computes  the  perceived  noise  level  in  PNdB  at  each 
observer  angle  from  the  input  1/3-octave  spectra.  The  calling  sequence  is  as 
follows : 


CALL  PNLC  (SS,  FAC,  PNDB,  OASPL) 
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where  SS  is  the  input  array  of  either  1/ 3-oc.tave  or  octave  SPL  values,  FAC  is 
a  constant  equal  to  0.15  for  1/3-octave  and  0.3  for  octave  levels,  PNDB  is 
the  output  PNL,  and  OASPL  is  the  conventional  overall  level.  The  method  used 
to  calculate  PNL  is  taken  from  reference  (17).  The  OASPL  output  from  sub¬ 
routine  PNLC  is  discarded  because  it  only  computes  the  summation  for  the 
first  24  values  of  SS.  This  is  sometimes  insufficient  for  model  scale 
conditions  where  the  frequency  range  of  interest  can  cover  as  many  as  33  1/3- 
octave  frequency  bands. 


4.3.11  TPNLC 


Subroutine  TPNLC  determines  a  pure-tone  correction  factor  to  the  PNL 
value  as  a  function  of  the  1/3-octave  SPL  spectrum.  The  calling  sequence  is 
as  follows: 


CALL  TPNLC  (SPL,  PTCOR) 


where  SPL  is  the  input  1/3-octave  spectrum  and  PTCOR  is  the  correction  to  be 
applied  to  PNL  to  account  for  the  presence  of  tones  in  the  spectrum.  Sub¬ 
routine  TPNLC  reads  in  SPL  and  returns  PTCOR.  The  tone  correction  and 
detection  procedure  is  based  on  the  method  proposed  in  reference  (18). 


4.4  PROGRAM  USAGE  AND  LOGIC 


A  complete  description  of  the  program  input  variables  and  input  format 
is  given  in  Appendix  A.  A  list  of  notes  and  suggestions  on  running  the 
program  is  also  included.  A  description  of  program  output  format,  including 
warning  flags  and  diagnostics,  is  given  in  Appendix  B.  A  sample  case  listing 
(including  input  data  card  images)  is  given  in  Appendix  C  for  a  7-tube 
suppressor  nozzle,  one  of  the  data-theory  comparison  cases  presented  in 
Reference  (3).  A  complete  FORTRAN  source  listing  of  the  program  logic  is 
given  in  Appendix  D. 

Program  users  should  be  completely  familiar  with  Appendix  A,  since  there 
are  many  pitfalls  which  can  be  avoided  by  giving  attention  to  the  recom¬ 
mendations  presented  therein.  The  program  flexibility  permits  analysis  of 
nozzle  planforms  of  any  imaginable  shape,  so  long  as  certain  input  rules  and 
guidelines  are  followed.  When  non-axisymmetric  nozzles  are  run,  a  completely 
three-dimensional,  turbulent,  compressible  flow  field  analysis  is  performed, 
and  input  mistakes  can  be  costly  in  terms  of  computer  processor  time.  The 
user  should  make  initial  checkout  runs  for  complex  nozzles,  running  just  one 
or  two  axial  slices  at  first,  to  ensure  that  all  input  is  as  desired,  before 
running  a  complete  jet  plume. 

The  program  is  designed  to  serve  as  a  diagnostic  tool,  in  addition  to 
functioning  in  the  standard  jet  noise  prediction  mode.  Individual  slice 
calculations  can  be  made  by  suitable  input  selection,  running  each  slice  (or 
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5.0  CONCLUDING  REMARKS 


A  theoretical  analysis  and  associated  computer  mode]  for  the  solution  of 
the  flow  field  and  the  acoustic  characteristics  of  nozzles  of  arbitrary  geo¬ 
metry  is  presented.  A  large  number  of  cases  have  been  calculated  with  this 
model  and  compared  with  appropriate  experimental  data  (Reference  3).  The 
computerized  procedure  presented  herein  provides  a  reasonably  accurate  method 
of  predicting  the  aeroacoustic  characteristics  of  a  wide  variety  of  exhaust 
nozzles  over  the  range  of  flow  conditions  of  interest,  at  least  for  those 
conditions  and  or  observer  angles  where  shock-cell  noise  does  not  dominate  the 
spectrum. 


APPENDIX  A 


DESCRIPTION  OF  INPUT 


The  input  data  is  supplied  through  NAMELIST  input  format,  with  the 
exception  of  the  alphanumeric  title  data  card,  which  precedes  the  input 
NAMELIST  data.  Any  number  of  successive  cases  can  be  run  consecutively, 
limited  only  by  the  user's  execution  time  available.  Each  successive  case 
requires  a  title  card  (80  -  character  label  in  columns  1  -  80) ,  followed  by 
the  INPUT  NAMELIST.  The  data  from  preceding  cases  remains  in  storage,  so 
only  those  variables  which  are  to  be  changed  from  the  preceding  case  input 
value  need  be  included  in  the  INPUT  file  of  succeeding  cases. 

A  suggested  input  preparation  format  is  given  in  Table  A-l.  Those 
variables  marked  by  an  asterisk  (*)  have  preset  values  built  into  the  program, 
and  need  not  be  input  unless  the  user  desires  to  over-ride  the  preset  value 
with  a  different  one.  The  definitions  of  each  of  the  input  variables  given  in 
Table  A-l  are  listed  in  Table  A-2.  Again,  preset  variables  are  marked  by  an 
asterisk  (*).  The  values  of  those  variables  which  are  preset  are  given  in 
Table  A-3.  The  format  of  Table  A-2  is  such  that  a  note  number  (where  approp¬ 
riate)  is  given  for  each  variable  which  corresponds  to  the  note  number  in 
Section  A. 1  (Notes  on  Input).  These  notes  give  further  elaboration  on  how  to 
specify  and  prepare  the  input  data. 
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Table  A-l.  Suggested  Input  Format 


Column 

2 

(80  -  CHARACTER  TITLE  CARD,  COLUMNS  1-80) 

$INPUT 

KX*  a  _ ,  NEST  =  _ ,  LPHI  *  =  _ ,  ISYM  =  _ 

I  QUIT*  =  _ ,  NN*  =  _ ,  NCASE*  =  _ ,  NBREF  * 

NPRINT*  =  _ ,  NCBDY  =  _ , 

N0V  = 

X  =  _ 

DSIG  =  _ 

BETAIN*  = 

DELTIN*  = 

AMUIN*  = 

RMIN 


XE  =  0 

ALP0  =  _ 0 

LEAV  =  1 

NUM  =  1 

KN  1 


DEQ  = 
DS 

NCELL  = 


Table  A-l.  Suggested  Input  Format  (Concluded). 


Column 

2 

DALP(1 ,2)  = 

y 

y  y 

y  y 

DALP (1 , 3)  = 

y 

y  y 

y  y 

(etc.,  for  boundary  4,  5,  6, 

.  . .NEST) 

RA(1,  2)  = 

y 

y  y 

y  y 

RA(1,  3)  = 

y 

y  y 

y  y 

(etc.,  for  boundary  4,  5,  6,  .. 

.  .  .Nest) 

CM*  = 

,  CH*  = 

,  CMVR*  = 

,  CMMC*  = 

y 

GAM  = 

,  CP  = 

,  PS  = 

,  ALFA*  = 

y 

DTHM  *  = 

,  RU2M  *  = 

,  AK*  = 

,  BK*  = 

y 

STRFR*  = 

,  STRFX*  = 

,  AT0TAL  = 

y 

ALPHMC*  =  _ 

,  BETAMC  *  = 

y 

NUMANG  = 

,  DIST  = 

,  FMIN  *  = 

,  FMAX*  = 

ALPHT  *  = 

y  y 

y 

y  y 

XCBDY  = 

y  y 

y 

y  y 

y 

RCBDY  = 

y  y 

y 

y  y 

y 

$ 

(NEXT  CASE,  IF  ANY) 
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Table  A-2.  Input  Variable  Definitions. 


Variable 

KX* 

NEST 

LPHI 

ISYM 

IQUIT 

NN* 

NCASE* 

NBREF* 

NPRINT* 

NCBDY 

N0V 

X 

DSIG 

BETAIN* 

DELTIN* 


Note  Description 

Number  of  axial  stations  to  be  analyzed;  a 
maximum  of  24  stations  is  permitted. 

I  Number  of  closed  boundary  contours  defining 
the  nozzle  exit  geometry.  A  maximum  of  110 
is  permitted. 

7  Number  of  symmetric  leaves  (repeating  seg¬ 

ments  in  the  nozzle  exit  planform. 

Nozzle  symmetry  indicator;  ISYM  =  1  for  ax- 
symmetric  nozzles  or  completely  asymmetric 
nozzles,  =  0  otherwise. 

Maximum  number  of  radii  at  which  flow  field 
is  calculated  (^200). 

12  Acoustic  Calculation  option  indicator. 

Number  of  cases  to  be  run  consecutively. 
Reference  condition  boundary  number. 

13  Aerodynamic  station  printout  indicator. 

9  Number  of  centerbody  input  coordinate  points. 

A  maximum  of  40  is  permitted. 

Minimum  number  of  radii  at  which  flow  field 
is  to  be  calculated,  for  each  axial  station 
(KX  values  required). 

II  Axial  location  of  each  axial  station,  ft. 

(KX  values  required). 

11  Radial  step  size  to  be  taken  for  flow  field 

calculation  at  each  axial  station,  ft. 

(KX  values  required). 

15  Axial  shear  stress  turbulence  constant 

(KX  values  required). 

Azimuthal  shear  stress  turbulence  constant 
(KX  values  required). 
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Table  A-2.  Input  Variable  Definitions.  (Continued) 


Variable 

Note 

Description 

AMUIN* 

Azimuthal  velocity  gradient  turbulence  fre¬ 
quency  constant  (KX  values  required). 

RMIN 

9 

Minimum  radius  for  flow  field  calculation 
at  each  axial  station  (KC  values  required). 

XE 

8 

Axial  location  of  exit  plane  of  each 
boundary,  ft.  (NEST  values  required). 

ALP0 

2 

Reference  angle  aQ  from  which  the  coordi¬ 
nates  of  each  boundary  point  are  specified, 
radians  (NEST  values  required). 

LEAV 

1,4 

Number  of  symmetric  leaves  (repeating  seg¬ 
ments)  of  each  boundary  (NEST  values 
required) . 

NUM 

1,5 

Number  of  input  points  (coordinate  pairs) 
to  be  supplied  for  each  boundary  (NEST 
values  required) . 

KN 

1 

The  number  of  the  boundary  which  encloses 
a  given  boundary  (NEST  values  required) . 

DEQ 

16 

Equivalent  flow  area  diameter  of  each 
boundary,  ft.  (NEST  values  required). 

DS 

16 

Shock-cell  spacing  characteristic  dimension 
usually  hydraulic  diameter,  of  each  boun¬ 
dary,  ft.  (NEST  values  required). 

NCELL 

16 

Number  of  shock  cells  for  each  boundary 
element  (NEST  values  required). 

PT 

6 

Stagnation  pressure  Inside  each  boundary, 
lbf/ft  (NEST  values  required). 

TT 

6 

Stagnation  Temperature  inside  each  boun¬ 
dary  °R  (NEST  values  required). 

DALP ( I ,  J  ) 

2,3,5 

Angular  increment  Act  from  preceding  boun¬ 
dary  point  which  locates  the  given  boun¬ 
dary  point  I  on  boundary  J,  radians  (omit 
boundary  number  1,  ambient  field). 
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Table  A-2.  Input  Variable  Definitions  (Continued). 


Variable 

Note 

Description 

RA  (I,J) 

2,3,5 

Radial  coordinates  of  boundary  point  I  on 
boundary  J,  ft.  (omit  boundary  number  1, 
ambient  f ield) . 

CM* 

10 

Empirical  jet  momentum  diffusion  rate 
spreading  parameter  Cm. 

CH* 

10 

Ratio  of  enthalpy- to-momentum  spreading 
parameters  Ch/Cm. 

CMVR* 

10 

Momentum  spreading  parameter  velocity 
ratio  influence  coefficient. 

erne* 

10 

Momentum  spreading  parameter  Mach  number 
influence  coefficient. 

GAM 

Specific  heat  ratio  y  =  Cp/Cv. 

CP 

Specific  heat  at  constant  pressure  Cp, 
in  (f t-lbf ) / (slug  -  °R) . 

PS 

Ambient  static  pressure,  lbf/ft  . 

ALFA* 

Turbulence  characteristic  frequency  con¬ 
stant. 

DTHM  * 

7 

Maximum  allowable  increment  in  angular 
coordinate,  (d<J>)max>  for  flow  field  calcula¬ 
tion. 

RU2M* 

2 

Minimum  value  of  jet  momentum  flux,  (pU  )  .  , 
below  which  the  flow  is  not  calculated. 

AK* 

Sound  pressure  level  proportionality  con¬ 
stant  for  mixing  noise  calculation. 

BK* 

Sound  pressure  level  proportionality  con- 
for  dipole  density-gradient  noise  calcula- 
t  ion. 

STRFR  * 

11 

Radial  coordinate  stretching  factor  for  use 
of  automatic  mesh  calculation. 

STRFX  * 

11 

Axial  coordinate  stretching  factor  for  use 
of  automatic  mesh  calculation. 
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Table  A-2.  Input  Variable  Definitions  (Concluded). 


Variable 

Note 

Description 

AT0TAL 

2 

Nozzle  Total  exit  flow  area,  ft  . 

ALPHMC * 

14 

Convection  Mach  number  weighting  factor. 

BETAMC* 

14 

Convection  Mach  number  weighting  factor. 

NUMANG 

Arena  selection  indicator;  NUMANG  =  1  indi¬ 
cates  constant  radius  arc,  NUMANG  =  2 
indicates  sideline  parallel  to  the  jet 
axis. 

DIST 

Arc  or  sideline  distance,  ft. 

FMIN  * 

Minimum  frequency  for  which  acoustic  cal¬ 
culations  are  required,  Hz  (£50) ;  an 
integer  variable. 

FMAX  * 

Maximum  frequency  for  which  acoustic  calcu¬ 
lations  are  required,  Hz  (£100,000);  an 
integer  variable. 

ALPHT* 

Convective  amplification  factor  turbulence 
constant  at;  15  values  required,  one  for 
each  observer  angle  9j  from  0j  =  20°  to  160 
in  10°  increments. 

XCBDY 

9 

Centerbody  input  point  axial  coordinate, 
NCBDY  values  required. 

RCBDY 

9 

Centerbody  input  point  radial  coordinate, 
NCBDT  values  required. 
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Table  A-3.  Preset  Input  Values. 


Variable 

Value 

AK 

0.08 

ALFA 

1.0 

ALPHT 

15*  0.5 

ALPHMC 

0.5 

AMU  IN 

24*  0.2 

BETAIN 

24*  4.0 

BETAMC 

0.325 

BK 

0.0 

CH 

1.15 

CM 

0.075 

CMMC 

0.08 

CMVR 

0.25 

DELTIN 

24*  4.0 

DTHM 

0.1 

FMAX 

100000 

FMIN 

50 

I  QUIT 

50 

KX 

15 

LPHI 

9999 

NBREF 

2 

NCASE 

1 

NN 

0 

NPRINT 

1 

RU2M 

3.0 

STRFR 

0.01 

STRFX 

1.259921 
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A. 1  NOTES  ON  INPUT 


1.  The  jet  nozzle  geometry  is  specified  by  input  of  the  number  of 
component  boundaries,  NEST,  along  with  pairs  of  coordinates,  RA  and 
DALP,  for  each  boundary  element.  The  ambient  field  is  always 
treated  as  the  first  boundary  in  the  input  arrays  for  UE,  PT,  TT, 
LEAV,  NUM,  KN,  XE,  and  ALPO.  This  is  why  some  numbers  have  already 
been  filled  in  on  Table  A-l  in  the  first  column  for  these  arrays. 

A  nozzle  with  N  elements  has  NEST  =  N  +  1  boundaries. 

2.  The  steps  to  specifying  nozzle  geometry  input  are  as  follows,  re¬ 
ferring  to  Figure  6: 

a.  Obtain  sketch  or  drawing  of  nozzle  exit  cross  section  and 
select  a  coordinate  origin  which  is  optimum  from  the  stand¬ 
point  of  symmetry  and  boundary  point  specification, 

b.  Number  each  boundary,  reserving  boundary  Number  1  for  the 
ambient  field. 

c.  With  respect  to  the  coordinate  origin,  select  a  reference 
angular  location  for  each  boundary,  ALP0. 

d.  For  each  boundary,  select  points  represented  by  pairs  of 
coordinates.  The  coordinates  used  as  input  are  radius, 

RA(I,J),  and  angular  increment  from  the  preceding  point, 

DALP (I, J).  For  the  first  point,  DALP (I, J)  is  the  angular 
increment  from  the  reference  angle  ALP0.  The  index  I  is  the 
boundary  point  number,  and  the  index  J  is  the  boundary  number. 
Both  ALP0  and  DALP  are  to  be  input  in  radians,  and  RA  is  input 
in  feet. 

3.  The  last  point  on  a  given  boundary  should  be  located  at  ALP0  if 
the  boundary  has  only  one  leaf.  The  sum  of  all  DALP (I, J)  should 
equal  zero  if  the  boundary  has  only  one  leaf. 

4.  If  the  boundary  is  a  circle  about  the  origin,  only  one  point  on  the 
boundary  need  be  supplied,  and  the  value  of  LEAV  for  that  boundary 
is  set  equal  to  the  number  of  boundary  points  desired  on  the  circle. 

5.  The  program  uses  linear  interpolation  between  input  boundary  point?.. 
If  a  boundary  is  made  up  of  or  contains  straight  line  segments, 
only  the  end-points  of  the  straight  line  segments  need  bo  input. 

6.  The  variables  PT  and  TT  refer  to  stagnation  pressure  and  temperature 
at  the  exit  plane  inside  the  boundary  of  interest.  Setting  the 
first  value  of  PT  equal  to  PS  gives  a  static  ambient  field.  The 
first  value  of  PT  greater  than  PS  simulates  non-zero  flight  veloc¬ 
ity. 
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7. 


The  variable  LPHI  determines  what  angular  extent  of  the  flow  field 
needs  to  be  calculated.  If  the  nozzle  geometry  is  axisymmetric, 
setting  LPHI  equal  to  a  large  number  (such  that  2tt/LPHI  is  less 
than  DTHM)  forces  the  program  to  calculate  the  flow  field  at  only 
one  angular  location.  The  flow  field  for  a  nozzle  containing  two 
adjacent  circular  jets,  for  example,  has  LPHI  =  4,  since  the  flow 
is  the  same  each  quadrant.  Several  examples  of  how  boundary  para¬ 
meters  are  specified  are  shown  in  Figure  A-l. 

8.  The  program  can  currently  only  handle  coplannar  nozzles,  i.e., 
every  nozzle  element  must  terminate  at  the  same  axial  location. 
Therefore  XE  must  be  the  same  for  all  input  boundaries. 

9.  The  centerbody,  if  any,  is  input  through  coordinates  pairs  XCBDY(J), 
RCBDY(J),  where  1  <_  J  <_  NCBDY.  A  maximum  of  40  points  can  be 
input.  The  LSPFIT  subroutine  uses  this  input  to  interpolate  for 
finding  the  values  of  RMIN  at  each  axial  location  X.  The  LSPFIT 
routine  can  treat  line  segments,  both  straight  and  curved.  Typical 
examples  of  centerbody  coordinate  input  are  shown  in  Figure  A-2. 

If  there  is  no  centerbody,  the  user  can  avoid  automatic  computation 
of  the  potential  core  of  axisymmetric  nozzles  (which  has  no  impact 
on  mixing  noise)  by  specifying  RMIN  as  input,  but  with  NCBDY  =  0. 
This  option  causes  the  computation  to  begin  at  r  =  RMIN(KA),  where 
KA  is  the  axial  station  number. 

10.  The  input  value  of  CM  is  modified  for  velocity  ratio  and  Mach 
number  effects  by  the  re’ation 


(1  +  CMVR*VR) (1  +  CMMC*ACH) 

where  DBDX  is  the  modified  value  of  Cm,  and  VR  and  ACH  are  the 
velocity  ratio  and  Mach  number,  respectively,  of  a  given  boundary. 
The  heat  transport  spreading  parameter  is  then  calculated  from  the 
relation 


Ch  =  CH  *  DBDX 

The  values  of  CM,  CMMC,  CMVR  and  CH  recommended  and  preset  in  the 
program  are  given  in  Table  A-3.  These  values  can  be  changed  by  the 
user  to  reflect  experimental  evidence  if  so  desired, 

11.  The  axial  locations  of  the  axial  stations  can  be  input  by  the  array 
X(KA),  where  1  £  KA  £  KX.  The  radial  mesh  step  size  can  also  be 
input  by  the  array  DSIC(KA).  An  automatic  grid  selection  procedure 
has  been  devised  to  obviate  the  need  for  supplying  all  values  of 
X(KA)  and  DSIG(KA).  The  only  input  required  is  the  first  axial 
station  X(l),  and  the  grid  stretching  factors  STRFR  and  STRFX.  The 
grid  is  then  calculated  from  the  following  relations: 
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Point  Number 


Curved  Line  Segment 


Straight  Line  Segment 


Repeat  Coordinate 

NCBDY  =  7, 

XCBDY  X1,X2,X3,X4)X5)X6(-X5) ,X7, 

RCBDY  =  R1,R2,R3,R4, 0,0,0, 


(a)  Example  1  -  Curved  Centerbody 


NCBDY  =  10, 

XCBDY  =X1,X2,X3  (=X2) ,X,1,X5(=X4) ,X6,X7,X8,X9(=X8) ,X10, 
RCBDY  =  R1,R2,R3(  R2)  .K4.K5(:=R.1) -R6'R7'R8-R9(=R8)  ,RlO, 


(b)  Example  2  -  Segmented-Cone  Centerbody  with  Curved  Tip 


Figure  A-2.  Centerbody  Input  Coordinate  Examples. 
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X (KA)  =  STRFX  *  X(KA-l) 

DSIG(KA)  =  STRFR  *  X(KA) 

This  provides  a  grid  which  exhibits  larger  and  larger  step  sizes  as 
the  plume  is  developed  downstream.  Recommended  value  of  STRFR  and 
STRFX  are  preset  and  listed  in  Table  A-3. 

12.  The  variable  NN  determines  the  type  of  acoustic  calculation  desired. 
Normal  (preset)  operation  is  with  NN  =  0,  which  give  the  complete 
acoustic  calculation  described  in  Section  3.  The  user  may  desire 

to  perform  diagnostic  computations  to  assess  the  relative  impor¬ 
tance  of  convectioi. .  shielding,  etc.  By  selecting  the  appropriate 
value  of  NN,  the  various  components  of  the  acoustic  calculation  can 
be  switched  on  and  off  in  various  combinations.  Setting  NN  =  4 
gives  only  the  aerodynamic  calculation,  and  the  acoustic  calcula¬ 
tions  are  bypassed.  The  various  options  for  NN  are  listed  below: 

NN  =  0  -  complete  acoustic  calculation. 

NN  =  1  -  convective  amplif ication,  no  shielding. 

NN  =  2  -  no  convective  amplification,  no  shielding. 

NN  =  3  -  no  convective  amplification,  with  shielding. 

NN  =  4  -  no  acoustic  calculation,  aerodynamics  only. 

13.  The  printout  of  aerodynamic  flow  field  data  is  controlled  by  NPRINT. 
When  NPRINT  =  0,  no  aerodyanmic  printout  is  provided.  If  NPRINT  = 

1,  aerodynamic  printout  is  provided  at  every  axial  station.  If 
NPRINT  =  2,  aerodynamic  printout  is  provided  at  every  second  axial 
station  (i.e.,  KA  =  1,  3,  5,  7,  etc).  For  PRINT  =  3,  printout  is 
provided  at  every  third  station,  etc. 

14.  For  dual  flow  nozzles,  if  the  inner  stream  has  a  higher  velocity 
than  the  outer  stream,  use  ALPHMC  =0.5  and  BETAMC  =  0.325  (preset 
values).  These  variables  are  weighting  factors  in  the  convection 
Mach  number  calculation,  which  is  computed  from  the  relation 

MC  =  ALPHMC  *  MACH  +  BETAMC  *  EMACH 

where  MACH  is  the  local  acoustic  Mach  number  U/Ca  and  EMACH  is  the 
exit  plane  reference  Mach  number  Uj/Ca.  If  the  outer  stream  has  a 
higher  velocity  than  the  inner  stream,  use  ALPHMC  =0.5  and  BETAMC 
=  0.325/VR,  where  VR  =  (Uouter /Olnncr) j •  For  multielement  suppres¬ 
sor  nozzles,  VR  =  Uj/Um,  where  Um  is  the  post-merged  potential  core 
velocity.  Tf  Um  is  not  known,  BETAMC  =  0.2  to  0.25  is  usually  a 
good  approximation. 
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15.  For  dual  flow  nozzles,  input  BETAIN  =  4.0  (preset)  for  all  values 
of  X,  provided  the  inner  stream  velocity  is  higher  than  the  outer 
stream  velocity  at  the  exit  plane.  If  the  outer  stream  velocity  is 
higher  than  the  inner  stream  velocity  at  the  exit  plane,  input 
BETAIN  =  0  for  all  axial  stations  where  X  £  10*DEQ (NBREF) ,  and 
BETAIN  =  4.0  thereafter,  where  NBREF  is  the  outer  stream  boundary 
number.  For  multielement  nozzles,  input  BETAIN  =  0  for  axial 
distances  less  than  10*DEQ(1),  where  DEQ(l)  is  the  equivalent 
diameter  based  on  total  flow  area  at  the  exit  plane. 

16.  For  each  boundary  element  DEQ,  DS  and  NCELL  are  input.  The  first 
value  of  DEQ  is  the  total  flow  area  equivalent  diameter.  The  first 
value  of  NCELL  determines  whether  or  not  the  shock  cell  noise  is 
computed.  If  NCELL(l)  is  input  zero,  no  shock  noise  is  computed; 
for  NCELL(1)>0,  the  shock  cell  noise  routine  is  called.  The  shock 
noise  of  each  boundary  element  is  computed  separately  and  added  to 
the  total  noise.  If  any  boundary  has  a  value  of  NCELL  =  0,  that 
boundary  element  Is  bypassed  in  the  shock  noise  calculation.  It  is 
recommended  that  NCELL  =  8  be  used  for  each  element  unless  the  actual 
number  is  known. 

A. 2  EXAMPLE  CASE  INPUT  SELECTION 


To  illustrate  how  geometric  input  parameters  are  selected  for  a  complex 
nozzle  geometry,  an  example  is  presented,  taken  from  Reference  (1).  The 
example  nozzle  exit  geometry  is  shown  in  Figure  A-3.  Consideration  of  this 
figure  indicates  information  over  a  45°  sector  of  the  flow  field  will  be 
sufficient  to  describe  the  complete  flow  field.  This  is  one-eighth  of  a 
circle,  thus  LPHI  =  8.  Neither  axial  total  similarity  or  dissimilarity 
exists  so  ISYM  is  0.  Counting  the  number  of  closed  contours  indicates  a 
value  of  NEST  of  12,  where  one  is  included  for  the  ambient  or  external 
field.  Values  of  PT  and  TT  must  be  provided  for  the  exit  state  existing  just 
within  each  of  these  contours.  Values  of  XE,  ALPO,  LEAV,  NUM,  KN,  DEQ,  DS, 
and  NCELL  must  be  provided  for  all  the  contours  except  the  first  which  is  the 
boundary  at  infinity.  Values  of  these  parameters  for  the  contours  shown  in 
Figure  A-3  are  now  considered  in  the  following  discussion. 

Boundary  2:  Description  of  this  boundary  starting  at  45°  to  the  system 
baseline  is  convenient.  Thus  ALPO  =  rt/4  radians.  Since  each  90°  sector  of 
the  contour  is  identical  with  the  proceeding  one,  LEAV  =  4.  Since  the  pro¬ 
gram  assumes  straight  lines  to  exist  between  successive  boundary  points, 
description  of  this  boundary  is  possible  with  only  three  points  for  each 
quadrant.  These  are  P(l,2),  P(2,2),  and  P(3,2).  Each  point  is  described  by 
its  distance  from  the  system  origin  and  the  angle  between  the  line  joining  it 
with  the  origin  and  the  line  joining  the  proceeding  point  with  the  origin. 
Note  that  no  value  of  RA  is  given  for  the  point  P(0,2)  since  it  will  be 
identical  to  RA(3,2).  The  value  of  NUM  for  boundary  2  will  therefore  be  3. 

Boundary  3:  This  contour  has  eight  symmetric  leaves,  thus  LEAV  =  8. 

ALPO  of  0.0  is  as  convenient  as  any  other  value.  The  eight  points  indicated 
P(l,3)  through  P(8,3)  probably  are  sufficient  to  describe  the  boundary.  Thus 
NUM  =  8. 
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Boundary  4:  Since  this  is  a  circle  about  the  origin  it  can  be  divided 
into  a  convenient  number  of  leaves  and  only  one  point  given  for  each  (NUM  = 

1) .  If  a  hundred  boundary  points  are  desired;  set  LEAV  =  100,  DALP(1,4)=tt/50 
and  RA(l,4)  equal  to  the  circle  radius. 

Boundary  5  through  12:  Each  of  these  contours  must  be  described  indi¬ 
vidually  unless  certain  artifical  changes  are  made  in  the  arrangement.  A 
partial  representation  of  Boundary  9  is  shown  in  Figure  A-3.  Note  that  suc¬ 
cessive  points  on  the  boundary  are  obtained  by  progressing  around  the  bound¬ 
ary  in  a  counter-clockwise  fashion.  In  order  to  reduce  the  labor  of  repre¬ 
senting  each  circle  separately,  a  straight  line  can  be  drawn  connecting  each 
circle.  Two  contours  can  then  be  visualized,  one  consisting  of  the  outer 
halves  of  the  circles  and  the  lines,  the  other  consisting  of  the  inner  halves 
of  the  circles  and  the  lines.  Each  contour  has  eight  leaves  and  only  one 
need  be  represented  by  the  programmer.  Since  this  technique  requires  the 
computer  to  integrate  along  each  straight  line  twice  in  the  course  of  computa¬ 
tion,  it  will  definitely  increase  the  computational  time  over  the  method 
wherein  each  contour  is  represented  separately. 
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APPENDIX  B 


» 


OUTPUT  DESCRIPTION 


The  output  format  is  generally  self-explanatory.  The  input  data  is 
first  printed  out,  using  the  same  nomenclature  previously  defined  in  Table  2. 
Nozzle  exit  plane  flow  conditions  (static  temperature,  velocity,  Mach  number, 
momentum  flux,  and  enthalpy  flux)  are  then  printed  out  for  each  boundary 
contour. 

At  each  axial  location  specified,  the  radial  and  tangential  distribu¬ 
tions  of  flow  field  properties  are  printed  out.  After  the  flow  field 
information,  the  noise  characteristics  of  that  particular  axial  station  are 
then  listed. 

Following  all  of  the  axial  station  flow  field  data,  a  summary  table  of 
the  noise  characteristics  (SPL  spectra,  PNL,  PWL,  OASPL)  is  given. 

Appendix  C  contains  an  input  deck  card  listing  and  output  printout  for  a 
sample  case  run.  This  particular  case  is  for  a  7-tube  nozzle  presented  in 
Reference  2.  Only  a  portion  of  the  total  output  is  shown  for  brevity,  but 
the  formats  of  the  various  output  data  are  all  included. 

Two  warning  flags  are  built  into  the  program.  The  first  is  a  case 
termination  flag,  which  occurs  whenever  an  input  total  pressure  (PT)  is  less 
than  the  input  static  pressure  (PS).  The  flag  message  is  as  follows: 

****ERROR  -  MACH  NO.  SQUARE  IS  NOT  GREATER  THAN  ZERO  -  CASE  WILL 
TERMINATE**** 


The  second  flag  is  a  warning  detected  in  subroutine  SLICE,  which  occurs 
whenever  the  number  of  turning  points  (NTP)  is  found  to  be  greater  than  2. 
The  flag  message  is  as  follows: 


WARNING  -  NO.  OF  TURNING  POINTS  IS  GREATER  THAN  2  AT 
KA  =  _ ,  X  =  _ ,  ITH  =  _ ,  THETA  =  _ ,  NTP  = 


where  KA  is  the  axial  station  number,  X  is  the  axial  location,  ITH  is  the 
observer  angle  index,  THETA  is  the  observer  angle  in  degrees  (6p),  and  NTP  is 
the  number  of  turning  points  found.  The  two  outermost  turning  points  are 
used  and  those  inboard  of  these  two  are  discarded  in  such  cases,  since  the 
acoustic  shielding  model  can  only  accommodate  up  to  2  turning  points.  The 
noise  output  at  those  values  of  9j  where  this  warning  appears  should  be 


treated  as  suspect,  since  the  acoustic  shielding  effects  are  not  properly 
modeled.  This  is  most  likely  to  occur  in  the  initial  mixing  regions  of 
multitube  nozzles  where  multiple  peaks  in  the  aximuthally-averaged  velocity 
profiles  are  likely  to  occur. 
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APPENDIX  C 


SAMPLE  OUTPUT  LISTING 


An  example  case  of  a  7-tube  multielement  nozzle  is  described  herein, 
selected  from  one  of  the  data/theory  comparison  cases  presented  in  Reference 
(3).  The  nozzle  consists  of  a  hexagonal  array  of  0.875-inch  diameter  tubes, 
with  a  spacing/diameter  ratio  of  3.  The  acoustic  arena  is  a  9  ft.  radius 
arc.  The  geometry  is  illustrated  in  the  sketch  below. 


Observer 


Centerline 


d  =  0.875  in. 
D  -•  3.0  in. 


7  Boundary 
Number 


The  input  data  cards  for  this  case  are  listed  in  Table  C-l.  Note  that  all 
geometry  input  lengths  are  in  feet,  and  all  input  geometry  angles  are  in 
radians.  The  output  listing  for  this  case  follows  Table  C-l. 


Table  C-l.  Input  Data  Card  Listing  Sample  Case 


5R52  9  LI  1 D-C  6-7?  16. 4  71  INPUT  DATA  CARD  LISTING  --  M*G 

cwo  7-TURF  AS a?, I  Nr72Lt  -  9J*??00  F PS  -  Uj»160n  CtG 

$ INPU  T 

NEST=8,  LPHI=12,  1  SY  =0,  13011=100, 

RU2M=3,  DTHM=0.1,  PS  =  2  1  16/ 

ATOTAL=0. 029231,  D £3=8*0.0729167,  DS=8*0. 0729167,  N  CE  L  L  =  3*  3  , 

K  N  =  3  *  1  /  X  £  =  8  *  0  , 

GAM=1.35,  CP=6619, 

ALPJ  =0.  0/S. 96144, 0. 7  2  5  44  7  ,  1.  772 64  .  2. 8198/3. 86  70 /A. 9142, 0.  J, 

LEA  V  =  0/  6*1  ,24/  ,9111=1  ,6*24  ,  1,  KN  =  8  *  1  ,  X£=80, 

D  ALP (  1 , 2  )  = 

.033396,.'  45590, .054084, .059862, .063450/. 0651 68/ 

. 0631 63 /.C63A50/. 059 862,. 03408A, .045590/. 033596/ 

.U17C 39, -.005317,-. 334356/-. 068 98 4, -.103591, -.126562/ 

-.126 56?/-. 103591/-. 068984,-. C34336/-. 00531 7, .017039, 

D  alp ( 1 , 31  = 

.033596/. 0455 90, .054 084,. 059862,. 063 450, .0651 68/ 

.06 5 163/. 06 34 50/. 053862/. 05 4084, .045 5 90 ,.033596/ 

.0170 33/-. 00 5317,-. 334336/-. 068 984, -.1035 91 ,-.126562/ 
-.126562/-. 103591, -.068984,-. 034336,-. 005317, .017039, 

D  ALD ( 1 / 4 )  = 

.  03  3  5  96 /.[,  4  5  590/.  0  54  084,.  05  9862  , .063450, .  06  51  68/ 

.0651  63  ,.  363450,.  359862/. 054084  ,. 045590,. 03  35  96/ 

.01  7039, -.03531  7, -.034336/-. 063 93 4,-. 1035  91  ,-.  126562/ 

-. 126562,-. 103591, -.063984,-. 034336/-. 005317,. 01 70 39, 

DALPC 1/5 ) = 

.  0355 96 /. C  45590,.  0  54 084, .059862/.  063450/.  06  5168/ 

.065  1  63/.  06  345  0,.  0  59 862,.  05  4084  , .  045  590/.  03  35  96, 

.017039, -.005317, -.034356,-. 068 93 4, -.103591,-. 126562/ 

-.  1  26562,-.  103591  ,-.063984  /  -  .034536/-. 00531  7,.  01  7039, 

0  ALP ( 1 / 6 1  = 

.  0  3  3596/.  045  5  90/.  0 54 0 8 4 , . 0 5 9 862  , . 0 6 3 4 50 / . 06  51  68/ 

. 06 51 63 ,.C6 3450,. 059862,. 054084, .045 5 90/. 03 35 96/ 

.017039, -.005317,-. 034336,-.  06 8  93 4, -.103591, -.126562, 
-.126552/-. 103591, -.068984,-. 034536, -.005317,. 017339, 

0  ALP  < 1 / 7  1  = 

.0335 96,. 045590,. 054 084, .059862, .063450,. 06 51 68/ 

. 065 1 63 /. 06 34 50/. 0 59 862/. 05 4084, .045 5 90 ,.033596/ 

.01 7039, -.0053 17, -.034336,-. 068984, -.1035 91, -.12656 2/ 
-.126562/-. 10 3591, -.068 98 4, -.034336/-. 005317, .01 70 39, 

R  A  (1  ,  2)  = 

. 12392, . 13145, .13759,. 14212, .14439,. 14583, 

.14489, .14212, .157  59, .13145/.  12  39 2, .11529, 

.1059 6/. 09 6 46/. 03748, .07991/. 07476, .07292/ 

. 0747 6/. 07991, .08743,. 09 646/. 10 59 6/. 11529, 

R  A  <  1  , 3>  = 

.12392,.13145/.13759,.14212/.14489,.14583, 

.1448 9,.  1  4  2  1  ?, .157 59,.  13145/.  12  39 2, .11529, 

.  1059 6/. 09646/. 03743/. 07991  /. 07476/. 07292, 

. 07476/. 07991, .03748/. 09646/. 10596, .11529, 

RA< 1 ,4) = 

. 12392, .13145/.13759,. 14212, .14489, .14583, 

.  1448  9,. 1  421  2, .137 59 ,.13145/. 12  39 2, .1  1529, 

.  1  0596/.  09646 ,.03  748,.  07991  ,. 07 476/. 0  7292  , 
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Table  C-l.  Input  Data  Card  Listing  Sample  Case  (Concluded). 


5  R  3  2  9  01  10-06-77  16.471  *«*  INPUT  DATA  CARO  LISTING  --  M*G‘d 

.07476,. 07991, .087 48, .09 646, .10 59 6, .115 29, 

B  A  <  1  ,  5)  = 

.12392,.13145,.13?59,.14212,.14489,.14583, 

.14489, .14212, .15759, .13145, .12 39 2, .IIS  29, 

.1059  6.  .09646, .03748, .07991, .07476, .07292, 

.07476, .07991, .08743, .09646, .10 59 6, .11529, 

RAC  1  ,6)  = 

.12392, .13145, .137  59, .14212, .14489, .14583, 

.14489, .14212, .13759, .13145, .12  39 2,  .11529, 

. 1059 6,.  0 9646, .06748,. 07991, .074 7 6, .07292  , 

.07476, .37991, .03743, .  09 646, .10596, .11529, 

R  A  (  1  ,  7)  = 

. 12392,. 13145, .13759,. 14212, .14489, .14583, 

.1448 9.. 14212, .137  59  , .13145, .12392, .11529, 

.  10596,  .09646, .08748, .07991  , .07476,. 07292, 

.07476. .37991, .03748, .09 646, .10 59 6, .11529, 

DALP(1»8)=3.26lg,  RAC1, 31=3. 036453, 

AlPHMCsO.  5,  BE  T  A'A  c  -=3  .  2  5  , 

FMINslOO,  F*IAX  =80000,  NUMAN  G  =  1 ,  DIST  =  9.0, 

KX  =  34,  X  =  Cl  .  0  729167,  STRFR  =  0.01, 

0S;G  =  13*j.  33729167,  14*0,  NOC/=10*20,14*3, 

BETA  IN* 1 5*0.  C, 

N  P  R  1  N  T  = 

N  C  A  S  E  =  1  , 


PT=21 16. 7*575?,  TT*SU0,7*160s 
S 


M  f. 


4  Aft  I 


COMPUTATION  Of  ACwO-ACOi)'»T  ll  P^'PEHIIES  Of-  SuJM*t SiO  NPZ2LLS  _ 

CASE  V)#  1  C^n  7- r I AP  =  2.  1  G'7?LF  -  VJ=??„0  F^S  -  TTJ=1*00  OfG-P 


KXs  ?4 

CH=  l.l^C 


NEST  =  6 

GAMs  1.3*^ 


If, PUT  DMA 


LPHls  12 

CP=  6619.0 


ISYMs  J 

r  t  hm=  .uoo 


NPF>  I  NT  =  6 

Pu?M  =  3.0000 


CM=  .075 
PS=  2116.0 


COMPUTATION 

-  -  slice  no. 

MESH  COK,T PQL 

X  -  - 

. ../  TUP9ULENCE 

CONSTANTS  ...  . 

-DELTA - 

- OSIG 

.  PMIN 

GOV  RfTA 

-MU  -- 

i 

.C7292 

.00729 

c.coooo 

20  0.00 

4.00 

.20 

_  2 

.?41*7 

..  .0:729  . 

0.00000 

2o..  a. oo 

.  ‘•.00  -  .._ 

.20  _ 

3 

.  1  l  575 

.0*729 

0.00000 

?■)  0.00 

4.00 

.20 

_  _  .  4  . 

.14583 

. 0  C  7  2^ 

0.0000c 

2)  .—0.00 

4,00 

.20. 

5 

.  l«374 

.or  729 

0.000™ 

20  0.00 

4.00 

.20 

6 

.23150 

iC( 729 

O.ooron 

20  0.0  0 

‘.00 

.20 

7 

.29167 

.Of  72v 

0.00000 

20  0.00 

4.00 

.20 

.36748 

.yi?v 

C.OOCOv 

?"  .  o.oo 

‘.00 

.20 

9 

.46299 

.0 '7?u 

0. ‘0000 

2  r.  o.co 

4.00 

.20 

1  3 

.58393 

.0 >729 

O.GOOOC 

20  -0.00 

.  ‘.00. 

.20 

1  1 

.73^95 

.<k  nc. 

0.0000c 

0  0.00 

4.00 

.20 

.  __  .  _  12 

_ .52598 

•  O'."4'')- 

r  .  ^  •>  0  0  f- 

0  0.00 

4  .  OA 

.20 

f  3 

)  .16667 

.*>11-7 

0  .  ;  0  0  C  ) 

'  o.co 

4  .on 

.20 

14 

1  .46941 

. 1  *♦  7 

0.00000 

0  ^0.00 

4.00 

.20 

t  5 

l  .^147 

c.v 

o  C .  C  n 

4.00 

.20 

16 

?. 99?  39 

.02373 

0 .  rot  0: 

0  .  .  4.00 

..  4.00 

.20 

17 

2.9990? 

.0244^ 

0  .  •*  0  c  0  c 

)  4.00 

4.00 

.20 

La 

3.7 J 193 

_.Q37 0  4 

0.00003 

0  4.00 

4.00 

_  _»20 

19 

4.66667 

.04667 

0.00000 

0  4.00 

4.00 

.20 

.  _  .20-.  ... 

5.87963  . 

...05680  . 

0.000QCL 

0  4.00 

— 4  •  QO. . 

...20 

21 

7.40787 

.07408 

0.00000 

0  4.00 

4.00 

.20 

-_-  ._  _22_ 

.  9.33333 

*09333 

_ 0 . Q  0  0  00  _ 

Q  _  4.00 

4.00 

.20. 

23 

11.75926 

.11759 

0.00000 

0  4.00 

4.00 

.20 

_ 24 

14*81SX4 

.14816  _ 

sSK'lBI  •  W  i  Ii  1 1 1 

0  _  _  4.00 

4.00 

...  .20 

-.  Jt£L_2)  =  _  ..  ALPOI  21=  5. 

9614  .LEAV.L.2l  = 

L  .  NUM(  ?)T  24 

KN(  2) =  1 

_ _ _ -OALPt  i. 

21-  .3336 

.  RAL.1.  2)5 _ . 

1239—  DALP.l  2,  -2) 

=  .... (LA 5% _ R6  (  .2,  2_>_5  -- .,  1  31 5 

PALP t  ^3«_21= _ .1541 _ PAt  3.  2)  = _ .137.6 _ QALP.  .1-4.1-21.= _ .Q599-.  pa  t  4i  2>  = _ .1*21 


-OALPl  5.  -2)  = 

— . 0635 

PJU-S.  2 1  =. 

_*14.49_  _DALP  1  _6* _21  =_ 

.0652 _ RA1  6-._2)  =  . 

•  1458.. _ _  - 

PALPI  7.  2) = 

.  .3652  _ 

RAL.7, 

_2)  = 

.1449 

DALP  ( _ 8*._2)_= 

-  . 1635- 

_ RAL-8.  .21  = 

.1421 _ 

OALP(  9 .  ?>  = 

.3599 

PA  (  9. 

21  = 

.1376 

OALPUO*  2)- 

.0591 

PA  1 10,  215  _ 

-^i-315 

nALPtU-^-ZJs 

-*3456  _ 

— RA (11« 

21  = 

-.1239 

DALP  (12*  21- 

.0336- 

_  RA 112.  2) s  - 

.1153 - 

DALPtl3«_2)  = 

_*ai7o  ... 

21  = 

•  1060 . 

DALP-C14  *,  21- 

-  .0.05,1 

R A ( 1 4  *  2)  = 

.0965 

DAI  °  1 15  «  2)  = 

RA ( 1 5 1 

_2l£ 

_ .1PJ5- 

PALP  (18*  2)  = 

-.0690 

RA ( 1 6 ♦  2) = 

.0799 

.  .  .  - _ DAL3  ’17.  2 )  = 

-.1036 

_.RA  IJ  7  • 

2)  = 

•  o?‘a 

DALPtia.  21- 

-.1266- 

__  RACiet-  23  = 

.0729  _ 

OALP ( 1 ° •  2)~ 

-.1266 

WA  (1*9. 

?)- 

.074P 

nf.LPiro.  2 > - 

-  .  1  C  36 

.PA  (20,  2)  5  . 

.0799  -  _  _ 

PALP  12  L«_2) 

.-.1690  _ 

PAl  21. 

21  =. 

.0875 

..  DALP(22,.  2)  = 

- .  Q2‘  3 

°A_122j_ 2J  * 

.0965 

DALP (23 ■  2) = 

-.0053 

PA ( 23* 

2)  = 

•  1  060 

PALP ( 24  *  21  = 

.0173 

RA (26,  ?) = 

.1153 _ 
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COMPUTATION  QF  AERQ-ACOUST IC  PROPERTIES  OF  SUPPPESSnR  NOZZLES 


CASE  NO, 

1 

CPD  7-TUBE  Ap»?,3 

N077LE  - 

VJ=??00  EPS  - 

T  T  J= 1 ftOO 

OEG-R 

xe< 

3)* 

0.00 

ALPO( 

3)  = 

.  7?54  LEAV(  31 

=  1 

NUM<  31= 

?4 

*N(  31= 

1 

— 

OALP<  ), 

31  = 

,03  3ft 

P  A  (  l. 

3)  = 

.  1?39 

dalp<  ?* 

31  = 

.0456 

PA  (  ?, 

31  * 

.1315 

DALP  <  3. 

3  >  = 

.OS41 

R A  (  3, 

3)  = 

,137ft 

D ALP (  4, 

3)  = 

.0590 

RA(  4, 

31  = 

.  14?1 

DALPi  5* 

3)  - 

.0635 

PA  <  5, 

31  = 

.  1449 

PALPI  6, 

31  = 

•  U66? 

PA(  ft, 

31  = 

.1458 

DALP (  7* 

3)  * 

.065? 

PA  (  7, 

3)  = 

.1449 

l)ALP<  8* 

31  = 

.0635 

PA  (  8, 

31  = 

.  14?1 

DALP (  9, 

3)  = 

•  0S‘>9 

RA  (  9, 

31  = 

.  \  37ft 

Dai  p n  o. 

31  = 

.054  1 

RA  ( 10, 

31  = 

.1315 

OALPUI* 

31  = 

•  0456 

RAUlf 

31* 

.1239 

DaL  p ( )?, 

31  = 

,03  16 

Rail?, 

31  = 

.nsj 

DALP  <13* 

3)  = 

.0170 

PA(13* 

31  = 

•  1060 

DALP ( 14, 

3 )  = 

-.(’05  1 

RA  (  14, 

31  = 

.0965 

DALP  < IS* 

3 )  =  - 

,0  34.3 

P  A  (  1  5  , 

31  = 

.0875 

DALP  (  1ft* 

3)  = 

-,Oft90 

RA (1ft, 

31  = 

.0794 

OALP ( 17, 

3)  *  - 

.1036 

RA  ( 17, 

31  = 

.0748 

DALP (18, 

31  = 

I?ftft 

RA ( |8, 

31  = 

•  0  7?9 

0ALP<1P* 

3)  =- 

.  1  ?66 

RA( 19, 

31  = 

.0748 

flALP(?n, 

31  = 

-,103ft 

RA  <?0» 

31  = 

.0799 

DALP121. 

; 

.0690 

RA  ( ?  1 , 

31  = 

.0875 

DAI P(??, 

3 )  = 

-.0343 

R  A ( ?? , 

31  = 

.0  965 

DALP (23, 

3  >  =  - 

.9053 

RA(?3, 

31  = 

.lOftO 

DALP(?4, 

31  = 

.0170 

PM24. 

31  = 

.1153 

x£< 

4)* 

0.00 

ALPOI 

4)  = 

.7706  LEAV1  4) 

=  1 

NUM (  41= 

?4 

l*N<  41  = 

1 

OALP(  1. 

41  = 

.0336 

8A<  It 

41  = 

.1239 

D AL  P (  ?, 

41  = 

,0456 

RA  (  P * 

41  = 

.1315 

DALP (  3 , 

4)  = 

.0501 

RA<  3, 

4)  = 

.1376 

PALPI  4, 

4)  = 

.0599 

R A  (  4, 

41  = 

.14?1 

OALP (  St 

4)x 

,0635 

PA<  5, 

41  = 

.1449 

OALP (  ft. 

41  = 

•  06  5? 

RAT  ft. 

41  = 

.1458 

DALP!  7. 

4)  = 

.065? 

RA(  7. 

4)  = 

.1449 

(ialpi  s, 

4)  = 

.0635 

RA(  ft. 

41  = 

.  14?1 

DALP  (  9, 

4)  = 

.0599 

RA<  9, 

41  = 

.  1  37ft 

DALP  (  1  0, 

41  = 

.0541 

RA  ( 10, 

41  = 

.  1315 
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COMPUTATION  of  afro-acoustic  propfrtifs  of  SUPPRFSSOP  nozzles 


CASE  NO. 


CPD  7-turf  A  W  =  2 . 3  NOZZLE  -  VJ=2?00  FPS  -  TTJ=l600  OFG-P 


AXIAL  LOCATION  =  .67292  (X/OEO  =  1.00000 


M 

R 

ANGLE 

IJ 

OF  NS ITT 

TFMP. 

U/URL  7 

TURB. I  NT 

.  P/DEO 

9 

.05833 

17. 1A 

0  .  o  0 

.0022835 

541. 00 

0.00000 

0.00000 

.  8  0  0  0  u 

9 

,f5833 

21  .A3 

0 . 0  0 

.0022835 

540.00 

0.00000 

0 .00000 

.80000 

9 

.'.58  3  3 

25.71 

0  .  '•  8 

.0022835 

54 .00 

0.00000 

0.00000 

.80000 

9 

.".58  33 

30. 00 

0.'  0 

.0022835 

540.00 

0.00000 

0.00000 

.80000 

1  0 

.  0  6  5  6  3 

0.00 

203. 37 

.00) 6986 

725.94 

.09246 

.0A216 

.90000 

10 

. G  658 3 

A.  29 

188.  ".0 

.0017420 

707.84 

.07638 

. 0  37  1  C 

.90000 

1  0 

.085(33 

8.57 

95. 3 A 

.5018676 

660.25 

.04335 

.02447 

.90000 

1C 

.06583 

12.66 

33. A9 

.(  020323 

606.73 

.01523 

.01067 

.90000 

1  0 

.  r-6S63 

17.1a 

S  .  8  A 

.'■0216  16 

669.91 

.00266 

.00325 

.90  000 

10 

.06563 

21  .A3 

0.0  0 

. ; 022835 

64  .00 

0,00000 

0.00000 

.900  0  0 

10 

.0  656  3 

26.71 

0 .  A 

.0022835 

54 " . o  C 

0.00000 

O.OooOC 

.90000 

1C 

.  0656  3 

3  3.00 

o. 

. 1C  228 35 

54i.CC 

0.00000 

0.000O0 

.90300 

1 1 

.0739? 

o.oe 

1  378.  A5 

. 1011598 

1 06  3.19 

.6267? 

.13779 

1 .00000 

1 1 

.0729? 

3.75 

128?. 33 

.0011754 

1 349.07 

.5830? 

.  13990 

1.00000 

11 

.0729? 

7.50 

1015.37 

.C01?  396 

994.75 

.46165 

. 1 3664 

1  .00000 

11 

.0729? 

11.25 

608.69 

.0013846 

896.59 

.27674 

.  10947 

1  .uoooo 

1 1 

.0729? 

IS. 00 

2 A?. 76 

.0016504 

747.13 

. 11037 

.06941 

1 .00000 

1 1 

.07292 

18.75 

57.25 

.  ,10  1  9685 

629.60 

.02600 

,01884 

1  .00000 

1 1 

.0729? 

22.50 

6  .  A  A 

.0021637 

569.90 

.0  029  3 

.00379 

1  .coooo 

11 

.0729? 

26.25 

0.00 

. u0??835 

54. .00 

0.00000 

0.00000 

1  .00000 

1 1 

. "729? 

30.00 

O.ro 

.0022835 

541.00 

0.00000 

0. 00000 

1  .  0O000 

1? 

.08021 

0.00 

21A5.7A 

.1010196 

1209.4? 

.97558 

.05248 

1  .  10000 

1? 

.08021 

3.75 

212A.5A 

.1010266 

1?M  .14 

.96594 

. 06461 

1.10000 

1? 

.08021 

7.50 

2025. 1a 

.0010525 

1  171.60 

.92075 

.10460 

1.10000 

1  2 

.08021 

11.25 

1717. A9 

.0011  OS? 

1  1  15.68 

.78087 

.1585? 

t  .  10000 

1? 

.08021 

1 5  •  C  0 

1059.87 

.0012270 

1004.97 

.48188 

. 16340 

1.1  0000 

1? 

.08021 

18.75 

363. 5A 

.0115447 

798.27 

. 16074 

.08533 

1 .10000 

1? 

.08021 

22.50 

53.70 

.0019665 

627.  f  3 

.02441 

.01880 

1  .  10000 

1? 

.08021 

26.25 

2.21 

.0021721 

567.69 

.00100 

.00334 

1.10000 

1? 

.08021 

30.00 

0.05 

.00226  .5 

641.00 

O.OOuOO 

0.00000 

1  .  1  0000 

13 

.08750 

0.00 

2199.19 

. " 0  09950 

1239.30 

.99988 

.0021  3 

1 .20000 

13 

.08750 

1.33 

2198. 94 

. 000996? 

1239.&7 

.99977 

.00350 

1  .20000 

13 

.  0 8750 

8.67 

2196.64 

.0009967 

1237.12 

.9987? 

.01148 

1 .20000 

13 

.08750 

li.OO 

2175.80 

.0010076 

1223.80 

.98925 

.04213 

1  .20000 

13 

.08750 

1  3.33 

2032.46 

.00 1 0508 

1171-47 

.92408 

.  1  189(1 

1 .20000 

13 

.08750 

18.67 

1500.16 

.0011379 

108  3.6\ 

.68206 

.18610 

1 .20000 

13 

.08750 

2C.00 

616.35 

.0013828 

891 . ?A 

.28023 

.  1  3  0  29 

1  .20000 

13 

.08750 

23.33 

110.21 

.0018367 

671.34 

.0501  1 

.03566 

1.20000 

13 

. 08750 

26.67 

7.85 

.n0?156? 

572.14 

.00367 

.00463 

I .20000 

13 

.  08750 

30.00 

0.0  0 

.">.j??S35 

540,00 

0.00000 

0.00000 

1.20000 

1  A 

.  "94  79 

8.00 

2 1 99.45 

.0009948 

1239.57 

1 .00000 

.00140 

1 .30000 

1  A 

. 0  9 A  79 

3.33 

2199. 45 

.,1009948 

1239.57 

1.00000 

.00161 

1 .30000 

lA 

. 89579 

8.67 

2199. AA 

.(  009948 

1239.55 

.  99999 

.00067 

1  .  3  0  0  0  0 

1 A 

•  09A79 

10.00 

2198.75 

.0009963 

1238,90 

.99968 

.00583 

1  .30000 

14 

.  09479 

13.33 

2)75.18 

.8010078 

1 223.48 

.98897 

. 04684 

1  .  30000 

l  A 

. 09A79 

18.87 

1911 .43 

.001(1745 

11  A7.62 

.86905 

.15684 

1  .  30000 

1A 

•  09A  79 

20.00 

1022.06 

.0012373 

996.57 

.46469 

.17621 

1 .30000 

1  A 

. n 9479 

23.33 

2  0  H  .  7  n 

.0016920 

728.79 

. 09489 

. 06070 

l . Tinnn 
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COMPUTATION  OF  AFKO- ACOU"'  T  I  C  PPOPFpTIFS  OF  SUPPRF SSOP  U077L  F  S 

C ASF  NO,  1  CPD  7-TllPf  AR=7.  1  N077LF  -  VJ=??00  FPS  -  TTJ=1600  DFC-P 


AXIAL  LOCATION  =  ,v7?9?  (X/DFQ  =  l.OOOUO) 


M 

R 

anglf 

U 

OF  NS I T  Y 

TF  MP. 

U/UWEF 

TUP8. INT 

.  p/ofo 

14 

.  .',04  70 

76.6  7 

14.75 

.'071 165 

587.59 

.00671 

.00707 

1  .  30000 

14 

. 79479 

30.03 

0.  ?0 

.0077835 

546,00 

0.00000 

C. 00000 

1  .30000 

IS 

. 1070B 

0 . 0  3 

7 1 69 ,4S 

.  .1009946 

1739.57 

1  .COC'OO 

.00031 

1  .40000 

15 

.  1  0708 

3.00 

7199.45 

.0009948 

1739.57 

1.00000 

.00 1 88 

1  .40000 

IS 

.  1  0  7  C  8 

6.00 

8)99.45 

,f  009948 

1739.57 

1  .  0  0  0  0  0 

.00140 

1  .40000 

IS 

.  10708 

9.00 

7199.45 

. O0O994B 

1739.56 

1 .00000 

.00159 

1  .40000 

IS 

.  1  0?0« 

17.00 

7198.59 

.000Q954 

1  738.7? 

.99959 

.00773 

1  .40000 

IS 

.  1  0708 

15.00 

7161 .67 

.  y  C  1  0  1  35 

1716,6? 

.96280 

.06  160 

1  .40000 

IS 

.  1C?0° 

1  B .  o  0 

lBol .51 

. 0010971 

1  179.1? 

.61907 

. 1 7665 

1  .40000 

15 

,  1070« 

71.00 

6  3  4  «  46 

.  i*OI  ?9 1  7 

954.6) 

.37939 

.  16055 

1  .40000 

IS 

.  1  0708 

74 .00 

146.71 

.0017736 

695.16 

.  0667  0 

.  08563 

1  .40000 

IS 

.  10708 

77.00 

9.45 

.0,071459 

574.63 

.00630 

.00516 

1  .40000 

IS 

. 1 0708 

7  3.00 

C.C  0 

. 0078835 

S40.CA 

0. 000 00 

0. 000(  0 

1.40000 

16 

. 1 093« 

0.00 

7199.45 

.0009948 

1239.57 

1.00000 

.00104 

1  .50000 

16 

.10938 

8.0  0 

7199.45 

. 0009946 

1739.57 

1  .OOOCO 

.0027? 

1.50000 

16 

.10938 

6.0  0 

7199. 4S 

.4009948 

1739.57 

1.00000 

.0077? 

1  .50000 

16 

.10938 

p.or 

7199.45 

.0009948 

1739.56 

1 . 00 J00 

.00)8? 

1.50000 

16 

.10938 

17.  CO 

7196.81 

. 0009953 

1736,05 

.99971 

.00679 

1,50000 

16 

. 1 093« 

1 S .  o  A 

7161.18 

.3010137 

1216.41 

.98260 

.06199 

1 .50000 

16 

.10938 

1  H  .  0  C 

1 7S2  .  C  p 

.0010996 

1  121 .44 

. 79660 

.18152 

1 .50000 

16 

. 1 0938 

71,0  0 

710. 4S 

.0013394 

9?C .63 

.32301 

.14531 

1  .50000 

16 

.  1  0938 

74.Q3 

98.76 

.0016815 

662.41 

.04468 

.03259 

1  .50000 

16 

.10938 

77. CO 

4.71 

.0071636 

564.70 

.0019? 

.00363 

1 .50000 

16 

. 1 0938 

30.00 

0 . 0  0 

.0072835 

SBC. CO 

O.OCOUO 

0.00000 

1.50000 

17 

.  1  1667 

0.00 

7199.45 

.0009946 

1239.57 

1 .00000 

.00706 

1  .60000 

17 

. 11667 

7.73 

7199. 4S 

.0009946 

1739.57 

1.00000 

.00734 

1 .60000 

17 

.11667 

S.4S 

7199. 45 

.0009946 

173Q.57 

1  .ooooc 

.00178 

1 .60000 

17 

.  1  1667 

«.  18 

7199.45 

. 0009948 

1739.56 

1 .00000 

.00)57 

1 .60000 

17 

.  1  1667 

10.91 

219R.91 

.0009957 

1239.04 

.999  75 

.00557 

1  .60000 

17 

.1 1667 

13.64 

7174.93 

.0010080 

1773.34 

.98885 

.04751 

1  .60000 

17 

.11687 

16.76 

1906.45 

.0010753 

1146.73 

•  866  7  4 

•1567c 

1  .60000 

17 

.  1  1687 

IP.  09 

1 079. A3 

.0017337 

999.94 

.47241 

. 1 7474 

1.60000 

17 

. 1 1 687 

71 .87 

735.07 

.  no  1  6,608 

747,44 

.  1  068  7 

•  0  65  3 1 

1.60000 

17 

. 1 1 887 

74. SS 

?3.71 

.0070856 

591.7? 

.00942 

.00875 

1  .60000 

17 

. 1 1 887 

77.  27 

O.nn 

.0077835 

540.00 

C. 00000 

C  .  0  0  0  0  0 

1 .6000C 

17 

.  1  1867 

3C.0  0 

0  .  L  0 

.0077635 

543.00 

0.00000 

0.00000 

1  .60000 

1« 

. 17306 

C  .  OC 

7199.45 

•0009946 

1739.57 

1 . D0O00 

.00  154 

1  .  70000 

18 

. 1 7396 

7.73 

7 1 95 ,45 

. 0009948 

1239.57 

1 .00000 

.00141 

1  .70000 

IB 

. 17396 

5.  AS 

7199.45 

.0009948 

1739.56 

1 .00000 

.00162 

1 .70080 

IB 

.17396 

8.18 

7199. 1 3 

.3009950 

1739.75 

.99965 

•  0  04  36 

1 .70000 

IB 

.17396 

10.91 

7186.71 

.  0  o  I  n  o  l  1 

1731 .46 

.99512 

.  076  1  7 

1  .  70000 

IB 

.17396 

13.64 

7051.20 

.001 0465 

1176.33 

.93259 

.114  0  3 

1.70000 

1  B 

.17398 

16.06 

1433. 87 

.0011478 

1374.76 

.65192 

. 1 6468 

1 .70000 

18 

.  1  7396 

19.09 

478. 50 

.0014574 

646. : 9 

.21756 

.10756 

1.70000 

IB 

.17398 

71  .82 

62.84 

.0019401 

635.56 

.07857 

.07146 

1  .70000 

18 

.17398 

74. 5S 

7.80 

.0071965 

561.37 

.00127 

.00342 

1 . 70000 

18 

.17398 

77.27 

O.nn 

. 00778  35 

54O.0C 

0.00000 

0.00000 

1  .70000 

1  8 

. 17396 

7n.no 

A  .  r, 

. 00776 35 

543.cn 

o .  r.  n  n  o  o 

n .  nonnr, 

1  .  7 n (> nil 

71 
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COMPUTATION  OF  AERO-ACOIJSTIC  PROPERTIES  Of  SUPPRESSOR  NOZZLE'S 
CASE  NO,  1  CPO  7-TURf  AP=2.3  NOZ7LE  -  VJ=2?00  FPS  -  TTJ=1600  DEG-R 


AXIAL  LOCATION  =  .1729?  (X/DEQ  =  1.000*01 


M 

R 

ANGLE 

U 

PE  NS  ITT 

TEMP. 

U/UREF 

TURB. 1NT 

.  R/DEU 

19 

.13125 

3.00 

2199.19 

.•1009950 

1239.3  ;■ 

.99988 

.00325 

1.80000 

19 

.13125 

2 . 5  A 

2198.99 

.1009951 

1239.12 

.99979 

.004)? 

1.60000 

19 

.  1  3  IPS 

5. OP 

2197.42 

. 3009962 

1237.75 

.99908 

.00926 

1 . 80  0 0 0 

19 

. 1 3125 

7.53 

2185.40 

.3010030 

1229.37 

.99361 

.02954 

1  .80000 

19 

.13125 

i  o.or 

2106.46 

.0010321 

1194,73 

.95772 

.08446 

1  .80000 

19 

. 13125 

1  2.50 

1776.53 

. bC 1 0963 

1  124.8" 

.60771 

.15885 

1  .BCOOO 

19 

.13125 

15.00 

1  038.2C 

.0012327 

li.Oo.33 

.47203 

.1618? 

1  .80000 

19 

.13125 

17.50 

318.73 

.3015743 

783.25 

.14491 

.07710 

1.80000 

19 

•13125 

2-j.or 

47.64 

.3019859 

62C. 91 

.02166 

.f.  1638 

1.80000 

19 

, 1 3 1 ?S 

22.5C 

2.88 

.0021988 

560.79 

.00131 

.00320 

1 .60000 

P 

.13125 

25.00 

0.00 

.0022835 

540.01 

0.00000 

0.00000 

1.80000 

19 

.13125 

27.  SC- 

O.CC 

.0022835 

546.00 

o.mcoo 

0.00000 

1.80000 

19 

.13125 

30.  90 

0.  ,1 

.0022835 

54c. 00 

0.00000 

V.  00000 

1.60000 

20 

. 1 3854 

*.oo 

2145.68 

.0010196 

1209.39 

.97555 

.05255 

1.90000 

20 

.13854 

2.  50 

2129.85 

.1010249 

12C3.09 

.96835 

.06116 

1  .90000 

20 

.13854 

S.OC 

2061.96 

.0010438 

1181.30 

.93749 

.06886 

1  .90000 

21 

.13854 

7.50 

1867.77 

.0010820 

1139.60 

.84920 

.13235 

1 .90000 

20 

.13854 

10.00 

1476.12 

.0011483 

1 0  7  3 , 86 

.65295 

.16129 

1  .90000 

20 

.13854 

12.50 

756.51 

.0013212 

933.33 

.34395 

.129)3 

1  .90000 

20 

. 1 7854 

15.00 

231.75 

.0016626 

741.66 

. 1 C537 

.05730 

1  .90000 

20 

. 1 3854 

17.53 

97.77 

.0020185 

61 i.90 

.  C  1  7  1  7 

.01295 

1.90000 

?0 

.13854 

20. OC 

2. SC 

.0022007 

560.32 

.00114 

.00315 

1  .90000 

20 

. 1 3854 

22.50 

0.00 

.0022835 

540.00 

0.00000 

0.00000 

1 .90000 

20 

.13854 

25.00 

0 . 0  o 

.0022835 

54  0.00 

0.00000 

0.00000 

1  .90000 

?0 

.13854 

27.50 

o.oo 

.0022835 

540.00 

0.00000 

0.00000 

1  .90000 

20 

. 1 3854 

30.00 

O.C" 

.0022835 

540.00 

0. 00000 

0.00000 

1.90000 

21 

.14583 

0.00 

1 778.1  | 

.00 1 1599 

1063.10 

.62657 

.13781 

2.00000 

21 

.14583 

2.31 

1303.87 

.3011715 

1 052.5? 

.59281 

. 1 3796 

2.00000 

21 

.14587 

4. (S? 

1100.68 

.0012170 

1013.17 

.50043 

.13495 

2.00000 

PI 

.14583 

6.92 

773.97 

.1013145 

938.05 

.35189 

.1 1671 

2.00000 

p  1 

.14583 

9.23 

414.00 

.0015002 

821 .91 

.18823 

.08055 

2.00000 

?1 

.14583 

1  l  .54 

154.16 

.0017639 

699.08 

.07009 

.03861 

2.00000 

21 

.14583 

1  3.85 

37.12 

.0020219 

610.3? 

.01688 

.01201 

2.00000 

21 

. 14587 

16.15 

4.97 

.0021826 

564.95 

. 00226 

.00328 

2.00000 

21 

.14587 

18.46 

0.00 

.3022835 

54a.  00 

0.00000 

0.00000 

P.00000 

?1 

.14583 

20.77 

0.00 

.0022835 

54J.00 

0.00000 

0.00000 

2.00000 

PI 

.  14583 

27.08 

0.00 

.0022835 

541.00 

0.00000 

0.00000 

2.00000 

PI 

.14583 

25.38 

0.00 

.0022635 

540.00 

0.00000 

0.00000 

2.00000 

PI 

.  14587 

27.69 

0. 3U 

.0022835 

540.00 

0.00000 

0.00000 

2.00000 

PI 

.14583 

30.00 

0.00 

.0022835 

540,00 

0.00000 

0.00000 

2.00000 

?? 

.15313 

0.00 

203.27 

.0016988 

725.85 

.0924? 

.04220 

2.  10000 

?? 

.15313 

2.14 

181.15 

.0017236 

715.40 

.08236 

.03875 

2. 1 0000 

p? 

.15313 

4.29 

131.99 

.0018010 

684.65 

.06001 

.03049 

2.1 0000 

?? 

.15313 

6.43 

73.29 

.0019140 

644,25 

.03332 

.01912 

2.10000 

pp 

.15313 

8.57 

31.17 

.0020427 

603.66 

.01417 

.00955 

2.10000 

pp 

.15313 

IP. 71 

9.03 

.0021524 

572,88 

.00410 

.00396 

2.10000 

pp 

.15313 

12.86 

1.19 

.  0  0? 1 963 

561,44 

.00050 

.00274 

2.10000 

pp 

.15317 

15.00 

O.AO 

.0022835 

540.00 

0.00000 

0.00000 

2.1 0000 

72 
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COMPUTATION 

OF  AERO- 

ACOUSTIC 

PP0PER1 IFS 

OF  SUPPRESSOR  N077LES 

CASF 

NO.  1 

c»n 

7-TURF  AR=? 

.3  N077LE 

-  V  J  =  , 

2?00  FPS  - 

TTJ=160i 

AXIAL 

LOCATION 

=  •  0  729?  (X/DEO  =  1  .000061 

M 

P 

ANOl  F. 

IJ 

HENS  I TY 

TEMP. 

U/UPEF 

TURR. I  NT 

.  P/DE  0 

2? 

.15313 

17.14 

0.40 

.0022835 

540. CO 

0.00000 

0 . 0  0  0  0  u 

2.1 0000 

2? 

.15313 

1  9.29 

0  .  1;  0 

.  00??835 

54  Ci  ,06 

O.OCCOO 

O.GOCOC 

2 . 1  000  0 

2? 

.15313 

?l  .43 

0.00 

.30?  ?H  35 

54  4 ,00 

0.00000 

0.00006 

2.  10000 

2? 

.15313 

23.57 

0  .  C  C: 

.C02283S 

54  C.O1' 

0  .  OOOoO 

0.00000 

2.  10000 

2? 

.15313 

?5 . 7  1 

O.1-? 

.C92283S 

54 0.0 V 

0.00000 

0.00000 

2.10000 

2? 

.15313 

27.86 

0. 30 

.0022835 

540.00 

0.00000 

0.00000 

2.10000 

2? 

.15313 

30.0  0 

O.'.O 

. 3022835 

54 ',.05 

0. 000O0 

0.00000 

2.10000 

23 

.1604? 

3.00 

4.61 

.0021890 

563,29 

.  0C2 1 0 

.00280 

2.20000 

23 

.1604? 

2.14 

3.77 

. 40  ? 1 444 

561 .03 

.0017? 

.00285 

2.20000 

23 

.  1604? 

4  .  ?9 

2.16 

.  60220  76 

858.56 

.00098 

.00276 

2.20000 

?3 

.1604? 

6.43 

0.40 

.0022835 

54  0,0.) 

0.00000 

0.00000 

2.20000 

23 

.1604? 

8 .57 

o.oo 

.C02283S 

840.00 

0. 00O00 

0.00000 

2.2C000 

?3 

.1604? 

10.71 

0 . 30 

.-30??H35 

540.0.) 

0.00000 

0.00000 

2.20000 

23 

.1604? 

12.86 

O.K 

. 0022835 

540.00 

y. ooooo 

0.00000 

2.20000 

23 

.  1604? 

15.04 

0.  OC 

.00? 2835 

54  0,00 

0.00000 

0.00000 

2.20000 

?3 

.1604? 

17.14 

o.r-o 

.0022835 

540.00 

0.00000 

0.00000 

2.20000 

23 

.1604? 

10.29 

O.'OC 

.0022835 

54  0.00 

0.00000 

0.00000 

2.20000 

?3 

.1604? 

21  .43 

0.30 

.0022835 

540.09 

o.oocyo 

y  .ooooo 

2.20000 

23 

.1604? 

23.57 

0.00 

.•<0228  35 

540.00 

0.00000 

o. ouooo 

2.20090 

?3 

.1604? 

?5 . 7  1 

O.rr 

.  c  022835 

84  r  .00 

0.00000 

0 .00000 

2.20000 

23 

.1604? 

?7 , 86 

0 . 0  o 

.6022835 

540.00 

0. 00O00 

0.00000 

2.20000 

?3 

.1604? 

30.00 

O.'.O 

.0022835 

54  (,.00 

0.00000 

0.00000 

2.20000 
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CIRCUMFERENTIALLY- AVER AGED  PARAMETERS 


NR 

RADIUS 

MACH  NO. 

TEMP. 

INTENSITY 

FREQUENCY 

1 

.0001 

1 • 9882 

2.2955 

. 669O6E- 1 2 

0. 

2 

.1000 

1.9662 

2.2955 

.28773E-06 

3. 

3 

.2000 

l .9662 

2.2955 

. 954 1 6E-0H 

C  . 

A 

.3000 

1 .9668 

2.2951 

.21384E-04 

7. 

S 

.8000 

1  .9182 

2.2398 

.161038*05 

2816. 

6 

.5000 

1.2278 

1  .9688 

. 1 7 1 98E ♦ 08 

25958. 

7 

.8000 

.  1809 

1.3852 

.510438*04 

16407. 

8 

.7000 

.0028 

1.0583 

.39885E-04 

4377. 

9 

.8000 

.0019 

1  .0  777 

.41234E-06 

524. 

10 

.9000 

.1277 

1.2755 

.1  ">  3  0  4  E  *04 

5139. 

11 

1 .0000 

.  8  929 

1.8257 

•  1  1 824E ♦ 08 

12312. 

1? 

1.1000 

1  .5573 

2.1133 

. 2°0 1 6E ♦ 08 

5707. 

13 

1  .2000 

1.7232 

2.2018 

.4221 3F*08 

2595. 

18 

1.3000 

1 .7788 

2.2287 

.460085*08 

1  676. 

15 

1 .8000 

1.7968 

2.2343 

.415825*08 

1326. 

16 

1.8000 

1 .8028 

2.2367 

.8291  RE *08 

1313. 

17 

1.8000 

1.7963 

2.2339 

.  386?4F*08 

1513. 

18 

1.7000 

1.7761 

2.2255 

•  4  354  3F  *08 

1888. 

19 

1.8000 

1.7236 

2.2018 

•  3C032F  *08 

2593. 

?0 

1.9000 

1 .5589 

2.1139 

•24184E*08 

4755. 

21 

2.0000 

.8929 

1  .8257 

. 14303E*08 

8933. 

2? 

2. 1000 

.1275 

1 .2750 

. 1S086E*08 

3436. 

23 

2.2000 

.0029 

1  .0  397 

.19514E-04 

733. 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATER  than 

2  AT 

KA  = 

1 

X  = 

.07292 

ITH=  10  TMFT  A  = 

no.oc 

NTP  = 

3 

warning 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATER  THAN 

2  AT 

ka= 

1 

X  = 

.07292 

ITH=  11  T  HE  T  A  = 

120.00 

NTP  = 

3 

warning 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATER  THAN 

2  AT 

KA  = 

1 

X  = 

.0729? 

ITH=  12  T  Hf  T  A  = 

1 30.00 

NTP  = 

3 

W  ARN I  NO 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATER  THAN 

2  At 

KA  = 

1 

X  = 

.0729? 

ITH=  13  THE  T  A  = 

140.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATFD  THAN 

2  AT 

KA  = 

1 

X  = 

.07292 

ITH=  14  THET  A  = 

150.00 

NTP  = 

3 

WARNING 

-  NO. 

OF 

TURNING 

POINTS  IS  GREATER  THAN 

2  AT 

KA  = 

1 

X  = 

.0729? 

I T  H=  15  THE  T  A  = 

160.00 

NTP  = 

3 
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X  <  1 1  = 

WARNING 
WA  =  2 


earning 

KA=  2 


.V  ARM  I  MG 
K  A=  2 


FARMING. 
*A=  2 


WARMING 

ka=  .2- 


x  (  ?)  = 

WARNING—' 
KA=  3 


warming 
KA  =  _„3_ 


warning 

KA  =  3 


WARNING 
K  A  -  3 


WAPNING 
«A  =  .3 


*(  3)  = 

WARNING 
KA  =  4  - 


WARNING 
WA*  4 


.0729  llfl I (  1)=  .?5899E*?1  FM<  1)=  .24)9E*01  UAVM  1)=  1847.59  UMAX (  1)=  2199.46 

NO.  OF  TURNING  POINTS  IS  GPFATFP  THAN  2  AT 
X®  .09187  1TH=  10  THFTA=  110.0C  MTP=  3  -  - - - - - 


NO.  .OF  TURNING  POINTS  IS  GREATER  THAN  2  AT 
X=  .09187  ITH=  11  THE  T  A=  120.00  NTP*  3 


NO.  OF  TURNING  POINTS  IS  GPFATFP  THAN  2  AT 
X  =  .09187  1 TH  =  12  THETA®  130.00  NTP.® 3 


NO  ._0E  TUPNIM6-  POINTS  IS  OPE  A  TER  THAN..?  AT  ... _ 

X=  .09187  I T H=  13  THE  T A®  140.00  NTP®  3 


MO.  OF  TURNING  POINTS  IS  GREATER  THAN  2  AT 
.  X=_  09187- . _1  T.H®.  18 _  tut TA=...  150.00  _ _ NIELS _ 1 


.0919  It8j(  21=  .?(38Sf«21  EM<  2  1  =  .?H?9f.01  UAVG  (—24= — li/»7^S8 - UMAX! 


NO.  OF  TURNING..  POINTS  IS  GREATER  THAN  2  AT 
X=  .  I  1 R-7S  1TH=  10  T  HE  1  A=  ll'i.CV 


NTP® 


NO.  OF  TURNING  POINTS  IS  GREATER  THAN  2  AT 
_ X  =.  ...  ..JL1575 _ ITH=.Jtl _ J.UEl  A.® I ?v .00 _ Ell£5 1 


NO.  OF  TURNING  POINTS  IS  GREATER  THAN  ?  AT  ......  _ 

X  =  . 1 1 S7S  ITH=  1?  T  HE  T  A=  1  30.00  NTP=  3 


NO.  OF  TURNING  POINTS  IS  GREATER  THAN  2  AT 
-X=  .11575 _  1TH=  13  THETA®  140.00 _ NTP-*  —3.. 


NO.  OF  TURNING  POINTS  IS  GREATER  THAN  2  AT  .  _ 

X*  .11575  IT  H«  14  THETA®  150.00  NTP®  3 


.1157  1)81  (  31=  .2*1  9?F  *21  FM  (  31=  ,?669F.*01  UAVGI  31=  1674.56  UMAX  (  31=  2199.45 

NO.  OF  TURNING  POINTS  IS  GREATER  THAN  ’  AT 
x=-  ,14583..  ITH=  10  —  THETA®  11C.00  NTP® _ 3 _ 


.NO.  OF  TURNING  POINTS.  IS. GREATER  THAN  2  AT  .  _  _ _ _ 

X«  • 1 4583  ITH=  U  THETA®  120.00  NTP®  3 
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WARNING  -  NO.  Of"  TURNING  O0INTS  IS  GREATER  THAN  ?  AT 

KA=  A  X  =  .14S83  I TH-  1?  THE  T  A  =  1  30.00  NTP  = 
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COMPUTATION  of  AERO-ACOuST IC  PROPERTIES  Of  SUPPRESSOR  NOZZLES 


CASE  NO.  1  CRO  7  —  TORE  AR  =  ?,3  NOZZLE  -  VJ=2?00  f PS  -  TTJ=1600  nEG-R 


AXIAL  LOCATION  =  .291ft7  (X/OEO  =  4.00000 


M 

P 

ANGLE 

II 

Of  NS  ITT 

TEMP. 

U/UPEE 

TURR.INT. 

R/DEO 

1 

.00001 

,'.C0 

2106.11 

.0010468 

1177.90 

.95756 

.00168 

.00010 

1 

.00001 

13.00 

2106.11 

.3010468 

1  177.9o 

.95756 

.00168 

.000)0 

1 

.A0001 

20.00 

2106.11 

.0010468 

1 1 77,99 

.95756 

.00  168 

.00010 

1 

.00001 

30.00 

2106.11 

.0010468 

117  7.9,' 

.95756 

.00168 

.00010 

2 

.0  0  729 

o.co 

2079.70 

.r  010551 

1  168.7  3 

.94555 

.05657 

.10000 

2 

.0  0729 

1  0.00 

2079.7? 

.  ■'010550 

1  168.73 

.94556 

.  05658 

.  10000 

2 

,00729 

20.00 

?n  79.7? 

.3010550 

I  168.73 

.94556 

.  05658 

.  inono 

2 

.00729 

30.  CO 

2)79.70 

.0010551 

1166.73 

.94555 

. 05658 

.  loooo 

3 

. 0 1458 

0.00 

1994.78 

.3010776 

1  144.29 

,9f  694 

.08729 

.20000 

3 

.0 1480 

7.50 

1 494 . 74 

.0010776 

1  144.27 

.9069? 

.08729 

.20000 

3 

.31488 

15.00 

1994.77 

.0010776 

1  144.28 

.90694 

.08730 

.20000 

3 

.01458 

22.50 

1994.73 

.'010776 

1144.27 

.9069? 

.08730 

.20000 

3 

.014SR 

30.00 

1994.77 

.0010776 

1144.28 

.90694 

.08730 

.20000 

4 

•  0  2 1 88 

0.00 

1838.10 

.0011099 

1111,01 

.83571 

.11429 

.30000 

4 

."2188 

7.80 

1838.  1  3 

.0011099 

1111,0? 

.8357? 

.1 1430 

.30000 

4 

.  0  ?  1  8  a 

15.03 

1838.04 

."011099 

1110.97 

.83568 

. 1 1431 

.30000 

4 

.  0  ?  1  8  8 

22.50 

1838.06 

.001 1099 

1110.95 

.83569 

. 1 1432 

.30000 

4 

•  C  ?  1  8  8 

33.00 

1837.99 

.no  111  00 

1110.92 

.83566 

.1143? 

.30000 

5 

.02917 

0.00 

1601.57 

.<•011489 

1073.31 

.72817 

.13158 

.40000 

5 

.  C29 1  7 

6.00 

1601.44 

.9011489 

1 " 7  3. 23 

.72811 

.1 3160 

.40000 

5 

.02917 

12.00 

1601.33 

."'011491 

1 373.07 

.72804 

.13166 

.40000 

5 

.02917 

18.00 

1601."? 

.0011493 

1372.85 

.7279? 

.13172 

.40000 

5 

.0291 7 

24.03 

1600.73 

.0011496 

1  372.64 

.72778 

.13176 

.40000 

5 

.32917 

30.00 

1600.68 

.0011496 

107?. 58 

.72776 

.13178 

.40000 

6 

.  03ft4ft 

0.00 

12Q8.70 

.0011959 

1031.04 

.59046 

.13285 

.50000 

ft 

.0364ft 

ft. 00 

1298.61 

.0011963 

1030.74 

.59043 

.13268 

.  5  0  0  o  0 

6 

. 0  364 ft 

12.00 

l?9ft.97 

.0011975 

1  29.74 

.58968 

.13295 

.50000 

6 

."3646 

18. CO 

1295.36 

.3011988 

1028.56 

.58895 

. 1 3329 

.50000 

ft 

.03646 

24.00 

1294.41 

.0011949 

1 027.64 

.58851 

.13355 

.50000 

ft 

. 0364ft 

30.00 

1293.49 

.0012004 

1 C27.20 

,58810 

.  1  3363 

.50000 

7 

.04378 

0.00 

980.18 

.0012487 

987,46 

.44565 

.11324 

.60000 

7 

.04375 

8.00 

978.21 

.0012502 

986.28 

.44475 

.11362 

.60000 

7 

.04375 

10.00 

973.63 

.1012537 

983.54 

.44267 

.  1  146? 

.60000 

7 

.04378 

15.0C 

967.62 

.0012584 

979.85 

.43994 

.11590 

.60000 

7 

. 04375 

20.no 

961.38 

.0012636 

975.87 

.43710 

.11703 

.60000 

7 

.04375 

25.00 

956.95 

.0012673 

97?. 99 

.43508 

.11782 

.60000 

7 

.34375 

39.05 

955.55 

.C012665 

97?.08 

.43445 

.1181? 

.60000 

8 

.05104 

0.00 

758.09 

.0012799 

963.41 

.34467 

.06822 

.70000 

8 

.05104 

5.00 

751.56 

.0012834 

9fto.81 

.34170 

.07114 

.700O0 

8 

.05104 

10.00 

733.93 

.0012928 

953. 8p 

.33369 

.07705 

.70000 

8 

.05104 

15.00 

710.08 

.0013064 

943.84 

.32285 

.0825? 

.70000 

8 

.05104 

20.00 

68ft. OR 

.3013216 

933.04 

.31193 

.08631 

.70000 

8 

.05104 

25.00 

668.47 

.3013335 

9?4.7? 

.30393 

.08843 

.70000 

a 

.05104 

30.00 

662.1 7 

.9013379 

9? 1 .68 

.30101 

.1)891  1 

.70000 

9 

.05833 

0.00 

758.05 

.0012799 

963.41 

. 34m65 

.06817 

.80000 

9 

.C5833 

4.29 

745.61 

.0012839 

960.44 

.33900 

.07124 

.80000 

9 

.05833 

8.57 

710.47 

.0012950 

95?. 1ft 

.3230? 

.07635 

.00000 

9 

.05833 

1  2.8ft 

658. 39 

.0013121 

939.74 

.29934 

.07855 

.80000 
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COMPtlT  AT  t  OP  AfRO-ACOllST  TC  PROPERTIES  OF  SUPPRESSOR  NOZZLES 


CASE  NO.  1  CRH  7-TMHp  AR=2.3  NOZZLE  -  VJ=2?00  EPS  -  TTJ=160G  DEG-W 


AXIAL  LOCATION  =  .29167  (X/OEQ  =  A.OOOO-il 


M 

R 

ANGLE 

U 

DE NS  ITT 

T  F  MP  . 

U/UREF 

TURH. INT 

.  R/OEQ 

9 

. ASA  33 

17.14 

597.68 

.(01 3330 

975.91 

.77163 

.0758? 

.80010 

9 

.05833 

71.43 

540.35 

•  CO  1 3531 

511.79 

.24567 

.16819 

.60000 

9 

.05833 

75.71 

498.87 

.0013676 

9(1.61 

.7768? 

.(5713 

.60000 

9 

.(5833 

30.00 

483.69 

.0013778 

698. 1 9 

.21991 

.05027 

.80000 

10 

. 16863 

0.00 

980.09 

.0017467 

587.45 

.4456 1 

.11371 

.90000 

10 

.06563 

4.79 

958.94 

.0017541 

983.75 

.43599 

.11533 

.90000 

10 

. C  656  3 

8.57 

898,45 

.901769? 

971.57 

.40849 

.  1  1  89? 

.90000 

10 

,06563 

12.86 

805.40 

.0017938 

953.09 

.36618 

.  1  1906 

.90000 

10 

. 16563 

17.14 

692,79 

.901 3744 

931.01 

.31498 

. 1 1248 

.90000 

10 

.06563 

71  .43 

579.98 

. 9013541 

51  .63 

.76369 

.09745 

.40000 

10 

.06563 

75.71 

49?.43 

.0013716 

895.03 

.72389 

.(7158 

.90000 

10 

.06563 

30.00 

458.84 

.101 3756 

896.38 

.70861 

.03069 

.90000 

11 

.17797 

3.90 

1798.64 

.0011960 

1031.04 

.59044 

.13253 

1.10000 

11 

.0779? 

3.75 

1778.97 

.0011997 

1077.81 

.58150 

.  1  3588 

1 .00000 

11 

.0779? 

7.50 

1713.17 

.0017170 

1C  17.4? 

.55158 

.14787 

1  .00000 

1 1 

.  0729? 

1  1  .25 

1109.8? 

.9012334 

999.73 

.50459 

.  14800 

1  .00000 

1 1 

.0779? 

15. CC 

977.73 

.0017630 

976.3? 

.44453 

.14777 

1.00000 

1 1 

.0779? 

18.75 

879.45 

.0017980 

950. CO 

.37717 

.13834 

1  .00000 

11 

.0729? 

77.50 

686.90 

.101 3  3  C  q 

976.59 

.31190 

.11965 

1.10000 

1 1 

.1779? 

76.25 

576.0? 

.0013596 

912.95 

.76139 

.0888? 

1.(1000 

1 1 

.0729? 

33.00 

533.54 

.<*01  3553 

909.80 

.74750 

.04985 

1  .00000 

1? 

.08071 

0.00 

1601.50 

.0011489 

1173.30 

.7281 3 

.13157 

1.10000 

17 

.08071 

3.75 

1577.0? 

.0011578 

I  ',64.60 

. 7 1 7u  1 

. I  3788 

1 .1  0000 

1  7 

.C8C71 

7.50 

1503.70 

.0011648 

1058.63 

.68367 

.15133 

1  .  1COOO 

I? 

.18171 

1  1  .75 

1383.35 

.001 1«48 

10 41.71 

.62895 

.16335 

1  .  1  0000 

17 

.98071 

15.00 

1719.6? 

.0017135 

i 1 6 .  n 

.55451 

.  16798 

1  .  1  0000 

17 

.98071 

18.75 

1078.14 

.0017597 

985.89 

.46745 

.16135 

1 .10000 

1? 

.08071 

77.50 

838.93 

.0017885 

956.97 

.3810? 

.14139 

1  .  1  0000 

1? 

.08021 

76.75 

688.35 

.0013137 

938.66 

.31296 

. 1 0494 

1  .  1  0000 

1? 

.08071 

30. OC 

679.41 

.0013701 

934.09 

.78617 

.05050 

1 .10000 

n 

. 98750 

t.OO 

1«38.04 

.0011 999 

n  n  .01 

.83568 

. 1 1478 

1 .70000 

n 

.98750 

3.33 

1817.98 

.0011135 

1  10  7.37 

.87656 

.17296 

1 .71000 

13 

. j8  750 

6.67 

1756.96 

.001)741 

1096.90 

.7988? 

.14180 

1 .70000 

13 

.08750 

1  0.00 

1653.43 

.901141? 

1  080.52 

.751 74 

.16110 

1.20000 

1  3 

.08750 

1  3.33 

1507.73 

.901 164? 

1059.16 

.68550 

.17533 

1.20000 

13 

.08750 

16.67 

1 376.8? 

.001193? 

1033.45 

.60325 

.18014 

1.20000 

13 

.08780 

7  9.00 

1  116.83 

.0017787 

1003.54 

.50778 

. 1 7740 

1 ,20000 

13 

.18750 

73.33 

915.38 

.0017641 

975.46 

.41619 

.15031 

1.20000 

13 

.98750 

76.67 

760.19 

.0017869 

958.?:, 

.34563 

.11(59 

1  .70000 

13 

. 08750 

30.00 

699.78 

.0017976 

953.93 

.31816 

.0394? 

1  .20000 

14 

.09479 

3.00 

1994. 74 

.0010776 

I 144.78 

.90693 

.08779 

1  . 3«000 

14 

.0,9479 

3.33 

1974.99 

.001082? 

1139.40 

.89795 

.10155 

1 .30000 

14 

, o 94  79 

6.67 

1914.9? 

.0019953 

1  175.76 

.87073 

.  17855 

1  .  30000 

14 

. 09479 

1  0.00 

18  07 .88 

.0011154 

1  H5.5? 

.82182 

.15554 

1.30000 

14 

.09479 

1  3.33 

1652.47 

.0011410 

1680.71 

.75131 

.  1  76h1 

1  .  30010 

14 

.0 9479 

16.67 

1451 .76 

.1011771 

1057.06 

.66005 

.  18671 

1  .30100 

14 

.09479 

70.00 

1718.1? 

.0017089 

1019.96 

.55383 

.  18708 

1 .30000 

1  4 

. 994  7  9 

7  1.33 

907.56 

•00174HS 

987.61 

.  i.i.67  3 

.  1  5998 

1  .  30000 

78 
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COMPUTATION  OP  AFWO-ACOlJST  K  PPOPEBtIES  OK  SUPPPFSSOP  N07ZLFS 


CASE  NO. 


CPn  7-TUHE  AP=2.3  N077LF 


V J=22G0  EPS  ~  TTJ=lf.OO  DEP 


AXIAL  LOCATION  =  .29167  (X/OEO  =  4.000C'1) 


M 

P 

ANGLE 

U 

14 

.09479 

26.67 

798.20 

14 

.06479 

30.00 

725.24 

15 

. 1 0208 

C  .  0  0 

2  C  7  9 . 6  8 

15 

•102Q8 

3.00 

2064.03 

15 

.10208 

6.00 

2015.26 

15 

. 1 0208 

9  .  U  0 

1928.64 

15 

. 10208 

12.00 

1 799.20 

15 

.10208 

15.00 

1624.99 

16 

.10208 

18.00 

1411.64 

15 

.10208 

21.00 

1172.55 

16 

.10208 

24  .or, 

942.79 

15 

.10208 

27.00 

766.69 

15 

. 1 0208 

3  >.00 

698.15 

16 

. 1 0938 

t.  .00 

2U6.C8 

16 

.10938 

3 ,  v  C 

2089.0? 

16 

.10938 

6.00 

2039.21 

16 

.10938 

9.00 

1948.1? 

16 

.10938 

12.00 

1810.36 

16 

.  1  0938 

1 5 . 0  * 

1623.35 

16 

.10938 

18. OC 

1  393.64 

16 

.10938 

21 .00 

1  135.88 

16 

.10938 

24.00 

wflH.64 

16 

.10938 

27.00 

698.70 

16 

.10938 

30.00 

624.17 

17 

.1 1667 

0.00 

2079.64 

17 

.11667 

2.73 

2064.86 

17 

.1 1667 

5.95 

2018.68 

17 

.11667 

-.18 

1937.3? 

17 

. 1 1667 

13.91 

1815.35 

17 

. 1 1667 

13.64 

1650.35 

17 

.11667 

16.36 

169S.39 

17 

.11667 

19.09 

1210.33 

17 

. 1 1667 

?1.«? 

965.46 

17 

. 1 1687 

28.55 

743.54 

17 

.11667 

27.27 

578.88 

IT 

.1 1667 

34. OC 

515.81 

IP 

.  12306 

O.OC 

1 994.58 

18 

.12396, 

2.73 

1977.29 

IP 

.12396 

5.45 

1  9  24 .4, 1 

18 

.12396 

4.18 

1831.25 

IP 

.12396 

10.91 

1695. 74 

IP 

.1 ?396 

13.64 

1516.33 

IP 

. 1 ?  396 

16.36 

l  702.89 

IP 

.12396 

19.09 

1057.75 

IP 

.12396 

21.82 

818.59 

IP 

.12396 

?4 .55 

606.74 

IP 

.12396 

27. 27 

452.0? 

18 

. 1 2396 

3,1.00 

39? .4? 

density 

TFMP. 

.0012757 

966 .58 

.0012871 

960.98 

.0010551 

1  16H.7? 

.0010596 

1  1  63.69 

.0010726 

1  149.60 

.(,010923 

U?8.86 

.•>011168 

1  IO.OP 

.0011454 

1 076.54 

.0011788 

1046.01 

.1012181 

1012.27 

.,01259? 

979.23 

.'>01287? 

967.97 

.0012947 

952.81 

.(010469 

1 177.88 

. 0010520 

1 172.16 

.0010665 

1 156.23 

.(01088? 

1137.09 

.0011161 

1 1 05,8  0 

.  >-.01  1464 

) "75.63 

.0011837 

1941,70 

.0012300 

1 002,50 

.0 012808 

962,75 

.001 3180 

975.53 

.0013288 

927,94 

.0010551 

1 168,69 

.0010594 

1167.93 

.0010717 

1160.53 

.0010906 

1  1 70,61 

.0 01 1180 

1  106.47 

.001 1426 

1079.16 

.001 1768 

1  )4  7 ,8? 

•  7  0 1 2  ?  0 1 

1  010.61 

.001275? 

966,97 

.0013331 

9?4 ,94 

.0013749 

89m, 87 

.0013875 

688.71 

.0010777 

1144,15 

.0010810 

1179.85 

.0010975 

1127.67 

.'■011116 

1  109.29 

.2011 354 

ll.P6.06 

.0011664 

1 ‘58.07 

.0012047 

1 *23.87 

.001258? 

979.99 

. 101 726? 

9?9, 78 

.00)3991 

881.3? 

.0019540 

848,04 

.0014709 

M  38 . 70 

U/UPEF  TDPH.IM.  P/DEO 


36291 
3297  3 
94554 
93843 
91626 
B7  687 
8  1  8  C  2 
7368  1 
64181 
,5331 1 
42865 
,  14858 
,3174? 
,95755 
,95)20 
.92714 
,88573 
.82310 
.  7  3  6 1>  7 
.63363 
.51644 
.40893 
.31767 
.28378 
.94553 
.93881 
.91790 
.8808? 
.82536 
.75034 
.65716 
.55029 
.43895 
.33806 
.26319 
.23434 
.  9  (•  6  8  5 
,89899 
.87477 
.83259 
.77080 
.68941 
.59237 
.48091 
.37218 
.27586 
.20551 
. 1 784? 


.11775 

.  00465 
.05655 
.07855 
.10976 
.14000 
.16625 
.18425 
.19010 
.18135 
.15690 

. 1 146? 

.0400? 
.00131 
.06982 
.10518 
.1377? 
.16571 
, 1 8464 
.18980 
.1796? 
.15401 
.1121? 
.05158 
.05661 
.07639 
.1  0664 
.13436 
.15985 
. 1 7  m  34 
. 1 8636 
.  1  8208 
.  16525 
.  1  3778 
.09931 
.05403 
.08734 
.09823 
.12066 
. 14455 
.16479 
.17735 
.17931 
.16939 
.149)8 
.12107 
.08895 
. 0497? 


1 .30000 
1  .30000 

1 .40000 

1 .40000 
1  .  4  f.  n  l>  0 
1 .4  CO  60 

1.40000 

1 .40000 

1 .40000 
1  .  4  o  0  0  0 

1 .40000 

1.40000 
1.40000 

1.50000 
1  .50000 
1  .50000 
1  .50000 

1 .50000 
1  .50000 

1 .50000 
1  .SoOOO 

1 .50000 
1  .5000 j 
1  .50000 
I  .  6  0  0  0  0 

1.60000 
1  .60000 
1  .60000 
1  .60000 
1  .60000 
1  .60000 
1  . 6  0  0  C  0 

1 .60000 
1  .60000 
1  .60000 

1 .60000 
1  .70000 

1 .70000 

1 .70000 

1 .70000 

1 .70000 

1 .70000 

1 .70000 

1.70000 

1 .70000 

1.70000 
1  .70000 
1  .  70(100 


ft  o  ft 


M 


ft  ft  ft 
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r,  K 

fO^PPTAT  !OkJ  OP  AFriO- ACOUSTIC  PROPFRT1FS  OF  SlJPPRfSSOP  NOZZLES 
CAST  mo.  1  CPU  7-TIIRf  A  R  =  ?  .  3  NOZZLE  -  VJ=2?0G  FPS  -  TTJ=1600  PEG-R 


AXIAL  LOCATION  =  .39167  (x/PEO  =  4.000UO) 


M 

i: 

ANGLE 

U 

no  N$1H 

temp. 

ll/UREF 

TU9H.lt-  T 

.  H/6fQ 

IP 

.  niPR 

'.00 

1837.39 

. 001 1108 

1 1 1  '.4? 

.63539 

.1 1439 

1.60000 

1  9 

.1 31?P 

P  .  5  8 

1HPJ.44 

.011136 

1  1(  7.37 

.  8?  76H 

.  1200  7 

1  .6(1000 

ip 

•  I  TIPS 

8.00 

1 768.89 

.00 1 1??8 

1  nop,?? 

•  6  C  424 

.  1  3359 

1  .  HO000 

IP 

.1  TIPS 

7.9: 

1661.4? 

.001 1 378 

K-c  3.77 

.  7 6447 

.14907 

1  .80000 

IP 

. 1 71  PS 

k.  or 

1556.9? 

.001 1547 

1  664.21 

.  7  C  7  8  7 

.16204 

1 . 6  0  0  0 

IP 

.  131P5 

IP.  80 

1  3r'7.48 

. C  0 1 1866 

1  >  9  V  .  1 4 

.63537 

.16931 

1 .80000 

ip 

•  1  3 1 ?5 

15.0.7 

1  P  O  6 . 7  0 

.O01PP45 

1  '06.97 

.54864 

.  1686  8 

1 . 8  a  c  0  0 

IP 

. 1  TIPS 

17.50 

997.54 

.00) P 764 

966,03 

.45354 

.15935 

1  .60000 

IP 

. 1 31  PS 

P'  .00 

787.4? 

.00174/9 

91 m,?4 

.35819 

.  14P4  9 

1.80000 

IP 

.1 71?P 

PP.50 

597.76 

.  0  C  1  4  ?  1  7 

n67, 

.26996 

.12005 

1.60060 

IP 

. 1 31PP 

3  5. 07 

471 .?? 

.0015017 

8P 1  .  1  2 

. 1 9606 

.  09427 

1.80000 

IP 

. 131  PS 

P  7.5C 

316.89 

•  C  0  1  c  5  9  5 

74:1 ,68 

.14408 

.06595 

1 .84000 

IP 

.1 31?c 

7 1; .  0  0 

P73.44 

.3015771 

751,86 

.  124  33 

. 04087 

1 .80000 

?3 

. 1 3854 

7.00 

1598.74 

.:  0  n  8 1 6 

1.79,72 

.  7?688 

.13194 

1.90000 

?c 

. 1  3854 

P.90 

1  579.89 

.501 1568 

1  067.79 

.71831 

.135)2 

1 .90000 

PO 

. 1 3«S4 

8.0C 

15?3.?5 

,001 1 644 

1 1,58.99 

.69256 

.14277 

1.9O000 

PO 

. 1 3854 

7. 50 

1  4?9.  ?  7 

,701  l  8  8 

1 944.23 

.64483 

.15099 

1  .90000 

PO 

. 1 3894 

1  .CO 

1701.45 

. 7  0 1  PC  8? 

1323.13 

.55190 

.1562? 

1 .9(1000 

PO 

. 1 3«84 

IP. 50 

1 1 36.77 

. CO  124 0  7 

993.84 

.5)664 

.15577 

1 . 9  7  0  0  0 

PO 

. 1 3694 

15.00 

954.54 

.  0  0  1  ?  8  9  6 

956,17 

.43349 

. 14844 

1.46000 

PO 

. 1 3«94 

17.5  ' 

766.  15 

. ;01 3837 

911.36 

.34834 

. 1 3488 

1 . 40000 

?o 

. 1 3«S4 

Pf.OO 

545.1-  8 

. 00 143PC 

861.07 

.  2 6  6  C  1 

. 1 1 549 

1  .90000 

PC 

. 1 3894 

PP.57 

4  ?5 . 4  ? 

.  7  0 1 5  ?  3  n 

8  0,65 

.1934? 

. 0935? 

1.90000 

PO 

. 1 3894 

PS.l'C 

?96.dS 

.001613? 

764.30 

.  1  3496 

.07079 

1 .40000 

PO 

. 1 7894 

P7 . 5  7 

P08.C  1 

.7016788 

734.48 

. ' 9457 

. 04  8  p6 

1 .9001  0 

PC 

. 1 7854 

30. 0C 

1 74.P6 

.'016987 

7P5.87 

.07923 

.03019 

1  .90000 

?i 

.14583 

C  .00 

1 ?4  7  «  46 

."012079 

1  '  2  7.84 

.58535 

.13417 

2.50000 

?i 

. 14587 

P.31 

1P73.P3 

.0012109 

1  18.35 

.57888 

.  1353? 

2. 00 000 

?i 

.  1  4587 

4.63 

1  P?1 .46 

.701PP16 

n  09.39 

.55538 

.13771 

P. 00000 

?i 

.14583 

6.PP 

1 140. ?c 

.001 ?4  0  0 

994.44 

.51840 

. 1 3956 

P.00000 

?i 

.14543 

P.P  1 

1033.50 

.  -0  IP67? 

473,08 

.46989 

. 1 3924 

P.00000 

?i 

.14583 

1  1  .54 

905.48 

.001 3048 

944.99 

.41169 

. 1 3506 

P.00000 

pi 

.14587 

)  3.85 

764.63 

. " 01 7539 

917.74 

.34764 

.  12645 

p.ooooo 

?i 

. 14587 

16.15 

619.74 

.0014161 

87c .74 

.28177 

.11370 

P.00000 

?i 

. 14587 

18.46 

481.43 

.001491 1 

626.87 

.P] 889 

.09777 

P.OOOOO 

?i 

. 14583 

P0.77 

357.91 

.0015764 

78 2.20 

.16273 

.060  12 

P.OOOOO 

?i 

.  14587 

P  7 . 0  8 

P54.87 

.0(116664 

7  34.96 

.11574 

.06236 

P.OOOOO 

?i 

.  14543 

?5. 38 

175.01 

.O0175C 1 

7  0  7*  .  5  8 

.07957 

•  04584 

P.OOOOO 

?i 

.14587 

P  7.68 

1  P 1 . 64 

.7018779 

68?. 01 

.05530 

.0309? 

P.OOOOO 

?i 

.14587 

37.00 

101 .49 

.  >01 H?49 

675.70 

.  04614 

.020  1  1 

P.OOOOO 

?? 

.15317 

0 . 0  0 

937.77 

.' 01 ?946 

95? .48 

.42637 

.1  1983 

?  .  1  0  0  0  0 

?? 

.15317 

P.14 

9  P  4  •  9  5 

.  0  0  1  ?  9  8  7 

944.49 

.42053 

.  11979 

P. 1 0000 

?? 

.15317 

4.PP 

887.81 

. :  0 1 3 1 7  t 

941  ,C4 

.40  365 

.11943 

P.  1 0000 

p? 

.18317 

6.43 

8PH.O8 

."01 7308 

9?6.59 

.37649 

.11798 

2. 1 0000 

?? 

.157)7 

4.57 

749.54 

.001 7598 

906.8/) 

.34079 

.  1  1460 

P  .  1  0  0  0  0 

?? 

.15313 

17.71 

656.99 

.7013990 

661.41 

.29670 

.10870 

P. 1 000 0 

?? 

.1531 3 

IP. 86 

556.37 

.0014486 

651.23 

.28296 
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COMPIITAT ION  Of  AFRO-ACOUSTIC  PROPEPTIFS  OF  SUPPRESSOR  N07ZLL S 
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DAGE  1 7 


COMPUTATION  OF  AFRO-ACOUSTIC  PROPERTIES  OK  SUPPRESSOR  NOZZLES 


CASE  NO', 


CPD  7-TUBE  AR=2.3  NOZZLE  -  VJ=2200  EPS  -  TTJ=1600  DEG-R 


AXIAL  LOCATION  =  1,16667  (X/OEQ  =  15.99049) 


M 


R 


ANGLE  U 


DENSITY  TEMP, 


U/UREF  TURB.INT.  R/DE(j 


1 

1 

1 


3 

3 

3 

3 

4 
4 
4 
4 

4 

5 
5 

5 

5 

5 

5 

6 
6 
6 
6 
6 
6 
7 
7 
7 
7 
7 
7 

7 

8 
8 
a 
8 
a 
a 
a 
9 
9 
9 
9 


.COCCI 

0.3C 

1238.43 

.''012299 

.00001 

1  (■ .  0  0 

1238.43 

.  -012299 

.00001 

23.00 

1238.43 

, 0  C 1 2299 

.00001 

30.00 

1238.43 

.0012299 

.01167 

0.00 

1236.43 

.0012310 

.01167 

10.03 

1236.43 

.0012310 

.01 167 

20.30 

1236.44 

.0012310 

.01167 

30.00 

1236.44 

.0012310 

.02333 

O.Ol) 

1230.9c 

.0012345 

.*2333 

7.5o 

1230.89 

.0012345 

.02333 

IS. CO 

1230.90 

.0012345 

.02333 

22.50 

1230.42 

.0012345 

.02333 

31.03 

1230.98 

.0012344 

.03500 

0.00 

1221.34 

. uO 1 24C5 

.03500 

7.50 

1221.47 

.0012404 

.33500 

15. CO 

1221.34 

.0012404 

•  03503 

22.50 

1221.41 

.•*0124  04 

•  0350C 

30.  OC 

1221.32 

.0012404 

.04667 

0.00 

1206.88 

.0012493 

. n4667 

6.00 

1206.94 

.0012493 

.04667 

12.00 

1206.98 

,0012493 

.04667 

16,0-3 

1206.86 

, CO l 249? 

.04667 

24.00 

1206.74 

.001249? 

.04667 

30.  CO 

1206.84 

.001249? 

.05833 

6. 00 

1 187.6? 

.001261 3 

.0.58  33 

6.00 

1187.62 

.oni?6i? 

•  05833 

12.00 

1187.26 

.001261 1 

. C  58  33 

18.CC 

1187.  j 6 

.0012610 

.05833 

28.00 

1187.09 

.0012608 

.05833 

30.00 

1  187.07 

.0012608 

•C700C 

0.00 

1162.27 

.0012765 

.07000 

6.00 

1  162.31 

.0012764 

.07030 

H’.OO 

1161.97 

.0012762 

.07000 

15. OC 

1  161.48 

.0012759 

.07000 

20.00 

1161. 28 

.0012755 

.07000 

25.00 

1 160.97 

.0012753 

.07000 

30.00 

1160.9? 

.0012752 

.08167 

0.00 

1130.21 

.0012946 

.08167 

5.00 

1130.C5 

.0012944 

.08167 

10.00 

1129.59 

.0012940 

.08167 

16.00 

1128.75 

.0012935 

.08167 

20.00 

1128.07 

.0012928 

.08167 

25.00 

1127.39 

.0012924 

.08167 

30. 00 

1127. u7 

.0012923 

.09333 

3.00 

1090.28 

.0013155 

.39333 

8.29 

1090.08 

.0013153 

.09333 

8. 57 

1089.36 

.0013147 

.09333 

12.86 

1 088 . 33 

.00)3139 

1  f  0 ? . 6? 

.56306 

.00079 

.00016 

1 302.6? 

.56306 

.00079 

.00016 

1.02.6? 

.56306 

.00079 

.00016 

) 0 02, 6? 

. 563C 6 

.00079 

.00016 

1001.65 

.56215 

.01989 

.16000 

1  CO  1 .65 

.56215 

.01960 

.16000 

1001.66 

.56216 

.01961 

.16000 

1001 .66 

,5621b 

.01960 

.16000 

998.88 

.55964 

.02794 

.32000 

998.87 

.55963 

.02794 

. 32000 

998,88 

.55964 

.02797 

.32000 

998,88 

.55965 

.02797 

.32000 

998.9., 

.55967 

.02798 

.3P0U0 

994. C4 

.55529 

.03476 

.48000 

994,07 

.55535 

.03477 

.46000 

994. C7 

.5553? 

.0348? 

.48000 

994.08 

.55532 

.03489 

.48000 

994.07 

.55528 

.03491 

,48poo 

986.99 

.54872 

.04115 

.64Q00 

987.02 

.54875 

.04119 

.64000 

987.05 

.54876 

.04135 

.64000 

987.07 

.54871 

.04154 

.64000 

987.  in 

.54865 

,04169 

.64000 

987.1? 

.54870 

.04170 

.64000 

977.66 

.53996 

.04768 

.80000 

977.69 

.53996 

.04783 

.80000 

977.75 

.53980 

.04814 

.80000 

977.86 

.53971 

.04850 

.80000 

977.98 

.53972 

.04878 

.80000 

978,03 

.53971 

.04888 

.80000 

966. p 1 

.52843 

. 05460 

.96000 

966,08 

.52845 

.05478 

.96000 

966.24 

.52830 

.05520 

.96000 

966,44 

.52808 

.05571 

.96000 

966,70 

.52799 

.056)8 

.96000 

966.86 

.52784 

.05648 

.96000 

966,94 

.5278? 

.05655 

.96000 

952,48 

.51386 

.0618) 

1.12000 

952.61 

.51378 

.06208 

1.12000 

952.9? 

.51358 

.0627? 

1.12000 

953.32 

.51320 

.06347 

1.12000 

953.79 

.51289 

.06411 

1.12000 

954.10 

.51258 

.06453 

1.12000 

954.20 

.51243 

.06465 

1.12000 

937.34 

.49570 

.06894 

1 .28000 

937.50 

.49561 

.06919 

1.28000 

937.91 

.89529 

.06987 

1.28000 

936.49 

.4948? 

.Ii7f)7? 

1 .28000 

84 
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COMPUTATION  OF  AFRO-ACOUSTIC  PROPERTIES  OF  SUpPRF SSOR  N072LES 


CASE 

NO.  1 

CRD 

7-TtlHF  AR  =  2 . 

3  M07ZLE 

-  VJ=??00  FPS  - 

T T J= 1 600  DEG-R 

AXIAL 

LOCATION 

=  1  . 

16667  (X/DEO  =  15.99999) 

M 

R 

ANGLE 

U 

DENSITY 

TEMP. 

U/UREF 

TURft. 1NT 

.  R/DFO 

9 

.<'9:133 

17.14 

1087.31 

.<.013129 

939.17 

.49436 

.37151 

1  .28000 

9 

.09333 

21  .43 

I086.u4 

.001312? 

939.73 

.49378 

.07208 

1  .28000 

9 

.09333 

25.71 

10«5.25 

.0013115 

940.17 

.4934? 

.0724? 

1  .28000 

9 

.09333 

3  3.00 

1084.88 

.0013113 

940.31 

.49325 

.07253 

1.2P000 

10 

.iosoo 

?.00 

1042.41 

.0013387 

921 .08 

.47394 

.07553 

1  .44000 

10 

.10500 

4.29 

1041.87 

.0013385 

921.26 

.47369 

.07585 

1.64000 

10 

.10500 

8.57 

1040.72 

.0013376 

921.88 

.47317 

.07667 

1  .44000 

10 

.10500 

12.86 

1038.96 

.0013363 

922.76 

.47237 

.07767 

1  .64000 

10 

.10500 

17.14 

1037.13 

.0013348 

923.78 

.47154 

.07854 

1  .44000 

10 

.10500 

21.43 

1 035.17 

.0013336 

924,64 

,47u65 

.07914 

1  .44000 

10 

.10500 

25.71 

1033.91 

.0013326 

925, 3ft 

.97007 

.07946 

1  .64000 

10 

.10500 

37.00 

1 033.26 

.0013323 

925.51 

.46978 

.07954 

1  .64000 

11 

.11667 

0.00 

985.99 

.001364? 

903.86 

.44829 

.08105 

1  .60000 

11 

.1  1*67 

3.75 

985.47 

.0013639 

904.08 

.44805 

.08133 

1  .bOOOO 

11 

. 1 1**7 

7.50 

984.09 

.0013629 

904.73 

.4474? 

.08206 

1  .60000 

11 

. 1 1*67 

11.25 

981.9? 

.0013615 

905.70 

.44644 

.08298 

1  .60000 

11 

. 1 1667 

15.00 

979.34 

.0013598 

9<)6.P4 

.44526 

.08384 

1 .60000 

11 

. 1 1667 

18.75 

976.76 

.0013580 

908.00 

,44409 

.08448 

1  .60000 

11 

.11667 

22.50 

974.34 

.0013565 

908.99 

.44299 

.08488 

1.60000 

11 

.11667 

26.25 

972. 7C 

.0013555 

909.66 

.44225 

.08504 

1  .60000 

11 

.11667 

33.00 

972.31 

.0013551 

909.95 

.44207 

.08510 

1  .600(50 

l? 

.12833 

0.00 

921.76 

.0013917 

8P5.99 

.41909 

.08510 

1.76000 

1? 

.12833 

3.75 

920.97 

.001 3914 

886.24 

.41873 

.08538 

1.76000 

1? 

.12833 

7.50 

918.91 

.('0139C1 

887.0? 

.41779 

.08610 

1.76000 

1? 

.12833 

1  1  .25 

915.85 

.0013883 

888.21 

.41640 

.08701 

1.76000 

1? 

.12833 

15.00 

912.2? 

.0013861 

880,63 

.41475 

.08781 

1.7*000 

1? 

.12833 

18.75 

908. 30 

.0013838 

891.05 

.41297 

.08835 

1.76000 

1? 

.12833 

22.50 

904.86 

.0013820 

892.27 

.41140 

.08861 

1.76000 

1? 

.12833 

26.25 

9C2.56 

.0013807 

893.09 

.41036 

.08866 

1.76000 

12 

.12833 

30.00 

901.71 

.0013802 

893.38 

.40997 

.08866 

1.76000 

13 

.14000 

0.00 

850.11 

.0014216 

867,36 

.38651 

.08730 

1  ,9?000 

13 

.14000 

3.33 

849,37 

.0014212 

867,60 

.38617 

.08750 

1.9?0O0 

13 

.14000 

6.67 

847.30 

.0014201 

868,33 

.38523 

.08805 

1  .92000 

13 

.14000 

10. 00 

843.77 

.0014183 

869.38 

.38363 

.08876 

1  .92000 

13 

.14000 

13.33 

839.44 

.0014161 

870.7? 

.38166 

.08943 

1.92000 

13 

.14000 

16.67 

834,79 

.0014138 

872,18 

.37954 

.08988 

1.92000 

13 

. 140 00 

20.00 

830,35 

.0014115 

873.58 

.37753 

.09009 

1.92000 

13 

.14000 

23.33 

826.85 

.0014096 

874.76 

.37593 

.09009 

1.92000 

13 

.14000 

26.67 

824.45 

.0014084 

875.51 

.37484 

.09001 

1  .92000 

13 

.14000 

30.00 

623.73 

.0014079 

875.80 

.3745? 

.08998 

1.92000 

1* 

.15167 

0. 00 

772.49 

.0014541 

847.97 

.3512? 

.08750 

2.UM000 

16 

.15167 

3.33 

771.38 

.0014538 

848.19 

.35072 

.08767 

2.08000 

1* 

.15167 

6.67 

768,55 

.«0l4525 

848.91 

.34943 

.08813 

2.08000 

14 

.15167 

10.00 

764.35 

.0014506 

850.(6 

.34752 

.08871 

2 ,08000 

14 

.15167 

13.33 

759.01 

.0014482 

851.47 

.34509 

.08918 

2.08000 

14 

.15167 

16.67 

753.47 

.0014455 

853.05 

.34257 

.08943 

2.08000 

14 

.15167 

20.00 

747.99 

.0014430 

854.52 

.34308 

.08941 

2.08000 

14 

.15167 

23.33 

743.72 

.0614409 

855.79 

.33814 

. 68923 

2.0800ft 

85 
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COMPUTATION  OF  AERO-ACOUSTIC  PROPERTIES  OF  SUPPRESSOR  NOZZLES 


CASE 

NO.  1 

CRf) 

7-TURE  AR=?» 

3  NOZZLE 

-  V J=220O  EPS  - 

TT J=1 600  OEG-R 

AXIAL 

LOCATION 

=  1.1 

16667  (X/UEQ  =  15.99999) 

M 

R 

ANGLE 

U 

DENSITY 

TEMP. 

U/UREF 

TUPP.INT 

.  R/DEO 

14 

.15167 

?6.67 

740.75 

,00 1439S 

856.58 

.33679 

.08902 

2.08000 

14 

.15167 

30 . 00 

739.65 

.0014391 

856.84 

.33629 

.08893 

2.08000 

15 

.16313 

o.OC 

690.4? 

.0014898 

827.68 

.31391 

• GBS64 

2.24000 

15 

.16333 

3.00 

689.59 

.0014894 

827.88 

.31353 

.08575 

2.24000 

15 

.16333 

6.00 

686.80 

.C014885 

828.43 

.31226 

.08603 

2.24000 

15 

.16333 

9.00 

682.73 

.0014868 

829.35 

.3104) 

.08638 

2.24000 

15 

.16333 

12.00 

677.69 

.  o  0  14846 

83'j  .56 

.30812 

.08667 

2.24000 

15 

.16333 

15.00 

671.85 

.0 014822 

831.90 

.30546 

.08677 

2.24000 

15 

.  16333 

18.00 

665.83 

.0014799 

833.24 

.30272 

.08666 

2.24000 

15 

.16333 

?1 .00 

660.6? 

.0014775 

834.56 

.30035 

. 08640 

2.24000 

15 

.16333 

?4.oe 

666.27 

.0014757 

835.59 

.29838 

.08605 

2.24000 

15 

.16333 

27.  cc 

663.61 

.0014745 

836.28 

.2971? 

.08578 

2.24000 

15 

.16333 

30.00 

652.68 

.0014740 

836.57 

.29675 

.08568 

2.24000 

16 

.17500 

0.00 

606.81 

.0015289 

806.53 

.27589 

.08184 

2.40000 

16 

.17600 

3.0  3 

605.99 

.0015286 

806.75 

.27552 

.08193 

2.40000 

16 

.17500 

6.00 

602.90 

.0015276 

807.20 

.27411 

.08208 

2.40000 

16 

.17500 

9.00 

698.56 

.0015260 

808.06 

.27214 

.08227 

2.40000 

16 

.17500 

12.00 

592.66 

.0015240 

809.09 

.26946 

.08234 

2.40000 

16 

.17600 

15.00 

586.40 

.0015216 

810.36 

.26661 

.08225 

2.40000 

16 

.17500 

18.00 

579.83 

.0015192 

811,65 

.26362 

.08193 

2.40000 

16 

.17500 

21. OC 

574. CO 

.0015169 

812.90 

.26097 

.08149 

2.40000 

16 

.17500 

24.00 

569.29 

.0015150 

813.93 

.25683 

.08101 

2.40000 

16 

,17500 

27.00 

566.13 

.0015138 

814,56 

.25740 

.0806? 

2.40000 

16 

.17500 

30. CC 

565.13 

.0015133 

814.81 

,25b94 

.08049 

2.40000 

17 

,18667 

0.00 

524.23 

.0015714 

784,69 

.23834 

.07644 

2.56000 

17 

.18667 

2.73 

523.45 

.0015711 

784.83 

.23799 

.07647 

2.56000 

17 

.18667 

6.45 

520.84 

.0015705 

7R5.15 

.23680 

.07651 

2.56000 

17 

.18667 

8.18 

516.92 

.0015693 

785.74 

.23502 

.07655 

2.56000 

17 

.18667 

10.91 

511.71 

.0015678 

7R6.50 

.23266 

.07649 

2.56000 

17 

.18667 

1  3.64 

505.65 

.0015661 

787.38 

.22990 

.07627 

2.56000 

17 

.18667 

16.36 

499.45 

.0015640 

788,42 

.22708 

.07592 

2.56000 

17 

.18667 

19.09 

493.33 

.0015619 

789.45 

.22430 

.07540 

2.56000 

17 

.18667 

21 .8? 

487.93 

.00156)0 

790.4? 

.22184 

.07483 

2.56000 

17 

.18667 

24 . 66 

483.66 

.0015585 

791 ,?0 

.21990 

.07429 

2.56000 

17 

.18667 

27.27 

481 .02 

.0015574 

791 .74 

.21870 

.0739? 

2.56000 

17 

.18667 

30.00 

480.13 

.0015570 

791.94 

.21629 

.07378 

2.56000 

18 

.19833 

0.00 

444.85 

.0016173 

76?. 41 

.20226 

.06978 

2.72000 

18 

.19833 

2.73 

444.06 

.0016171 

762,53 

.20190 

.06978 

2.72000 

18 

.19833 

5.45 

441.44 

.0016166 

762.75 

.20070 

.06974 

2.72000 

18 

.19833 

8.18 

437.39 

.0016157 

763.17 

.19886 

.06964 

2.72000 

18 

.19833 

19.91 

43?. 15 

.0016145 

763.74 

.19648 

.0694? 

2.72000 

18 

.19833 

13.64 

426.10 

.0016130 

764.46 

.19373 

.06906 

2.72000 

18 

.19833 

16.36 

419.7? 

.0016113 

765.27 

.19083 

.06854 

2.72000 

18 

.19833 

19.  C-0 

413.45 

.0016095 

766.10 

.18798 

.06790 

2.72000 

18 

.19833 

21.82 

407.99 

.0016078 

766.94 

.18550 

.0672? 

2.72000 

18 

.19833 

24.65 

403.63 

.0016064 

767.6? 

.18352 

.06660 

2.72000 

18 

.19833 

27.27 

400.90 

.0016054 

768,08 

.18227 

.06616 

2.72000 

18 

.19833 

30.no 

399.90 

.0016051 

768. ?3 

.1818? 

.06600 

2.72000 

86 
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COMPUTATION  OF  Af PO-ACOUSTIC  PROPERTIES  OF  SUPPRESSOR  NOZZLES 


CASE 

NO,  1 

CRD 

7-TUfie  AR 

=2.3  MOZZLF 

-  V  J= 

2200  FPS 

-  TTJ= 1 600  DEG-R 

AXIAL 

LOCATION 

=  1. 

16667  (X/DEO  =  15. 

99999) 

M 

R 

ANGLE 

U 

DENSITY 

TFMP. 

U/UNEF 

TUPB.INT.  R/DEQ 

19 

.21000 

0.00 

370.73 

.001666? 

740.06 

.16855 

.06228 

2.88000 

19 

.21000 

2.50 

370.01 

.0016661 

74  o  .  1  0 

.16823 

.  C6225 

2.88000 

19 

.21000 

S.Of 

367.9? 

.0016658 

740.23 

.16728 

.06217 

2.88000 

19 

.21000 

7.50 

364.58 

.0016653 

740.45 

.16576 

.062G0 

2.88000 

19 

.21000 

10.  CC 

360.  1  1 

.0016647 

740,74 

.16373 

.0617? 

2.88000 

19 

.21000 

12.50 

354.93 

.0016637 

741.15 

.16137 

.06133 

2.88000 

19 

.21000 

15. OC 

349.35 

.0016625 

741 .68 

.15883 

.06081 

2.88000 

19 

.21000 

17.50 

343.58 

. CC 1 66 1 3 

742.23 

.15621 

.06017 

2.88000 

19 

.21000 

2C.00 

338,18 

.0016599 

742.84 

.15376 

.05948 

2.88000 

19 

.21000 

22.50 

333.44 

.0016587 

743.4-, 

.15160 

.05879 

2.88000 

19 

.21000 

2  S.Of 

329.79 

.0016576 

743.89 

. 14994 

.05821 

2.88000 

19 

.21000 

27.50 

327.48 

.0016569 

744.21 

.14889 

.05781 

2.88000 

19 

.21000 

30.00 

326.69 

.0016566 

744.32 

.14853 

.05767 

2.88000 

20 

.22167 

0.00 

303.31 

.0017173 

718.03 

.13790 

.05436 

3.04000 

20 

.22167 

2.50 

302.60 

.0017173 

718.02 

.13758 

.05431 

3.04000 

20 

.22167 

5. CO 

300.64 

.0017172 

718.06 

.13669 

.05418 

3.04000 

20 

.22167 

7.50 

297.45 

.0017171 

718.13 

. 1 3524 

.05394 

3.04000 

20 

.22167 

10.00 

293.36 

.0017167 

718,30 

.13338 

.05359 

3.04000 

20 

.22167 

12.50 

288.49 

.0017161 

718.5? 

.13116 

.05312 

3.04000 

20 

.22167 

15.00 

2«3.14 

.0017155 

718.8ft 

.12873 

.05252 

3.04000 

20 

.22167 

17.50 

277.74 

.0017145 

719.18 

.12628 

.05183 

3.04000 

20 

.  2  2 .1 6  7 

20.00 

272.60 

.0017136 

719.6,) 

.12394 

.05109 

3.04000 

20 

.22167 

22.50 

268.15 

.0017125 

720,03 

.12192 

.05038 

3.04000 

20 

.22167 

25. OC 

264.66 

.0017117 

7?C ,4ft 

.12033 

,04977 

3.04000 

20 

.221 67 

27.50 

262.47 

.0017110 

720,66 

.11933 

.04937 

3.04000 

20 

. 22 167 

30. OC 

261.72 

.0017108 

720,75 

.11899 

.0492? 

3.04000 

21 

.23333 

0.00 

243.54 

.0017699 

696,68 

.11073 

.04640 

3.20000 

21 

.23333 

2.31 

243.02 

.0017700 

696.67 

.11049 

.04635 

3.20000 

21 

.23333 

A. 6? 

241.46 

.0017701 

696,63 

.10978 

.04621 

3.20000 

21 

.23333 

6.92 

239.06 

.0017701 

696,62 

. 1C869 

.04598 

3.20000 

21 

.23333 

9.23 

235.80 

.0017701 

696,61 

.10721 

.04564 

3.20000 

21 

.23333 

11.54 

231.86 

.0017701 

696,6? 

.  1C54? 

.04520 

3.20000 

21 

.23333 

13.  85 

227.5? 

.0017699 

696,70 

.10344 

.04466 

3.20000 

21 

.23333 

16.15 

222.94 

.0017695 

696,83 

.10136 

,04404 

3.20000 

21 

.23333 

1 8.46 

218.46 

.0017690 

697,04 

.0993? 

.04336 

3.20000 

21 

.23333 

20.77 

214.3! 

.0017683 

697,31 

,09744 

.04268 

3.20000 

21 

.23333 

23.05 

210.74 

.0017676 

697,58 

.09582 

.04204 

3.20000 

21 

. 2333  3 

25.38 

?08.y0 

.0017670 

697.83 

.09457 

.04152 

3.20000 

21 

.23333 

27.69 

206.28 

.0017666 

698,00 

.09379 

.04118 

3.20000 

21 

.23333 

30.00 

205.69 

.0017664 

698.07 

.09352 

.04106 

3.20000 

22 

.24500 

0.00 

191.94 

.0018229 

676.44 

.08727 

.03874 

3.36000 

22 

•  24500 

2.14 

191.56 

.0018230 

676,41 

.08709 

,03869 

3.36000 

22 

.24500 

4.29 

190.41 

.0018231 

676.35 

.08657 

.03857 

3.36000 

22 

.24500 

6.43 

188.54 

.0018234 

676.25 

.08572 

.03835 

3.36000 

22 

.24500 

8.57 

186. 04 

.0018237 

676.15 

.08458 

.03804 

3.36000 

22 

.24500 

10.71 

183.00 

.0018240 

676.04 

.08320 

.03765 

3.360U0 

22 

.24500 

12.86 

179.58 

.'.'01824? 

675,97 

.08165 

.03718 

3.36000 

22 

.24500 

15.0ft 

175.91 

.001 824? 

675.94 

.07998 

.03664 

3.36000 

87 
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COMPUTATION  OF  AERO-ACOUSTIC  PROPERTIES  OF  SUPPRESSOR  NOZZLES 


CASE  NO,  1  CPO  7-TURF  A9=2.3  NOZ7LE  -  VJ=2200  FPS  -  TTJ=I600  OEG-R 


AXIAL  LOCATION  =  1.16667  (X/OEO  =  IS. 99999) 


M 

P 

ANGLE 

u 

DENSITY 

TEMP. 

U/UREF 

TUPB.INT 

.  R/DEU 

2? 

.24500 

17.14 

172.18 

.0018241 

675.97 

.07828 

.03605 

3.36000 

2? 

.24S00 

19.29 

168,57 

.0018239 

676.07 

.07664 

.03544 

3.36000 

2? 

.24500 

21.43 

165.28 

.0018235 

676.21 

.07515 

.03484 

3.36000 

2? 

.24500 

23.57 

162.49 

.0018231 

676.37 

.07388 

.03430 

3.36000 

2? 

.24500 

25.71 

l 60 . 36 

.0018226 

676.53 

.07291 

.03386 

3.36000 

2? 

.24500 

27.86 

159.03 

.0018223 

676.64 

.07231 

.03358 

3.36000 

22 

.24500 

3'J.OO 

158.58 

.0018222 

676.60 

.C7210 

.03348 

3.36000 

23 

.25667 

e.oo 

14R.57 

.0018751 

657.61 

.06755 

.03165 

3.52000 

23 

.25667 

2.14 

148.24 

.0018752 

657.58 

.06740 

.03160 

3.52000 

23 

.25667 

4.29 

147.26 

.0018754 

657.48 

.06695 

.03147 

3.52000 

23 

.25667 

6.43 

145.67 

. ‘010759 

657.34 

.C6623 

.03126 

3.52000 

23 

.25667 

0.57 

143.53 

.0018764 

657.16 

.06526 

.03096 

3.52000 

23 

.25667 

10.71 

140.94 

.0018769 

656.97 

.06408 

.03058 

3.52000 

23 

.25667 

12.86 

138. Cl 

.0018774 

656.79 

.06275 

.03013 

3.52000 

23 

.25667 

15.00 

134.07 

.0018778 

656.66 

.0613? 

.0296? 

3.52000 

23 

.25667 

17.14 

131.67 

.0018780 

656.58 

.05987 

.02906 

3.52000 

23 

.25667 

19.29 

128.57 

.0018780 

656.57 

.05«46 

.02849 

3.52000 

23 

.25667 

21.43 

125.74 

.0018779 

656.6? 

.05717 

.02794 

3.52000 

23 

.25667 

23.57 

123.33 

.0018776 

656.71 

.05607 

.02744 

3.52000 

23 

.25667 

25.71 

121.49 

.0018774 

656.81 

.05524 

.02704 

3.52000 

23 

.25667 

27.86 

120.34 

.0018771 

656.89 

.05471 

.02678 

3.52000 

23 

.25667 

30.00 

119.95 

.0018771 

656.92 

.05453 

.02669 

3.52000 

24 

.26033 

0.00 

112.88 

.0019255 

640.38 

.05132 

.02529 

3.68000 

24 

.  26833 

2.00 

112.64 

.0019256 

640.34 

.05121 

.02526 

3.68000 

24 

•  26533 

4.00 

111.93 

.0019260 

640.24 

.05089 

.02515 

3.68000 

24 

.26033 

6.0  C 

110.78 

•C019264 

640.08 

.05037 

.02497 

3.68000 

24 

.26833 

8.0.1 

109.2? 

.0019271 

639.87 

.04966 

.02473 

3.68000 

24 

.26833 

K'.OO 

107.32 

.3019278 

639.64 

.04879 

.0244? 

3.68000 

24 

.26033 

12.00 

105.13 

.0019285 

639.41 

.04780 

.02406 

3.68000 

24 

.26833 

14.00 

102.75 

.0019291 

639.19 

.04672 

.02364 

3.68000 

24 

.26833 

16.00 

100.27 

.0019296 

639.0? 

.04559 

.02319 

3.68000 

24 

.26833 

18. OC 

97.79 

.0019300 

638.90 

.04446 

.02272 

3.60000 

24 

.26033 

20.00 

95.44 

.0019302 

638,84 

.04339 

.02224 

3.68000 

24 

.26033 

22.00 

93.31 

.0019302 

638.83 

.04242 

.02179 

3.60000 

24 

,26833 

24.00 

91.51 

.0019301 

638.87 

.04161 

.02139 

3.68000 

24 

.  26833 

26.00 

90.16 

.0019299 

630.92 

.04099 

.02108 

3.68000 

24 

.26833 

28.00 

89.31 

.0019298 

638.97 

.04060 

.02088 

3.60000 

24 

.26833 

10.05 

89.0? 

.0019297 

638.99 

.04047 

.02081 

3.68000 

2S 

.28000 

c.oc 

84.19 

.0019734 

624.84 

.03828 

.01978 

3.84000 

25 

.28000 

2.00 

84.00 

.0019735 

624.80 

.03819 

.01975 

3.84000 

25 

.28000 

4. 00 

83.43 

.0019739 

624.69 

.03793 

.01966 

3.04000 

25 

.28000 

6.00 

82.50 

.0019745 

624.51 

.03751 

.01950 

3.04000 

25 

.28000 

0.00 

«1.?4 

.0019752 

624,28 

.03694 

.01928 

3.H40O0 

25 

.28000 

l-’.oc 

79.70 

.0019760 

624.01 

.03624 

.01901 

3.84000 

25 

.28000 

12.02 

77.94 

.1-019769 

623.74 

.03544 

.01868 

3.84000 

25 

.20000 

14. OC 

76. '.I 

•  00 1 9777 

623. 4R 

.03456 

.0183? 

3.84000 

25 

.28000 

16. or 

74.00 

.0019785 

623.25 

.03365 

.0179? 

3.04000 

25 

.28000 

1 8. mi 

71  .99 

.0019790 

623.00 

.03273 

.01751 

3.04000 

88 
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COMPUTATION  OF  AERO-ACOUST IC  PROPERTIES  OF  SUPPRESSOR  NOZZLES 


CASE  NO,  1  CPD  7-TURe  AP=2,3  NOZZLE  -  VJ=?200  FPS  -  TTJ=1600  DEG-R 


AXIAL  LOCATION  =  l. lfe*(S7  <X/OEO  =  15.99990) 


M 

P 

angle 

U 

DENSITY 

TEMP. 

U/UREF 

TURR.lNT 

.  R/DEl) 

?s 

.28000 

20.00 

70.07 

.0019794 

622.97 

.03186 

.01710 

3.84000 

25 

.28000 

22.00 

68.33 

.0019795 

622.91 

.03107 

.01671 

3.b4000 

25 

• 280OC 

28.00 

66.87 

.0019795 

622.91 

.03040 

.01637 

3.84000 

25 

.28000 

26.00 

65.76 

.0019795 

622.93 

.02990 

.01611 

3.84000 

25 

.28000 

28.0  0 

6S./6 

.0019794 

622.96 

.02958 

.01593 

3.84000 

25 

•  28C00 

3C.0C 

64.63 

.  J019793 

622.97 

.02947 

.01587 

3.84000 

26 

.291*7 

v'.OP 

61.63 

.ii  020180 

611.04 

.02802 

.01514 

4.00000 

26 

.29167 

1.87 

61.50 

.0020181 

611.00 

.02796 

.01512 

4.00000 

26 

.291*7 

3.75 

61.11 

.0020185 

610.90 

.02778 

.01505 

4.00000 

26 

.291*7 

5.62 

60.46 

,n020190 

610.73 

.02749 

.01493 

4.00000 

26 

.291*7 

7.5.: 

59.59 

.',•020197 

610.52 

.02709 

.01476 

4.00000 

26 

.29167 

9.37 

58.51 

.0020206 

61 1.26 

.02660 

.01455 

4.00000 

2* 

.291*7 

11.25 

57.26 

.0020215 

609.99 

.02603 

.01431 

4.00000 

26 

.29167 

13.12 

65.88 

. J020224 

609.71 

.02540 

.01403 

4.00000 

26 

.29167 

15.00 

54.41 

.0020232 

609.46 

.02474 

.01373 

4.00000 

26 

.29167 

16.87 

52.91 

.0020240 

6C9.24 

.0240* 

.01341 

4.00000 

26 

.291*7 

18.75 

51.44 

.0020245 

609.06 

.02339 

.01308 

4.00000 

26 

.291*7 

2C.62 

50.05 

.0020249 

608.95 

.02276 

.01277 

4.00000 

26 

.291*7 

22.50 

48.81 

.0020252 

608,88 

.02219 

.01247 

4.00000 

26 

.29167 

28.37 

47.  *6 

.0020252 

698.86 

.0217? 

.01222 

4.00000 

26 

.291*7 

26.25 

46.98 

.(  020252 

608.86 

.02136 

.01202 

4.00000 

26 

.29167 

26.12 

46.49 

.C 020252 

608,87 

.02114 

.01189 

4.00000 

26 

.29167 

70.0'j 

46.32 

."020252 

608,88 

.02106 

.01185 

4.00000 

27 

.30333 

0.00 

44.28 

.0020586 

598.94 

.02013 

.01135 

4.16000 

27 

.70333 

1  .87 

44.18 

.0020589 

598.90 

.02009 

.01133 

4.16000 

27 

.30333 

3.75 

43.88 

.0020592 

598.80 

.01995 

.01127 

4.16000 

27 

.30333 

5.*2 

43.39 

.0020598 

598.63 

.01973 

.01117 

4.16000 

27 

.30333 

7.50 

42.7  2 

.0020606 

598,42 

.01942 

.01103 

4.16000 

27 

.30333 

9.37 

41.89 

.0020614 

598,16 

.01904 

.01086 

4.16000 

27 

.30333 

11.25 

40.93 

.0020624 

697,88 

.01861 

.010*6 

4.16000 

27 

.30333 

13.12 

39.87 

.0020634 

597,60 

.01813 

.01043 

4.16000 

27 

.30333 

16.00 

38.74 

.0020643 

597.34 

.01761 

.01018 

4.16000 

27 

.30333 

16.87 

37.58 

.0020651 

597,10 

.01709 

.00992 

4.16000 

27 

.30333 

18.75 

36.45 

.0020658 

596.91 

.01657 

.00966 

4.16000 

27 

.30333 

20.62 

35.37 

.002066? 

596,77 

.01608 

.00940 

4.16000 

27 

.30333 

22.50 

34.41 

.0020665 

596,69 

.01564 

.00916 

4.16000 

27 

.30333 

24.37 

33.60 

.0020667 

596,65 

.01528 

.00895 

4.16000 

27 

.30333 

26.25 

32.99 

.0020667 

596,63 

.01500 

.00879 

4.16000 

27 

.30333 

26.12 

32.61 

.0020667 

596.64 

.01482 

.00869 

4.16000 

27 

.30333 

30.00 

32.48 

.0020667 

596.64 

.01477 

.00865 

4.16000 

28 

.31500 

O.OC 

31.22 

.0020955 

588.45 

.01419 

.0083? 

4.32000 

28 

.31500 

1.76 

31.15 

.0020956 

588.42 

.01416 

.00831 

4.32000 

28 

.71500 

3.53 

30.95 

.0020959 

588,34 

.01407 

.00827 

4.32000 

28 

.71500 

5.29 

30.63 

.0020964 

588,19 

.01392 

.00819 

4.32000 

28 

.31500 

7.06 

30.18 

.0020970 

588.01 

.01372 

.00810 

4.32000 

28 

.31500 

8.82 

29.62 

.0020978 

587.78 

.01347 

.00798 

4.32000 

28 

.31500 

10.59 

28.97 

.0020987 

587,54 

.01317 

.00783 

4.32000 

28 

.71560 

12.75 

28.25 

.0020996 

5P7.2R 

.01284 

.00767 

4.32000 

89 
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APPENDIX  D 


PROGRAM  SOURCE  CODE  LISTING 


This  appendix  contains  the  FORTRAN  IV  source  code  listing  for  the 
aeroacoustic  prediction  model,  suitable  for  running  on  the  CDC  7600  computer. 
The  listing  of  subroutines  is  in  alphabetical  order,  as  follows: 


1.  MAIN  Program  (MGB) 

2 .  ARRCCOS 

3.  ATMOS 

4.  CRD 

5.  ERF 

6.  LSPFIT 

7 .  OUTPUT 

9.  PNLC 

9.  SHOCK 

10.  SLICE 

11 .  TPNLC 
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PPpGPAP  MGU  (  iNPGT.OUTPUT.TAPES^NPUT.TAPEf^OIJTPUT  ) 
M*r,*R  MANl*GLIt«F*PALSfl  JET  NOISE  PREDICTION 


DATA  EO/SO  10  0.1  ?S  ,  1  bO  ,  ?0  0 . 2S0 , 3  1 S .  4  0  0  .  SOlf  .  0 

i  i?sc . i6oo»?ooo.?soc. :iisv,4;  ao.5()co.e>io(;.pooo.  lyocu 
J.PSOOP,  31  so  a.4"-)^  O.S  A  000, 030  00. 80000*1  oocoo/ 
DATA  ALPHT/1Q*0.S/.  pKTAIN/2<.*4. j/t  DEL T I N/24 *4 .  ■  / 


FF  <K  )  =.i,o 

IF (flCH?.EO.0.C)  GO  TO  POO 
UE<K)=SQPT(ACH?°SV2) 

FF (K) =PU?F (k )*CP» (TT (K ) -TT ( I ) ) /UF (K) 
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NPO=NP»INT 

IF (NPPlNT.tr. 0)  WPI Tf 


'^fXjO^kT'Cf'-XCT'O— «  fO  if  .£>  r-  X(TO  —  «f\jm  f  lP  X  ^  X  aOMOjnji/^slN  XT' 

XXX  X  r  X  X  X  X  \T  C?  .-7*  777777700c  -i>  3COOCO  —  f\, 


x  x  x  x  axxxxxxa:  xxxxxxxctx  xxxxxxxxxxxxxxxxxxxxxxxxxxxx 
oocc  ccoccccc  ococccc  cccwowcc  cccc  ec/ccjcceo^cccoceoce. 

5221X25222  2  2  2  X22i22iXXX5Xii2.XX5.2iX2i2LX2XX2XXXX5.X5.  2 
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CARCCOS  ARC  COSINE  (PRINCIPAL  VALUE)  ARCCOS 

FUNCTION  APCCOS(X)  ARCCOS 

IF (X.GT.O.O)  GO  TO  S  ARCCOS 

IF(X.LT.O.O)  GO  TO  10  ARCCOS 
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